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SOURCES AND SINKS AT THE AXIS OF A 
ROTATING LIQUIDt 


By ROBERT R. LONG 
(The Johns Hopkins University, Baltimore, Maryland) 


Received 16 June 1955; revise received 27 March 1956] 


SUMMARY 

A solution is obtained for the flow of a rotating, frictionless, incompressible fluid 
due to a strong source or sink at the axis of rotation. The type of motion is con- 
trolled by the value of the Rossby number, Ro, a ratio of inertial and Coriolis forces. 
The solution resembles potential flow if Ro is very large. As Ro decreases the sink 
draws more and more from the axis of rotation until, at Ro 0-261, the fluid 
approaches the sink in a jet centred at the axis and with a radius about half that 
of the cylinder. No solution is obtained for smaller values of Ro. Ifa jet is postulated 
for small values of Ro, it is shown that its radius decreases as the cube root of the 
flux. Several experimental photographs are shown. They contain some of the 
featuresof thetheory. The jet type of motion becomes very pronounced for weak sinks. 
1. Theory 
[x a previous paper the author (1) has shown that a solution of the differen- 
tial equation 


cts = ots l cud 
J | Y ° 19 9 
ous sO Ug p~ ( ] ) 


> 


cx Cp* pp ep 

represents steady-state, axially symmetric flow of a rotating, frictionless, 
incompressible fluid; y% is the Stokes’s stream-function. The coordinate 
system is shown in Fig. 1. In the derivation of equation (1) it was assumed 
that at x © the fluid has a constant absolute angular velocity Q about 
the x-axis and a constant linear velocity u, parallel to the axis of rotation. 
The parameter o is 2Q up. 
The purpose of this note is to apply (1) to the problem of sources and 
sinks located on the axis of rotation. We suppose first that there is a sink 
at the point « = 0, p = Oin an infinitely long cylinder of radius b containing 
rotating liquid. We require the streamline pattern to be symmetric about 


the plane x = 0 and assume that at |x © the absolute velocity field is 


a constant angular velocity Q and constant velocity u, toward the sink. 
In view of the symmetry, the problem is equivalent to flow from the 
negative x-axis toward a hole at x = 0, p = 0 in a wall coincident with 
the plane x = 0. An appropriate solution of (1) is 
Uy p* 2_g2tr 
ub oF 4 p> Ave a “© J, (xp). 


» 
- x 


+ This research was sponsored by the Office of Naval Research under Contract 
N-onr-248(31). 
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386 ROBERT R. LONG 


Unless « = o the solution (2) yields uniform flow, uw), as x + —o, ff 
was also shown by the author (2) that w, the component of the absolute 
velocity tangent to a circle centred at the axis of rotation, is related to 


by the equation 
7 wp ous. (3) 


Equations (2) and (3) show that w—> Qp as x -> —o, or solid rotation. 














fr 


Fic. 1. Coordinate system. The x-axis is taken along the axis of rotation. 


y 


At the wall of the cylinder, the uniform flux through a plane perpen- 


dicular to the axis is —2zy. Hence at p = b, 4 = —u,6?/2. Therefore, 
ab = z,, where z, are the zeros of the Bessel function J,(z), the first two 
being z, = 3-8317, z. = 7:0156. The solution (2) becomes 
Up p> 

oP” % 2 Ro*)(zlb) J [> P 
ab — 4 p > A, pez —Ro “)2(x, SJ; (205 : (4) 

~- ) 
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We have substituted for (cb)~! the symbol Ro Up/2Qb. This is called 
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the Rossby number in meteorological literature (1). It is a measure of the 
importance of the rotational or Coriolis forces in rotating systems. The 
summation in (4) begins at »,, where n, is the first integer for which 
Ro > 1/z,. For smaller integers we would have oscillating functions of x 
in (4) and disturbed motion at infinity. 

We notice that & 0 on the axis p 0. We obtain a sink at x 0, 


0 by requiring that & has the same value at the plane x = 0 (except 


it the point x = 0, p = 0) as it has on the cylinder wall, namely — uw, b?/2. 
The resulting jump discontinuity of % yields the infinite velocity at (0, 0) 
required in a sink. Hence for zx 0.0 p oa 2, 
U, b? p~ ay p - 
( 1] p A,,J,| Zn}: (5) 
») 
» }? , b 


Multiplying both sides of (5) by J,(z,, p/b) and integrating with respect to 
o'b from 0 to 1, we find 


‘i. a (6) 


The solution becomes 


UU, p= - exp! (z2 — Ro-?)!a/b = 
’ of u, bp % p| n | BAC p/b). (7) 
» / y JS2(z ) n 
_ oa “Hn O\Wn 
Since the summation begins at n 1, the solution exists only if Ro exceeds 
l/z,, or Re 0-261. This corresponds to low or moderate values of the 


angular velocity of the system or to strong sinks. Obviously we may change 


the sign of w, and (7) will still satisfy the differential equation (1). With 
the proper choice of pressure, the resulting motion will then satisfy the 
primitive equations of motion and continuity for the steady symmetric 
flow of a frictionless, incompressible, rotating fluid. The pressure may be 


found from the equation 


Pp, w+e+u* — px(po) , Mo , 22" po (8) 
q 2 q 2 
where q is the constant density, p,, is the pressure at x —oo, and pz, is 


the Lagrangian distance of a streamline from the axis. Equation (7), with 
a negative value of uw», yields flow due to the emission of fluid from the 
point (0,0). The emitted fluid possesses vorticity so that this system may 
be regarded as a source combined with a vortex of a special type. The case 
of a source emitting irrotational fluid has been considered by Barua (3). 
In the latter problem the irrotational fluid moves along the axis in the 
form of a cylindrical jet. 

The flow patterns given by (7) are shown for several values of Ro in Fig. 2. 
We note that as Ro decreases from oo the sink draws more and more from 
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Fic. 2. Flow patterns computed from equation (7). The curves are isolines 


of po/b = (—2yb/u, b?)* (see Long, 1953). 





Ro:@ 


Ro=.286 


Ro=.261 


— —~w 


i it 


en er 





sO 


the a: 
Ro = 


LZ>- 


This | 
eylin 
the d 

If 


we fi 


wher 


For ¢ 


along 


whet 
wave 
wave 
these 
sink 
no fi 
a lor 
solut 
the 1 
sup} 
sink: 

Al 
0:26 
syst 
weal 
cirel 
patt 
of ni 


Whe 











SOURCES AND SINKS AT THE AXIS OF A ROTATING LIQUID 389 
the axis of rotation. This tendency is a maximum for the critical value of 
Ro = 0-261, in which case the fluid approaches the sink in a ‘jet’. As 


>? 


© the velocity tends to 


Jo(Z, p |. (9) 


J5(2) 





This is a maximum in the vicinity of the axis, becoming negative near the 
evlinder wall. Although the condition that u— u, at x % is violated, 
the difficulty is avoided by choosing Ro slightly larger than 0-261. 

If we solve equation (1) for free stationary waves in an infinite cylinder 


we find, typically, 


) 


us “oP pB ei Ale, 5) (10) 
2 7 ‘bh 
o2\ J 
where k (o° i) (11) 
b, 

For a given value of m we may interpret (1C) and (11) as a wave moving 

long the axis of a cylinder at a speed ¢ u, relative to the liquid. Then 

$()°)2 (12) 

2° + 4ar*b?/A2 

where A, is the wavelength. The speed decreases with n and increases with 

wavelength. The maximum speed of any possible wave of this type is a 

wave of infinite length with no internal nodal surfaces, i.e. z,, = z,. Applying 

these results to the problem of the sink, if ob < z,, the fluid approaches the 

sink at great distances with a speed greater than that of any wave. Hence 

no free wave can advance against this current. If ob exceeds z, slightly, 

along wave can move upstream against the current. The steady-state 

solution at ob = z, shows that this has happened and the wave has altered 

the upstream velocity distribution to that of equation (9). This reasoning 

supports the assumption that the flow at x © is uniform for strong 
sinks. 


Although a steady-state solution exists for all values of Ro greater than 
261, we find that closed circulations begin to form in some parts of the 
system when Ro is close to the critical value. These circulations remain 
weak (Fig. 2) even at Ro (-286. Very near Ro = 0-261 the closed 
circulation strengthens, spreads upstream, and develops finally into the 
pattern of Ro = 0-261 shown in Fig. 2. The significance of the appearance 
of negative values of u follows from the equation 


| cus l C 
u (wp). (13) 


pcp op € p 


When wu is negative the absolute angular momentum decreases outward 
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locally, and we would expect this to be unstable. This suggests the possi- 
bility that weak sinks may lead to flows that are essentially turbulent 
and unsteady, and that no physically meaningful steady-state solution 
may exist. Since we may write 

FP 


Ro Yemky (14) 


where F is the volume flux through the sink, we are Jed to the further 
possibility that no steady-state flow exists for a sink of any finite strength 
in a completely unbounded rotating fluid (b 1). In this case, however. 
we may postulate the existence of a quasi-steady jet of radius b, moving 
toward the sink. If so, the only pertinent quantities, F, Q, bo, lead to 


only one non-dimensional number, say 
Ro’ 3 (15) 


which must therefore be an absolute constant. The magnitude of the 
constant may be inferred from a consideration of the jet existing theoreti- 
cally at Ro = 0-261. Taking the boundary of the jet at the value of p 
where |\du/dp| is a maximum, b,/b ~ 0-48 and 

3 

Ro’ Ro”. ~ 2:3. (16) 

a7) 
This is probably a crude estimate since the theoretical jet has a radius that 
is not small compared to that of the cylinder and therefore does not resemble 
a jet in an unlimited fluid.t 


2. Experimental results 

An experimental study was made of the effect of a sink at the axis of 
rotation of a tall cylinder. A cylinder of water of radius ) = 5-7 em. and 
height 125 em. is positioned on a rotating turntable and brought up to 
solid rotation. Water is then extracted at the bottom at the axis of rotation 
through a hole of radius 0-66 em. The flux and rotation rate may be varied 
to obtain different values of Ro = F/2Q7b*. The flow pattern is recorded 
by streak photography, using a camera rotating with the turntable, with 
line of sight perpendicular to a plane of light passed through the axis of 
rotation. Aluminium particles are used as tracers. 

The experimental photographs are shown in Figs. 3—7 in order of decreas- 
ing Rossby number. Figs. 3 and 4, with Ro values of 0 and 0-281 respec- 
tively, show the effect of a moderate amount of rotation. The sink draws 
more strongly from regions near the axis of rotation in Fig. 4, and there is 


| Sir Geoffrey Taylor, in a private communication to the author, remarked that Barua 
obtained Ro’ 0-31 for a jet due to a source at the axis of a rotating fluid. 
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a definite irdication of the eddies near the walls, predicted by the theory 
and shown in the 0-286 drawing of Fig. 2. No abrupt change of flow pattern 
at the critical value of the Rossby number. 


in the experiments 


occurs 





Fic. 3. Flow toward a sink. No rotation, Ro 0. 





Fic. 4. Flow toward a sink. Ro 0-281. 





Fic. 5. Flow toward a sink. Ro 0-203. 


Instead, the tendency for a jet formation seems to increase in a continuous 
manner. This is not surprising in view of the fact that the experiment is 


only quasi-steady. This point is illustrated in Fig. 5. The Rossby number, 
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0-203, is well below the critical value. A well-defined jet is being formed 
in the vicinity of the sink, apparently by upstream propagation of the | 1. R- I 








waves seen in the photograph. These waves have progressed a short distance | J. 
OS saci 
To 
3, 8. N 
} 9 
‘ 
Fic. 6. Flow toward a sink. Ro = 0-100. 
> 
i 
Fic. 7. Flow toward a sink. Ro = 0-006. 
ahead of the sink and the approaching flow is only affected in this vicinity. 
For smaller Rossby numbers the formation of the jet is more rapid. In 
Fig. 6 it exists far ahead of the sink. 
~ ‘ 


Fig. 7 illustrates the flow at an extremely low value of Ro = 0-V061, 
The jet is concentrated very near the axis of rotation and is spinning very 
rapidly. A single tracer particle has performed five revolutions about the 
axis in the } sec. time exposure of the photograph, yielding an angular 
velocity about 60 times the basic rotation (27 r.p.m.). 

An effort may be made to recompute the constant in equation (16) from 4 
experimental observation. Thus in Fig. 6, if we take the diameter of the 
jet as, perhaps, } the diameter of the vessel, equation (16) yields Ro’ ~ 1-6. 
Similar calculations for experiments with lower Rossby numbers also lead 
to values around 1-0. In all cases, however, the diameter of the jet is 
difficult to estimate in the photographs and Ro’ is very sensitive to varia- 


tions in this estimate. 
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THE SHAPE OF THE NAPPE OF A THIN WATERFALL 


By G. N. LANCE and E. C. DELAND 


(Department of Mathematics, The University, Southampton, and 
Department of Engineering, University of California, Los Angeles) 


[Received 27 October 1955] 


SUMMARY 
The shape of the nappe of a thin waterfall is obtained by solving the non-linear 
differential equation which describes the shape, on a mechanical differential analyser, | 
All forces acting on the waterfall, namely, gravity, surface tension, and air-pressure, 
are incorporated in the analysis. Previous writers have considered only the force of 
gravity. It is shown that the shape is extremely sensitive to air-pressure. Through- 
out, a comparison is made with the water-bell phenomenon which has been discussed 


previously and which is the three-dimensional analogue of the present work. 


1. Introduction 

Equations have been derived giving the shape of the nappe of a waterfall 
by several authors. R. V. Southwell and G. Vaisey (1) obtained the shape 
of a thick fall by relaxation methods. F. W. Blaisdell (2) and A. Fathy 
and M. 8. Amin (3) assumed the fall to be thin, an assumption that is also 
made here; however, all these authors assumed that the surface tension 
and air-pressure forces were negligible. Nevertheless, it has been shown 
in connexion with the water-bell problem by G. N. Lance and R. L. Perry (4) 
and by G. N. Lance and E. C. Deland (5) that air-pressure plays an impor- 
tant role in determining the shape of the nappe. 

The present paper gives the derivation of the differential equation whose 
solution gives the shape of the nappe of a waterfall. The depth of water 
flowing over is assumed to be small so that the velocity of all particles is the 
same. The direction of projection need not be horizontal. The differential 
equation, which is non-linear, contains three parameters, one involves 7' the 
surface tension and one involves the air-pressure difference II — p, where I] 
is the air-pressure on the upper nappe and p that on the lower nappe. Such 
a pressure difference could be produced either by air trapped below the 
nappe or by a wind blowing across the face of the fall. 

A first integral of the second order differential equation can readily be 
obtained but a complete solution involves incomplete elliptic integrals of 
the third kind. It was decided, therefore, to solve the equation on the 
differential analyser at the University of California, Los Angeles. This 
enables us to obtain curves for a wide variety of parameters. 


(Quart. Journ. Mech. and Applied Math., Vol. IX, Pt. 4 (1956)] 
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SHAPE OF THE NAPPE OF A THIN WATERFALL 395 
2. Derivation of the differential equation of the nappe 
In Fig. 1 the origin of coordinates is the point from which the water is 
projected, t he x-axis is horizontal, x being positive in the direction of flow, 
the y coordinate is measured positive downwards. 


O ; 7 _ x 





Outside, Pressures [I 


inside, Pressure=P 


YY 
1. The coordinate system. 
Let Y be the mass of water flowing per unit 'ength of fall per unit time; 
then if p is the density of water, v is the velocity of the water, and h is the 


thickness of the fall at a general point of its path, the equation of continuity 


may be written 


a phv. (1) 
The tai sential equation of motion is 
hos og sin us h os pudv ds, (2) 
where ds is the are length, ys is the slope of the curve, and g is the acceleration 
due to gravity. The normal equation of motion may be written 
(U1 ] 2T/R) phy cos us v oh R, (3) 
where 7’ is the surface tension of water, FP is the radius of curvature, and I] 
and he air-pressures on upper and lower nappes respectively. 
Now can be integrated and gives the energy equation 
ye ve 2gy, (4) 
where the velocity of projection. Also, using (1) to eliminate hp from 
Oo), We et ion 
. 2T di (11 p) g COS os dis (5) 
v n o 
V / V - ds x 
and’ may be eliminated from (5). We introduce the following three para- 


meters 


(Il—p)/Qvp. B = 2T/Qvy, and y = g/t. 
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The equation which results, namely 
y"[(1+2yy)!—B] = a(1+y’?)!+-y(1+-y)(1+-2yy)+4, (6) 
is the differential equation of the nappe, the prime denotes differentiation 
with respect to x. Since the thickness h is small the solution of (6) gives 
the shape of either nappe. 
A first integral of this equation can be obtained since x does not appear 
explicitly in the equation; we find 
(1+y’?)-![(1+2yy)!—B] = const—a(1+2yy)/2y. 
The constant is chosen so that y’ = tana, when y = 0. 
In the special case of horizontal projection we let 4, = 0. Hence 


(1+ y)*[(1+ 2yy)!—B] = (1—B)cos yYyo—ary. (7) 
The expression for x may be written down 
uv 
{ 1+2yt)!—B \? 4 
x ( ( <a -1| “dt. (8) 
7 (1—B)cos ys, —at 
0 


However, it can be shown that the integral on the right-hand side of (8) in- 
volves elliptic integrals of the third kind. For this reason we do not pursue 
this line of attack further. 


3. Discussion of the first integral 
Several interesting results may be obtained directly from (7). The curve 
has maxima or minima when y’ = 0. The values of y at which y’ = 0 are 


iven by 7 , 
oe (1+2yy)! = +{(1—B)cosy)+B—ay}. (9) 


Since v = v)(1-+-2yy)! is necessarily positive, it is required to find the 
roots of (9) which make the right-hand side positive. In the special case of 
horizontal projection (4, = 0) these are 

y 0 and | x +y— 2aB+-(a?+-y?-4 Z2ay— 4aBy)* | x2, (10) 

The water will be falling (or rising) vertically when y’ is infinite, this 

occurs when ay = (1—f)cos¢, in the general case and when 
XY 1—p (11) 
in the case of horizontal projection. 

The integral curve will contain a singular point when y’ is indeterminate, 
which occurs when (10) and (11) are satisfied simultaneously. The value 
of « for such a curve is given by 


(1—B)/« [ ow +-y—2aB+ (a?+-y?-+ 2ay- 4aBy)* | x? (12) 
or ¥ 2y/(1-+B). (13) 
The y-coordinate of the singular point is —(1—f?)/2y. It can easily be 


shown that near this point the shape of the curve is approximated to by 
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SHAPE OF THE 


the vertex of a semi-cubical parabola whose tangent is vertical. Since 


(6) 8 <l and} 0 the cusp always occurs above the horizontal; also for a 
tion given 8 and y there is just one value of «, namely, that given by (13), which 
rives gives a curve with a cusp. 

Thus the present case is easier to solve than the three-dimensional one 

Dear (4and 5) because in the latter case no first integral can be found and hence 
the value of « which gives the singular curve is unknown. 

4. Differential analyser solution of equations 
Since a complete solution cannot be obtained in terms of elementary 
functions it was decided to use the differential analyser for this purpose. 

i The values of 8 and y are obtained from reasonable assumptions about the 

magnitude of Q and vy. Then solutions were obtained for varying values 
| ofa. The curves of particular interest are those for which « takes the 

8 particular values 

(i) « = 0: natural waterfall ; 

in- ? ii) a 2y/(1+ 8): cuspidal curve. 

” Some typical curves are shown as Figs. 2, 3, and 4. The curves of Fig. 2 
are obtained using Q 13-26 g./sec./unit length and v, = 90 cm./sec., 
with the result that 8 = j 0-12. The values of « chosen lie between 

0-2 and 0-4, with the critical value corresponding to (ii) above, i.e. 

()-216, giving the cuspidal solution. 

“55 The curves of Fig. 3 differ from those of Fig. 2 in that the velocity of 

9 projection is v, = 30 em./sec., so that B = 0-362 and y 1-09. In the 

he present case — 2-0 x 0-4 and the cusp occurs for « —1-60. It can 

of easily be verified that the maxima and minima of the curves occur at the 
levels given by (10); in particular they all touch the line y = 0. It should 

0) be noted that the integral curves are oscillatory in the x-direction and that 
the period depends on the value of «. When « < 0 the period decreases as 

"= increases. This periodicity could also have been foreseen from (8) since 
the elliptic functions are periodic. 

Fig. 4, in which the parameters have the same values as in Fig. 2, shows 
the effects of projection at 1S. 

" 5. Conclusions 

One important difference between the water-bells (5) and the present 
problem is the fact that in the former the ‘loops’ occurred below the axis, 
| whereas in the latter all the looped curves are above and bounded by the 

axis. Hence it would be surprising if the unstable cusp-shaped curves 

S actually appear on a waterfall, even if the pressure difference could be made 


great enough 
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Previously, it was impossible to derive theoretically the value of « 
required to give the ‘singular’ curve; it is clear, however, that such a curve 
would intersect the water surface. In the falls case the singular curve does 


not intersect the water sul face. 
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SUMMARY 


This paper deals with a two-dimensional aerofoil moving in an inviscid, incom. | 


pressible fluid bounded by a plane wall. The solution is obtained by extending thy 


method of Green (1) developed for a stationary aerofoil in a uniform stream near 


a plane boundary, in which the results are given in series form. The results given 
apply to a fairly general type of aerofoil, but special attention is paid to those of 
the flat plate, symmetrical, and circular are types in order to assess the effects of 
thickness and camber. 


1. Introduction 

THE problem of finding the lift and moment forces acting on a stationary 
two-dimensional aerofoil in a uniform stream of inviscid incompressible 
fluid bounded by a plane wall has been treated quite extensively by several 
authors (Green (1, 2), Tomotika and co-authors (3), Havelock (4), and 
Fujikawa (5)), but the case of an aerofoil moving in a general manner in 
such a bounded region has not been discussed. The most suitable method 
of solution was found to be an extension of that developed by Green (1) for 
the stationary aerofoil in a uniform stream, in which results are given in 
series form. 


A fairly general type of aerofoil is considered, and the complex potential 


function appropriate for this is determined, allowance being made for a 


circulation around the aerofoil. The pressure distribution throughout the 
fluid and on the aerofoil boundary may be determined by using a modifi- 
cation of the usual pressure equation, and this enables the forces and couple 
acting to be calculated by using the ordinary contour integrals containing 
the pressure. The results thus obtained may be expressed in series form, 
and this is done for three different types of aerofoil. 


2. The conditions governing the complex potential 2 

Consider the two-dimensional motion due to a cylinder of aerofoil cross- 
section moving in a general manner in the neighbourhood of a plane 
boundary extending to infinity, the fluid being assumed at rest at infinity 
where the presence of the aerofoil is unnoticed. Since the motion is purely 


two-dimensional, we may simply consider the flow in a plane defined by ' 


a normal cross-section of the cylinder. 


[Quart. Journ. Mech, and Applied Math., Vol. IX, Pt. 4 (1956) 
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To specify the motion, we define two sets of moving rectangular Cartesian 


axes, the (x, y)-axes fixed in the aerofoil and moving with it in all respects, 


and the _y')-axes whose origin is also fixed inside the aerofoil section 


but which do not rotate with the aerofoil and so remain at all times parallel 
to a set of axes fixed in space. The orientation of these (x’, y’)-axes is so 
chosen that their real axis and the boundary wall are parallel. This enables 


us to write the boundary wall as y’ b at any time (b being a function 

of the time 7), and the force components perpendicular to and parallel to 
the wall as Y’ and X’ respectively. 

The (a, y)- and (x’, y’)-axes are chosen to have a common origin, and this 

is taken to be inside the boundary curve C of the aerofoil. Thus, defining 

r+ rv’ +iy’, the relation z’ = ze®’ holds at all times, 6 being the 


angle between the (2, y)- and (z’, y’)-axes at any time. 

The problem to be solved is essentially one in which two prescribed 
boundary conditions have to be satisfied, one on the fixed wall and the 
other on the boundary of the moving aerofoil. Also the stream function % 
must reduce to zero at infinity, where the fluid is at rest. Clearly we require 
that / be constant on the plane boundary. Further, if the aerofoil possesses 
velocity components w’, v’ along the 2’- and y’-axes respectively, and has 
an angular velocity w about the origin of coordinates, then it is well known 
that the boundary condition on the aerofoil to be satisfied by ¢ is 


ee 9 


ub u'y’+0'x' + 4w(a*+y’). 

The boundary conditions are thus known, and must be satisfied simul- 
taneously by w. Defining the complex potential Q by Q = ¢+-1, it is 
therefore clear that Q must also satisfy two boundary conditions, and in 
addition must reduce to zero at infinity. Bearing all these conditions in 
mind, we follow Green (1) in choosing Q so that conditions both on the 
rigid wall and at infinity are automatically satisfied. The condition holding 
on the aerofoil boundary is then used to determine the various arbitrary 


| 


coefficients occurring in the expression assumed for Q. 


3. The formulation of Q 
To define the aerofoil considered here, we adopt the transformation 
ly 

“a 
it > ant (3.1) 
which maps the boundary of the unit circle |¢ | in the t-plane on to the 
boundary of the aerofoil in the z-plane. The region outside the aerofoil C 
transtorms conformally to the region inside the unit circle y in the ¢t-plane, 


while description of C in a positive sense corresponds to description of y 


In a negative sense. 


pd 
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Assuming circulation K = 27A, around the aerofoil, a suitable complex 


potential function is 

| on A a 
Q = 1A, {log z’ — log(z’ + 21b)} + > [4 nt 3.5 
z 02 g( vb)} l \(2’)” (z | 2ib)”| (0 


to 


n=1 
which automatically makes the boundary wall y’ b a streamline, and 
vanishes at infinity. (Throughout this paper a bar placed over a quantity 
will be used to denote the complex conjugate of that quantity.) The method 
of solution is now similar to that given by Green, and we rewrite (3.2) in 
the form 


Q = iA,logz’ 4 > St + Byte’) (3.3) 
n=1 ig 


omitting terms independent of z’, which do not contribute to the forces 
acting. Here 


B,, = "a Agb > "a, A, bm") (3.4) 
r=1 
gntl In+-r—1\ a" 7 
and °_ = ———, = —— (n,r>1). (3.5) 
n.2” n  )Qnsr 


To satisfy the boundary condition on the aerofoil we need to express 0 
as a function of t. Using (3.1) and omitting terms independent of t we have 
the results 


(z')" = > "B, i, = logz’ = —logt+ 5 Bt, = (2’)kr > vy, t”, (3.6) 
0 


n=r n na=-—? 


where the coefficients °8,,, "B,,, "y,, are easily expressed in terms of a,,. Thus 
in terms of t, Q becomes 


wo 
Q —iA,logt + 2, (Cnt + D,t-"), (3.7) 
r= 
where 
n o x 
¥ — 4,0 - % " > 92 Q\ 
Ci, = 0 B,, Ay r Zz "By, A, " p "Yr B,, D, > "y n B,. (v.5) 
r=1 r=1 r=n 
The boundary condition on the aerofoil is equivalent to 
im{Q+ #2’ (t)—hiwz'(t)z’(#)} = constant (3.9) 
on the aerofoil t# = 1, where w’ = u’+iv’ is the complex velocity of the 


origin of coordinates, referred to the z’-axes. 
Now on the aerofoil 
, - ~ x > 
2(t)2'@) = Sb, + 3 5,8 ‘| 
; (3.10) 
where RK =— > de 
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nple a Hence, equating coefficients of t” to the conjugates of coefficients of 7”, 
39) gives the relations 

2 C,+@'ay4,e% = D,+twb, (n> 1)) 

C,+7'a,e” D, +-iwb, + w"e iG, \ 


We now need to find the coefficients A,,, and following Green we assume 











nti i 
tho that 
aes A > "A,b-* (n>0). (3.12) 
me) I k=0 
Thus from (3.4) and (3.8) we obtain the expansions 
B p 19 °A; > mr." A, b rib n+k) (3.13) 
OTCeS ie v ? 1 
L k 
02 0 N \ } N . a, 0 k_ 
4 ( > iB, °Ay PM 1;.;5 ie a 2 Aj, 
’ k=0 r=1 k=1 m=1 
k-—1 r 
Ss ' > ny 0 k-r k 
> "Y n “k-r 1, m0 
) } k=2 r=1 m=1 (o ] 4) 
D _ S No” { b-? 
SS 42 k-1 r 
. 2 : ma n k fr k 
la > = > y _" Xj r A, m b 
n+1r=n™m™=n 
The boundary conditions (3.11) must be true for every value of b, so that 
(3.6 using (3.14) and equating coefficients of corresponding powers of b~', we 
obtain the relations 
Thus a. : — 
6B, 9A, +18, Ag+" va, iwb, +w'e-"G, 
; : , (3.15) 
Y 06 °A, > "Ay tweed, 41 iwb, (n> 2) | 
r=1 
? i s 
04, Sv 7B { y 1y,9A, }-Y-1 % Ay (n 1)) (3.16) 
8 r=1 | O (nm 1) 
k-1 fr = 
\ » 7 0 P ~ May Mm k-r 
{ P a A, 3 XQ A, m > > "Yn Xi—r A, m 
r=1m=1 
q — 
| > Xg- f > 3 yr . Xfe—r ‘ 1, m (A 2 h n) ° 
1 ne t i n 
| in An (2 2) 
0) k 2) 
| © 38 
| (3.17) 
: 10 By means of these equations the coefficients "A, forn > 1 can be expressed 


| successively in terms of the circulation coefficients °A;, and then the 


potential function Q is completely defined in terms of these coefficients. 
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4. The forces and couple on the aerofoil in terms of Q 
The force components Y’, X’ and the moment of the couple [ about th 
origin are given by the integrals 
Y’+7X’ | p dz’, r = fe | pz’ dz’ (4,] 
Gg 04 
taken round the boundary curve C of the aerofoil, » being the pressure at 
any point of the surface. Here the dashes denote that the particular 
quantities in question are referred to the (a’, y’)-axes. The pressure distri- 
bution throughout the fluid of density p is given by 
eQ @O 
Pp constant — $p(u?-+-v?)- bol 4 =) : 
, C 7 CT 
where wu, v are the components of the velocity of the fluid referred to a set 
of axes fixed in space. But the axes (a, y) to which Q may be referred are 
moving axes, and if we make our fixed axes coincide instantaneously with 
the position of these moving axes, then 


eQ ; . . re CQ eQ 
= —u+iv, ie. w+v? = - ie 
Cz Gz OF 
; a 
and since z zel? Ze-'?, we also have 
. a. CQ eQ 
ur v" 7 a=° 
C2 Cz 


Again, by employing well-known methods (cf. (6)), we find that 


é a2. a a .an.,.i 
— = — (we-”-+ wz) - 1 w 1a — 1. uy La - 
c 


a 


. Q, .@Q 


OT oz co ob ow’ Ww Cw 
where the first two terms take account of the moving axes, and the remain- 
ing terms allow for the fact that Q is a function of 0, b, w’, @’, and w, which 
may all vary with time. Allowance has been made for a variable circulation 
K by assuming it can be expressed as 
K K (w’, 7’, 0, w, b). 
Thus the pressure equation may be written in terms of ¢ as 

a constant —‘ = “\S a) (w’e-*94 ioe) (5 | 
p ct ct\dt dt ot \dt 


=) 40° .. €Q)/dz\-} c re yee ere we = ited 
—(we wz2z)—|—. + | w + vy 1 ;t+w —+wo— || 
ct \ dt o6 cb Cw cw CW 
and then the forces are given by 


Y’+7i1X’ = [ pdz’ = — | p 


@ 
( y 


. Gz 
and the couple by I —re | pz dt. 


} 





MOVE 
By 


© 


and ¢ 
These 


where 


wher 


Ne 


abou 


whe! 


and 


whe 


Alse 


whe 


t 
0ce 
qui 
inte 
on: 
wit 
has 
poi 
at 1 








MOVING AEROFOIL IN NEIGHBOURHOOD OF PLANE BOUNDARY 405 
By using expression (3.7) for Q it is now possible to calculate the forces 
ut the and couple. In the first place we need several expansions about t = 0. 


These are 





} / l : A 
| f? N c. te 
at n=0 ™ 
— } 
ir where 
p1cu 9 2 2 9 
“"e « <a, Ig Ay Ag 5 (4.2) 
dic Co Cy () Cy 29 C3 =". C4 - - 3 
ai az ap 
, 2 21 ¢ 3 
f lp +-A5Q3) 15 Ag+ 4) a5 64) Ay A444 a aad 
ao ag 
l ~~ 
S t p2 g,, t” | 
PC | C n=0 (4.3) 
7 W “ 
where Zz aC... 
? 0 
Now on y, —=€%Y l)a@,t-", and so we obtain the expansions 
Ys a 
bou f () 
dz\ . . , = , 7 
ms) dt > d,t " 
dt / n=0 n=1 
(4.4) 
where d S (7 | a... d,, > (n rr ] a, ot | 
7 0 r=0 
und =| A) Sh w+ > ht") 
mM dt} Ll n=0 at | ais 
y] Ps (4.5) 
where ] > (r—1)¢,,.,4-, h', > (n+r grin se | 
=0 r=0 
dz\-1 dz et < y 
Also 2 $4 S tt T t_t-*) 
\dt al n—=0 n= 1 . 
(4.6) 
7 
where ] S } ] L, » h, a > a, vx | 
on r 0 r 1 
Using these results we may proceed with the evaluation of the integrals 
occurring in the expressions for the forces and couple. The integration is 
juite straightforward, and the integrals have poles at t = 0. When the 
integrand has a factor (dz/dt)-1, we note that (dz/dt)-! is regular within and 
ny provided that all the singular points of the transformation z z(t) lie 
thin the original aerofoil contour C. In the usual case where the aerofoil 
has a shat Dp ti viling edge, the transformation ceases to be conformal at this 


point, so that any integral involving (dz/dt)-! must include a partial residue 


t the trailing edge. This is avoided in the ordinary way by choosing the 











circulation so that éQ/ét 


usual Joukowski condition. 


By employing these methods we find expressions for the forces 


couple in the form: 
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Y’+4X’ ; a - 
——-- = —8A, > ne, (C,+D,) 
pate 0 n=1 : 
oO =] ; : - aoe | 
‘ > Cy 1) > r(n T r)(¢ »< n+r i D, Dy +r)| 
n=1 r=1 
. ("“< 7 eT 
— >, 1) 2 7(n—r)C, D, r| — te Ag+ 
n=2 r=1 
we —iAgd, + ¥ nd,,5C, 
n=1 
a > nd, 2 D, D,d;)+ 
n=2 
+iw(tA,(a,—h}) > C,hnii—D, hy 
n=1 
— +1 D,,) T C1 ay) T 
C C C oe C 
+- |W 2 u , u + W ~ 4 
8 ob ou ou Cw 
losa (= = = 
)22Ag(ay + > a4,)—a,(C,+ D,) 
nn = - 7 
Dd na, 44(C,, -D,,)| 
n=1 
and 
LE e y “ ry | 
= — im| S In| 2, r(n—r)C, D,_, 
Jtn—tl . r ~ Y Ue | 
ee eh + 2 Cn In +-1—Dy hn 1)} — 
nm 


—)> Inj iA,n(C,+D 


vt n 
n=1 


=8 
Loar oiOls, r 
wD’ ¢ (Aga, Da, | > 


—iw te 1 (by - - 1) 
n=1 


lo; 
x {2A 
| ai 


Vv b 
7a n 
1 


U 


)+ Dd r(n I r)(C,C, 4,+D, Dre) + 


y n(C,l,—D,1,+C, b,—D. b.))— 








- 0 at the trailing edge. This is, of course, the 


(4.7) 


(4.8) 


and 
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MOVING 


By using equations (3.11) and (3.14) in conjunction with (4.8) and (4.7) it 
is possible to express the forces and couple as power series in b-, i.e. we can 
write 
(4.9) 


Y’4+iX’ = ¥ (¥;,+iX,04, 


> [.. b-*, 
U k=0 


This process will not be carried out for the general aerofoil defined by 
3.1) as the expressions become so long and complicated that little of 
relevance can be deduced from them. Accordingly we shall consider later 
the expression of the forces and couple as power series, when three specific 
types of aerofoil will be considered in order to assess the effects of thickness 


and cambe1 


5. The Joukowski condition 

To discuss fully the forces acting on the aerofoil, it is desirable to obtain 
a value for the circulation. The usual procedure is adopted of choosing 
the circulation so that the fluid leaves the trailing edge smoothly, and 
since the transformation (3.1) can always be chosen so that the trailing 
Oatt —1l, 


edge corresponds to ¢ |, we have the condition ¢Q/ét 


From (3.7) this gives 


A-  ° 


A, ww e~-“"a 2 


’ : . 7/18 
])’ n(D,, dD, : iwh,, wea, +1) 0, 
1 


using the boundary condition (3.11). Thus we have: 


1)"nb,,, 


) 0 0 y= 
n=1 
1, Udo € ‘ Ao Ly 
i 
( \ l = = 2 ] 0 
r fot z \ l) am XG A, m ses . Xo A y-m) 
1 fh 1 r - 
{ 1 = n= k-r _m, m r 
> > z \ n(n) n Mer 1, m Y-n “k-r 7 
1 
(k > 2) 


equating to zero the coefficient of each power of b-!. 
The complex potential function Q, and hence the forces and couple, are 


now uniquely determined 


6. The flat plate aerofoil 
Here the transformation (3.1) assumes a particularly simple form, and 


if the plat 


- is instantaneously inclined to the plane boundary at an angle @, 
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then our transformations are 


z= —(t-!+-42), z’ = zet?, 


(6.1) 


l being the length of the plate. Expanding the forces and couple in the | 


form (4.9) we find that: 


i [2 


* = lu'(u’ sin @—v’ cos 6) — — cos 0{2(w’ sin 9—#’ cos 8) +l? 
8 L 


y; . ; ‘ oo ; l = 
ee (w’ sin @—v’ cos ay! (2u sin 6—v’ cos #@) + w-—(1- sin2@) + 
7p 4 4 
sy oo x bi. 
+ oj sin @—#’ cos 6)+-a lain 20 
1 | 2 
Se aere 9 a ae 
== 16 sin @—v’ cos @){u'(1+-3 sin?@)— 2v’ sin 26} + 
7p ) 
»(t\5 ee 
+ 4w* i) (4—3sin*6@)sin 6 — 
wlt |, si MOV Jail ciao -n20)) 
oa ju’ (1—9sin?6+- 4 sin4@)+ 3v’ sin 26(9— 4 sin?@)}— 
) 
iy, Pe secs ss sis 
—4w0|{—) cos?@—2{-) (wv sin @—Z@’ cos 6)(3—2 sin20)cos 6 
4 4 
xe li : P. er 
—? — —ly'(u' sin@—v’ cos @) + —sin 6{2(w’ sin 8@—Z@’ cos 0)+-la} 
7p 5 
X; Pos ; ; l : 
a. (w’ sin @—v cos 6)| v’ —w— cos 6}sin @— 
TT p 4 4 
]\3 
~(3) sin?{4(2’ sin @—0’ cos @)+-la} 
Xo ‘a ; y\® . ‘ 
= ——v'(w’ sin @—v’ cos Acos*?—bu*{ sin 6 sin 26-4 
Tp 16 4 


P\*.... 5 ’ 
| 4(3) (usin @ cos 6(1— 4 sin?@) + v’(5— 4 sin 20)sin26) 4 


)"— ae - (1\5 
+2(3) (w’ sin 6—7#’ cos 6)(3 —2 sin?@)sin 6-4 2i(() cos @ sin 26 








» (6.2) ) 
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and 
’ 12 ] 
6. [) (u' cos @ v’ ain OY (eu sin # v’ cos @)—w \ ~ 
a 
In the - 7 o7 > l . 
ol _ (w’ sin 6—v' cos ¢)— 6a 
16 4 
[3 . 
t’ sin @—v’ cos 8)(w’ sin 20— v’ cos 20) 
0 16 
1\ ae 1 or - ; 
fev? cos fi sin @+ 4{ (2 sin @— 7’ cos @)sin 6 
) } 4 4 
2“) (w’ sin @—v’ cos 0 - (6.4) 
| ap \4) 
fu’ cos (1+ 10 sin?@)—v’ sin 6(9— 10 sin?6)} 
w| fu’ cos 6(5—34sin7@)+-v’ sin 6(37 — 34 sin?@)!- 
4 
6 > 
6a? | sin 20 | } (Ww sin 8 6’ cos 0)(5—2 sin?6)- 
I\e 
ay (5 ?sin°6@) 
} 
7. The symmetrical aerofoil of small thickness 
Consideration of this aerofoil enables us to assess the effect of thickness 
onthe results. The appropriate coefficients a, have been given by Fujikawa 
5) in the form 
] 
l l l 
(1+. ay € As (1—e) 
4 4 ss 4 
(7.1) 
[ Ez: 
) l e(1—2 A, e a O (n 4) 
{ 4 
where ¢ is small and real, and / is the aerofoil chord. It may easily be shown 
that the thickness ratio of the aerofoil, defined as usual by the ratio of 
: F : ’ 3V3 ce 
maximum thickness to chord length, is given by ; e(1—e). For most 
purposes it is sufficiently accurate to neglect the square of «, but we shall 


retain it for purposes of comparison with Fujikawa. 
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The power series expansions for the forces thus become: 


Y; Non! ei 
® — [(1+«—e?)u’(w’ sin 8—v’ cos 0) — 
7™p 
wig in2 ‘si 
—e“w {w'(15—8 sin?@)+- 40’ sin 26}— 
4 
1 3 
20 €(2— 3e)sin 8 
4 
5) - -/; . 7 24)’ j 
: i} [4(1-++e)(%’ sin @—#’ cos 6)cos 0+ 2¢? w’ sin 204 
+i’(3-+-4 sin?@) + 2a] cos (1+ de—§ e)] 
ry ]2 


—— (1+ 2e—e?)(w’ sin @—v’ cos 6)(2u’ sin @—v’ cos 6) 


1\4 ' 
- Seu (2—sin?0) 
4 


]3 
—W 1(1+2e)(1 


sat 1 sin?4)(w’ sin 8 
) 


v’ cos @) 
—e* wu’ sin 0(11—3 sin?@)+ e?v’ cos 6(7 


1\8 ; ; 
+-2(1+2e+ (3) ld’ cin’ ®’ cos 0) ‘alk ts 20 


| 4| 


3sin?@)'+ 


YF ..; . ae ee ros 
2 — 16 sin @—v’ cos @)fu’(1-+-3 sin?@)— 2v’ sin 26-+- 
7p ) 


+ dew'(7+ 18 sin?@)—5e v’ sin 20+ Je?u’(3+-4 sin?6)} — 


it... ‘ ; , , 
2u( |) [2a'(1—9 sin?6-+- 4 sin*@)-+ v’ sin 26(9—4 sin?@)+ 


-ef{w'(1—32 sin?@+-4 sin*@)+-v’ sin 260(17—2 sin?@)!+- 


+e*{u'(59—8 sin?@)sin?6—v’ sin 26(53 — 40 sin2@)']+- 


a", ; 
- w*(3) sin 6[4(4—3 sin?@) + e(21—8 sin?@)— 
—e?(99—80 sin?) ]— 
1\5 
a3) cos 6{4 cos?6 + 8e—e?(15—32 sin26)! — 
a 
-(;) [(2%’ sin 6—#’ cos @)cos 6{2(3—2 sin?@)+ 1le}4 


+€*{ — 32’ sin 20(1—8 sin?@) + #’(5—27 sin?0+ 24 sin’6)}] 
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é e”)v’(w’ sin @—v’ cos #)- 
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De| w?(2— 3e)cos 6-4 
\4) 


;) | 4(1+-€)(u%’ sin @—v’ cos 8)sin 0— e?(7 — 4 sin?0) a’ — 


2¢7’ sin 26+ 2H] sin 0(1+ de—e?) | 
l\* 
2e—e*)(w’ sin —v’ cos 6)v’ sin 8 42" sin 26 
4, 
[8 a ' 
(w sin 6—v’ cos 6)sin 26— 
2¢ w’ cos 8 cos 26+ v’ sin 4(3 2 sin’)! 
«*{3u’ cos 6 cos 26+-v’ sin 6(1— 6 sin?6)}]- 
i } are oe 2A > oe al 2 9 coin2 \ 
—(u sinéd—v cos f - 6—ecos 26 + . (3—2s1n*@)} — 
16 a 
1\4( . «2 
tei| ! sin20 € cos 26 (3-4 2 sin26)| 
$} | 2 
bis ) i> | 29 on 9 oin29 
uw siné v cos 6) cos*é + —(7—2 sin*@)+ 
») 


Ee“ 


+ —(3+4sin?6 + 
4 ) | 
l 4 ° ° pie ° 

2en| ;] | w’ sin 26(1— 4 sin?6)-+ 2v’ sin?0(5—4 sin?@) + 


wu’ . . 9 , 9 . 
é! > sin 26(5-+-4sin*@é)+-v’ cos?6(3-4 4sin29)\ + 


| 2 
e*} — 5w’ sin 40+ 2v’(6— 27 sin?0-+-20 sin) | — 


({t\5 
| J cos 6{12 sin2@+-«(5-+-8 sin20) +e 


9 


(17—80sin?@)'!+ 


|=) sin 6{4 cos*# 2e—e?(31 —32 sin26)}4 


]\ 4 
(3) (% sin @—’ cos @)sin 6{2(3—2 sin?@)+-5e} 


eu (24-19 sin*6— 24 sin4@)+ $0’e? sin 20(19— 24 sin?6)] 
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8. The circular arc aerofoil of small camber 


Green (1) has obtained the coefficients a, for a circular are aerofoil of 
chord / inclined to the boundary wall at an angle 6, and subtending ap | 





angle 4x at its centre. The camber o of this aerofoil is, of course, given | 


by o = }tana. Applying our restriction that the camber is small, we may 
write «, = tana and neglect third and higher powers of ¢,. Thus the 


coefficients a,, become: 


dg = — —e,(1+2e,), a, el, 


_— 
_ 
a) 


After suitable reduction the expansions for the forces Y’ and X’ become: 


re ; ’ ] 3 
Yo = lu'{(u’ sin @—v’ cos #)+,(u’ cos 0+’ sin 0)! + 2a ] €, cos # 
7p 4 
/? ‘ ° 5.80) 1 of 
+ w il sin 26—v’(3—2 sin?) + 2e,{w’(3—2sin20)+ v’ sin 26] 
8 
(w’ sin 6—#’ cos @)cos 6-4 as cos 6-4 
i " a , 
te, 46 sin 20+-2’(2—sin?0)+-a sin26| — 
¥ 4 
-7{3w' sin 20—3¢'(3—2 sin?0)+- wl cos 0) } 
¢ . jos : . 
— = —-—|(u' sin @—v’ cos 6)(2w’ sin 6—v’ cos 6) +- 
Tp 4 


+e,(2u’ sin 20—3w’v’ cos 20—v’? sin 20)+ 
P 
-}(4u'? 2 cos?0— 3w’v’ sin 26-4 Dv’? sin?) — 
3 
w 76 L@ sin 6@— v’ cos 6)(1+-sin?6-+ 2e, sin 26)- 
) 


+ eF{w’ sin 0(15—10 sin?6)—v’ cos 0(9— 10 sin24)"] 


> l . * 9 
a €,{sin 20+ 2e€,(3—2 sin?6)} 


T\3 1\ . 
ee l w' sin @—#’ cos 0-+- w® }sin 20 
4 4 


1 de, (3—2sin*@)sin 6— 20’ cos8@+ w — (2 sin20)| 

| 4 
ee ee oe ee 
— ‘cos 6 cos 20-+-3’(5— 4 sin26)sin 0 al asad’ 
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Vy w’ sin @—v’ cos 6){u’(1+-3 sin?@)— 2v’ sin 26' + 
ofoil o¢| 7 16 
ling ay L os 6(1+- 10 sin?@)—w'v’ sin 0(17 — 26 sin?@) 4 
» Liver + 2v"* cos 6(3—8 sin?6)}+- 
ve may ' 2 . 9 iy " a 
YY | Le#{w'? sin 6(13—8 sin?@)—w'v’ cos (17 — 40 sin?@) — 
us the 
4v’? sin #(7 —8 sin*@)}]— 
2w{—) [2u’(1—9sin?é+-4 sin4@)+ v’(9—4 sin?@)sin 20+ 
\4 
8] €,,—w’ sin 26(11—6 sin?6@) + $v’(31— 60 sin?0+- 24 sin‘@)} 
) 
é 20 — 21 sin*@+- 12 sin*@)-+- $v’ sin 20(23 — 12 sin2@)+] + 
Come ti=| -| (4—3 sin*@)sin 6+-4e,(13—4sin?6)cos@ 2e} sin é]— 
1} 
| } | 2(u' sin @—Z’ cos @)(3—2 sin?6)cos 6+1a cos? 
) } 
€; 4— 9 sin*6-+-8 sin4@)+-40’ sin 26(5—8 sin?6) — 
— tal (3—4sin?6)sin 6+ 
€* ’ sin 26(23—8 sin?6)-+-¢’(20— 27 sin2@- 8 sin4@) — 
3wlcos 6} | 
and 
X = \3 
= t sin@G—v cos@) €,(U cosé-+-v sing); 2<,0°(;) sin@ 
2 sin*@)—v’ sin 26+ 2e,{w’ sin 20-+-v’(1+-2 sin?@)}] 
| t¥5 
}=(¢% sin@6—v cos 6@+-w--]}sin 6- 
Q | >) 
| 4 : iu “),./ . 2A ® ‘ | 
“\* sin 2 2v'(1+-sin*@)—w ; Cos 0) — 
u'(1+-2 sin?@)— 30’ sin 20-4 kissin |, 
™ 4 sin @—v COs @)sin 6+-€,(u sin 26—v' cos 26) 
e} cos 6(u’ cos 6 L-w’ sin é)|— 
8.2 —J|(w' sin #—v’ cos 8)(cos 0@— te, sin @)sin @ 
5u’(1+-4sin2@)cos 6+-v’ sin (3 20 sin%9)}| + 
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l\4 a : 
- te, 0°(3) {eS 2 sin?4)-4 2e, sin 26} — 





at Pr ” _ 
—— | |{w' sin 6—*#' cos 6+ \sin2@- 
16 4) 
+e, sin 6| a’ sin 206—2’ cos 20+ a — cos 8) + 
4 
-|. F| —w (1+ 8sin*8)sin 6 --#’ cos 6(3-+- 8 sin?@) 
l tan 
—w—(3- 4sin26)| 
4 
X% ae ver ; , 
—2 — —__y'| (u’ sin 0—v’ cos 0)cos?6 + (w’(1—6 sin?@)cos 64 
7p 16 2 
+-»’(5—6 sin?4)sin 6} — 
oe {3u'(1—8sin2@)sin 6+-v’(5+-24 sin2@)cos 6! | +- 
4 j 


i\47 , ; ; ; : 
2u( 7) E (1—4sin?@)sin 20+-2v'(5— 4 sin?6)sin26- 


€ , a oe ° te ~ * 9 
*{—u’'(10+-38 sin?@— 48 sin4@) + 2v’ sin 20(5— 12 sin?6)} 
41 


€j{2u’(2—3 sin?6)sin 20+ v’(5+- 22 sin?6- 12sin*9)}| 7 


eft? ; , P ° 

w* J [6 sin @ sin 20-+-e,(1—4 sin?@)sin 6— 2? cos @]+- 
4 

9» 


+2(3} [(%’ sin @—#’ cos 6)(3—2 sin?@)sin 6+ 11 cos @ sin 26-4 


+ fe,{—w' sin 20(1+ 8 sin?@) + 2%’(5-+-8 sin20)cos20 


1@1(3+-4sin?0)cos 64 
19 sin?@—8 sin‘) + 


+ 20 (15—8sin*6)sin 20—l sin 6} | 


9. Discussion of the results for the three aerofoils 
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We may now discuss some of the more interesting cases of motion, making 


the stipulation that the velocity component of the origin in a direction 
perpendicular to the aerofoil chord shall be small compared with the 
component along the aerofoil chord. This is in order that the aerofoil will 


not become a bluff obstacle, as it is known that the Joukowski condition 
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sives results in good agreement with experimental results only when such 


a stipulation is made. We therefore require that (w’ sin @ 


compared with (u’ cos 6+-v’ sin 6). 


v’ cos 8) be small 


Again, the expansions for forces and 


couple are actually written in powers of //b, and hence are more likely to 


be convergent for 1/b 


|, so that we shall consider only values of b greater 


than the chord length. Further, any conclusions we may draw from the 


results must be treated with some caution as only three terms of the series 


ire given 


searing all these conditions in mind, we consider the following cases: 


(ase g iu constant. vy’ OU, Ww 0 


This is equivalent to the simple stationary aerofoil problem which has 


been discussed by various authors. It is easily seen that the expressions 


for this case obtained here agree in all respects with those obtained by 


Tomotika, Green, Havelock, and Fujikawa, so that there is no point in 


reproducing their discussions of this case. 


It is believed that the following cases include some new results. 


Case 2. v’ constant. u’ UO. w i) 


If the aerofoil moves directly towards the wall, we write @ 


x), 


: being small to conform with our stipulation regarding the velocity com- 


ponent perpendicular to the aerofoil chord. For the flat plate we then have: 


il) *sin*a{ 008 )  O(b-) 
X l/l l /] . | 
\ iN) |cos , ela) sin2a+ O(b-3), 
| 1 /1\cos 2a mit 

| : 1 COS « saz) LOsin“«) | 


Xo 


O(b-3), 


aplv’* sin « 


9 


i mp—v *sin2a 
0 = . - 
PS 





(9.1) 


It follows that one of the effects of the wall is to induce a small force Y’ 


ting on the plate and directed towards the wall. The sign of the force 


is inde pendent oi the sign « 


approaches nearer the wall 


or negative 
For the symmetrical aerofoil we find 
} ai ‘ 
- AVAL *sin*a) 1 + 2e—.*) 
X vt. / : ; rive 
X’ AVAL jad aici raz) ss 


f x, and Y’ increases in magnitude as the plate 


Further, X’/X}, T/T, 


| for all small positive 


and both expressions increase as b decreases. 


3?) J O(b | 


il | 
(1+ 3€)COS | _O(b >) | 
(9 ,) 


so that the effect of thickness is to increase both Y’ and X’/ Xj in magnitude. 
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The expressions for the circular are become: 





ll = a 
— = —-|-]v*| (sina+e, COs «)* + —| sina COS ~-+- 
7p 4\b, b | 
, ; 
as + sin «(5—8 sin?x) +3(] —8sin?x)cos | | |. O(b-3) 
» en 
xX’ l Vis 
—_ — 1+-1(cosa—e,Si — sin?~— 
x, | (cos a—e, sin x) 5 vals) | in? 
(2¢,(1—4sin*a)cos a ef sin a(29—24sin2x))] , O(b-3) : 
(sin a~-+-e€, Cos x) 





Thus for small positive values of «, the force Y’ induced by the wall is 
considerably increased in magnitude by the camber. Considering X'/X, 
we see that X4/X4 is affected very little by ¢,, while, since sin « and e, cosa 








are of the same order of magnitude for small «, X3/X% is considerably} 


diminished by the camber terms, and may even be changed in sign. | 
Hence camber tends to diminish X’/X{, below the values it assumes for the 
flat plate. 

These conclusions are best illustrated by considering a numerical case. 


Thus for « = 5°, for the flat plate we have 


\ 


a ply’? 


= 1+0-2491(;) -0-0005(; F icine 
Xo b b 
1+0-2472() 0.02891, Pb eas 
. b v 


while for the symmetrical 


l \ 
—_{0-001901+-0-001894 —]} -+L..., 
b b} ° 


aerofoil with e = 0-1, i.e. thickness ratio 0:117 


} = 2. 0-002266 0.002368" S ee 
mplv : b b 

~ 1+0-2715(, -o-0156(;) 

a - hb] - \b, 


and for the circular are with «, = 0-1, i.e. camber equal to 0-05, 


r’ l —— 
_ — —_{0-008725- 
aplv’ b 


md = 14 0-2469(7) -0-0558 : en 
Xo b, b 


ee 
0-005166 beiccas 
b 
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Case 3. uw’ constant, v’ constant, w 0 


For the flat plate the results are: 








| ) | ‘ ih ] 
ie | {2 sin 0 cos 6) a 
} . } u b 
a | = D7’ . 1\2 
-/(1+-3sin*6é) — -sin 261 ( 
16| u J\b 
| Y=n, sin @—v’' cos 
; 
. 
Y ria l 2p 
sin 0 cos*t 
i \7 TAY 
seantis (9.4) 
X'/X i. solu’ (au’ sin @—v’ cos @) 
no . l _ sin 26 v cos 26 
' = 
| }\/ ucost » sing 
- l/l "| u’ cos 6(1+ 10 sin?@)—v’ sin 6(9— 10 sin?@) 
(') , é. F 
32\b) | wu’ cos 0+-v’ sin @ 
[? : — 
TO sin 0 v cos l)(u cosé+yv sin @) 
{ 
there is no velocity component perpendicular to the plate, i.e. if 
sin 4 cos @) = 0, then both the forces and couple vanish (at least to 
order of approximation 
Consider the case where @ is small, (w’ sin@—v’ cos @) small compared 
cos @+-v' sin @), and @, v’, uw’, (w’sin@—v’ cos @) all positive. Then 
is b decreases from infinity, X’/ X94 first decreases and then tends to increase, 
0-117 | possibly to a value above unity. With Y’/Y, we may compare the 
esponding case with 2 0, and we see that there will be less initial 
rease of YY, and a greater final increase, as 6 decreases from infinity. 
Since (w’ sin 20—v’ cos 2 2(u’ sin @—v’ cos #)cos 6+v’ > 0, T/T) also 
ibits an initial decrease and final increase as the plate approaches 
er the W 
Since 0 (by assumption) the aerofoil is actually moving away from 
wall, and so it is perhaps better to state our conclusions in terms of 
If then we assume that we commence with some minimum 
le of we may say that Y’/¥4, X’/X9, and [’/T, all decrease from 
es which may be greater than unity to values less than unity, and 


me the value unity as the distance of the plate 
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The forces on the symmetrical aerofoil are: 





- ] : y’ l 1 /l\*{ . 
—— = | —_(]-+-e—e*){ 2sin 0 — cos 6) BIE : 1+-3sin?6) + 
eee ag 
2.2 fe e? a v’ | : e*\ . | 
| 5+ 12sin2@) —— (3+ 8 sin20) — 2 L «i ge —— nin 99 
3° aa i ) 
" oa . (9.5) 
xX Bs 1 /l\? . 
1 —i(1+-e—e?)_sin 0+ | 10820 Se — 3¢7)4 
x 4( e—e ); rely (cos*6+- 3e— ye") 
Y), = zplu’(u’ sin @—v’ cos 0)(1-++«—e2) 
Xo mpl(1+-e—e*)v'(w’ sin 6@—v’ cos @) 


Under the same assumptions as for the plate case, the aerofoil thickness 
will emphasize the tendency for X’/X{, to decrease and then increase as 
l/b increases. Writing Y'/Y, in the form 
Yy’ 1 /l as I - € 

Z . +} 8(1-+-e—e*)sin 64 I --3sin?@) + — (5+ 12 sin*@) 
Ys 16 b b\ Y 


9 


€ 


4 


(3-4 8 sin20)! F - Ma | ¢—e?) il} 1 Se s)sin ol cos 0+ O(b-) 
} 4u'| b\ ‘. 2 1b 

we see that, for all //b, the term involving w’ is numerically increased by the 

small thickness parameter ¢, and this will tend to compensate for the 

decrease produced by « in the terms independent of v’, for small //b. Hence 

it seems likely that Y’/¥4 will behave much as it does for the flat plate. 
Neglecting ¢{ for a first approximation, the expressions for the circular 

arc become: 


y’ a. | 7 
— = |] — lo sin 6 “eos 0-4 €,{2 cos 6-4 ’ sin | 
Y, 4\b} | u u | 


L(Y") 1-4-3 sin?) — =” sin 20 
ahs +5 sin* sin 2¢ 
' 16 , u’ 


€,{u’? cos 6(1—4 sin?@) + 


+w'v'(1] 


12 sin?@)sin 6— 2v’2(3— 4 sin?@)cos 6! 
2u'{(w' sin @—v’ cos @)+-e,(u’ cos @+-v’ sin 6)! 
+ 0O(b-*) . (9.6) 





Yo = mplu’{(w’ sin @—v' cos 0) + ,(w’ cos 0+’ sin 6)} 


.; l 
— = 1—}(sin@+., cos 6)- + 
t vl t( 1 ; i 
| (L)*P gog89— 2ex{t ¢086(1+-4 sin*®) —v’ sin (3—4sin*#) 
16 (w’ sin @—v’ cos @)+-e,(u’ cos 0+’ sin @) 


b 
+ 0O(b-%) 


Xo = —7plv'{(u' sin @—v’ cos 0) +-€,(w’ cos 0+ v’ sin @)} 
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iptions we find that the camber tends to keep 


0 and _ : bel 


yw | for a considerably greater range of values 
aerofoil. 
10°, the expressions for the flat plate become 


> 


0-07 175(") 0-06602(") 
0-0 1303(;) 0-06063(( ) 
l 

5) 


d/} 





l 
0-07450(- | 0-O3800 
) 


0-1 we then have 








]\2 
0-O8148 0-08194| 
} b, 
— oe ft 
004745 0-07576| 
4h h 
lar are e, Q-1 are 


l 
0-1242(7) 0-05: 
} 


0-068] i(") 0-00045(-) Sern 


———. 
= | 
——— 


istant and non-zero 
parison with the case of an aircraft taking 
We re- 


small in comparison with w’, all these quantities 


after a certain minimum value of b. 


sin @—v'cos@) must be small compared with 
The be- 


of primary importance here, and for the 


ive (for a positive vertical force Y’). 


ce 1s 


1 /l ; l 
l | ) Ssing (1 3 sin26)| a. 
Nb l 
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vy’ Po l 
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Under the above assumptions the terms involving w will be negative jn = Thu 
sign, and will decrease in magnitude as b increases. It follows that the 


~ . . . . . Xx’ 
effect of introducing the rotation w is to lessen the magnitude of Y’. zpl 
By comparison with Case 3 it is possible to determine the variation of Y' 
with b, and we see that as the distance of the plate from the ground increases, 

Y’ will diminish to values less than Y{ and then exhibit the usual increase 
to attain the value Y. Thus for w’/v’ = 20, @ = 10°, 
y’ l 1\? 
= = 1—0-0748/_) + 0-0660 = 
Yi, \b b 
genre | 
wl [I : at wl\?/L\? 
anomaes (0-v641 -0-1015—) +.0-0213(- ) =} =p... 
u'\b}\ b u'} \b 
, ; ; . V 
For the symmetrical and circular are aerofoils, the effect of a small 
rotation w may be found by restricting w to be so small that the products a fr 
ew, €,w, etc., may be neglected. The forces will then be modified by w in by 
exactly the same manner as in the flat plate case. — 
of . 
Case 5. u’ constant, v’ variable, w = 0 foll 
For the flat plate the force X’ takes the form 
Xx’ rr ; L/l\ . ify... 
- —ly'(w’ sin @—v’ cos aI — ili) é+ () cos) — wh 
7p 4\D 6\b . 
f Wr 
ar if\. 1 /l\? oe ' 
—- 2 sin 26) ]- sin 6+ —|[-]} (3 2 sin20)| |. O(b-). (9.8 
8 | 4\b " 32\b | | 
- 
An interesting case is one in which the plate describes translational oscilla- ) 1 
rc 
tions in the direction of the y’-axis, and this may be obtained by writing s 
: . 8 
vp’ = vysin nt, vy and nv being constant. Since 
ma 
di det 
= Uv, Sin nT, ent 
ar 
> pe 
vy \ thi 
we have 6 = 611 cost), by constant. 
nb, of 


Clearly by is the mean distance of the plate centre from the wall. Sine 


Vo/nby < 1 we may obtain the expansions 


1 l / 9 : « 
oa Rick. Yo COS NT + O(b5 **, 
b bo nb, 
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Thus X’ may be written as 


x — L/l ivi 
| iad rata.) \- 


l/l 1 /l\* ee 
Ti sin 6/1 }sin 0 ( cos?6|sin nz - 
| L6\b, | 


bie his lL /t\? Pas ve 
nv, sin 26) |—]|sin @ | (3—2 sin?6) + —- cos nr + 
Do 32 by 


to) | ) 4n*be| 

ae ' L/e\. 1 /l\? ; 
—— s1n°@ sin 2n7 r2.cos 6/1 | }sin i cos 20\cos 2n7+ 
S109 7 | ai) 16\D9 


sin 24 cos 3n7t-+ O(b7*). 

We notice the presence in the variable part of X’ of terms which have 
a frequency twice that of the main motion; also there is a term (involving 
b>?) with a frequency three times that of the main motion. Further, it is 
evident that from X’ a thrust can be obtained on the plate. The mean value 
of X’ over an interval of time 27 n (the period of oscillation of the plate) 
follows immediately, and we have 
X l cos Gl] (\sino (~)"  O(b3), 
7 2 | 4\b, 16\b,) J 
where , X’ denotes the 1 n value of X’ over the specified interval of time. 
Writing ,,X, for the mean value of X4, then 


X 1/1 I \? - , ) 
] | \sin @ ale.) Bea wie sl apv5 cos 6. (9.9) 
Now ,,X' is positive the plate experiences not a resistance but a 
propelling force, assuming wu’ positive. This type of motion resembles a 
simple form of flapping flight, and from (9.9) the wall effect on this force X’ 
may be found. For @ > 0 and small, as b, decreases from infinity, ,,.X’ first 
decreases and then tends to rise above the isolated value ,,.Xo. Clearly the 
energy of the propulsion is derived from the work done in producing the 
periodic transverse os tions of the plate. Of course, in actual practice 
this mean force will be considerably modified by the influence of a system 
of free vortices accompanying the plate. 
For the symmetrical aerofoil, the corresponding case yields 
\ = 
\ | r( 1+. af sin @ 
L /t\4{ Oa ee an - oI (},—-8 ( 
cate ! (3-+-2 sin*@) ; (5 saa O(b5 ). (9.10) 
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The effect of thickness is to diminish ,,X’/,,Xo for small 1/bj, but as}, and 
decreases the ratio may eventually be increased by the thickness. r 
Similarly for the circular are we find a 
* . 7p = 
xX’ . l 
m— = 1—}(sin#+e, cos @)—+ 
m0 bo 
1 /1\J 2e, sin 6 cos 20+-€?(1 8 sin*@)co: 36 
- ( a fe amcccichat B : | O(by*). (9.11) 
16\d, | $(sin @ €, Cos @) | 
Clearly the effect of the camber is to decrease the ratio ,,X’/,,Xo found for 
the flat plate. 
sy way of illustration, we find that for @ 5°, the flat plate case gives 
X’ aa 3 
a l 0-02179| | 0-06250( 
m9 bo bo 
while for the symmetrical aerofoil « = 0-1, 
xX’ 1 | . 
= I 0.02375 0-071 i4(; 
m9 bo bo 
and for the circular are e, = 0-1, 
X’ ea . 
“= = 1 0-04670/ | (0-06028 J O1 
m*0 bo bo, fre 
The lift force Y’ and couple may be treated in the same way as X’, ~ 
and for the flat plate we find - 
sn — 1/l\. 1 /1\? — 
u 2 sin OI | }sin 6 {_ (1 3 sin20)| — 
tpl | 2\b, 16\b, 
ad & : 5 a } 
4v9( cos?6{ 1 —— sin) — - 
bo thy - 
: 3/1\. 1 /1\? ; 
u v, cos @sin nz! 1 { }sin 0+ ( (1 7 sin2@)| - 
| } by 16 bo 
I ijt\. 
Vy N COS NT cos?6) | | }sin g ) 
| +\b, 
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t as and 
l L/l\. l (! 2 —e 
sin 26/1 sin @ 4 (1+ 10sin°6)} — 
aol \ ap, | o2 5) ) 
Kode ¥ l/l P | (;) 7 ee 
Up COS U0 sin 0 cos 20 (4 10 sin°@ sin 8 } 
| ‘ls,) 32\b, iid 
ae ; | (. l/l 1 /l\? a 
9 1] Luv, sin nT} sin 6! sin 6 4 }eos 20 (9—10sin*@)sin 6} — 
£\05 32 by J 
. Troe 1 /l\? : 
cos?! ] ( }sin } - | (1 | 10sin26)| + 
2\bo, 32 \0o J 
gives , tear 
NV, COS U COS? 7] i }sin 0+ 
16 | 4 bo) 
Lees: 5 2 ee ° | 
. | } (5—2sin*é) ——_, (u™ sin*6+- v6 cos 26)) + 
64 bo n“D5 : 
; l/l 1 /l\? ’ : 
v2 cos 6 cos 27! sin 6 | }cos 20 | (11—10sin 26)sin 6\ + 
£\ 9 32 bo j 
oa (w’ sin 30 sin 2n7-+- $v, cos 6 cos 20 cos 3n7)+ O(b5%). 
Once again we note the presence of terms of twice and three times the 


frequency of the main motion. Averaging as before, 


3X", Y 1/1 1 /1\?2 
me sin 61 5\p, }sim 9 ala) (1+ 3 sin*)| - 
l 3. , 
| \( _ sin 6)cos*6+ O(b, ), 
bo ld) 
m wal I/l\. 1 /l\? — 
et ee — eee wal lites, as 


L/l via ae ee 
v= COS U/ sind Os PY i Ss 24)s \1. O(b-3 
OF | i A }co 2 ae) ( 10 sin*@) ~_— M(b5 °°). 


For the sake of completeness we give the average values for Y’ for the 


symmetrical and circular are aerofoils. For the symmetrical aerofoil, 


Ie 1 ¢—e*)sin 6-4 
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and for the circular are aerofoil, 


Y’ i se 
m—_ = u’*(sin 6-+-e, cos 8) 


l 
mpl 2 





: 1\. 
_— ( J¢sin 6+-e, cos 6) + 
; bo 





7 2(1—4 sin? 3 6— 3—8sin?6)si 
rm u {(1+-3 sin’) —© 2(1—4sin 8)cos 0 e,(1 : 8 sin?0)sin 6) a 
16\b, | t (sin 6-+-e, cos @) | 


of 1 \| +. aye 
— ail }| (eos é—e, sin 6)?— 
"0 
] 
8 


l ; sas mos j 
( }s cos # sin 24+ 2e, cos 6(1—6 sin?@) — 4c? sin 6 cos 26° | + O(b5). 
Jo} ar 


The general force and moment expressions (6.2)—(6.4), (7.2), (7.3), (8.2), 
and (8.3) enable other cases of oscillatory or accelerated motion to be 


discussed if desired. 


10. Conclusion 

The calculations of this paper have all been based on the simple Jou- 
kowski trailing edge condition, and no system of free vortices has been 
considered. It is realized that for the majority of cases of motion considered 
above, the circulation around the cylinder will continually vary (e.g. with 
the distance b of the origin of coordinates from the wall), and hence 
some sort of assessment of the effects of a vortex wake would be 
desirable. However, the results already obtained with no vortex system 
are complicated, and introduction of the further complication of a vortex 
wake would tend to make the results almost incapable of profitable 
discussion. 

In spite of the approximation just mentioned, it is hoped that the 
results obtained will at least give some idea of the wall effects for various 
types of motion. 

It is felt that to deal adequately with the problem of the effect of a 
vortex wake, a separate paper is required. 

In conclusion, I should like to thank Dr. R. M. Morris for suggesting the 
problem, and for her helpful advice throughout the work. 
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SUMMARY 
\ thin dihedral wing in supersonic flow is given a small velocity in sideslip. By 
rized theory the surface pressure distribution is investigated for wings whose 
form is metrical about the line of undisturbed flow through the vertex. 
s invo dividing the problem into two parts, depending on whether the 
iding edges lie outside or inside the Mach ‘cone from their vertex. For 
case the force derivatives are determined and graphs are exhibited which give 
lue of these derivatives for wings of delta planform. For other wing shapes 
hs are shown which give the contribution from the interference pressure alone. 
ing it is found that to a good approximation the rolling derivative 
ries with t ine of the angle of dihedral and the yawing and sideslip derivatives 
v with t juare of the sine of this angle. As the values of the derivatives at 


known explicitly, this provides a useful method for obtaining the 


1, Introduction 

Ix a previous paper (1) (hereafter referred to as Part I), the effect of 
lihedral upon the lift and drag coefficients in supersonic flow was investi- 
gated for thin airfoils symmetrical about their centre chord. The present 
paper is intended to complete the determination of the force coefficients 
by the calculation of the aerodynamic derivatives with respect to sideslip. 
the derivative of the pitching moment is obtained in Appendix I and is 
given in terms of the results of Part I. As in Part I, the equations of motion 


we linearized on the basis of small perturbations imposed on an otherwise 
iniform main stream. Wings having swept-back leading edges are included 
in the analysis and two cases arise depending on whether the leading edges 
lie inside or outside the Mach cone originating from their vertex. These 
two cases are considered separately. 


Previous work upon the derivatives in sideslip has been devoted to the 
case when the angle of dihedral is small. A first order theory was first 
obtained by Robinson and Hunter-Tod (2) using a uniform distribution 
of sources. The problem was later treated by Multhopp (3) who used both 
the cone-field theory and conformal mapping to advantage. Since then, 
an attempt has been made by Nocilla (4) to develop a theory giving values 


Quart. Journ. Mech. and Applied Math., Vol. IX, Pt. 4 (1956)] 
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of the pressure distribution on the wing at larger values of dihedral by ar 
expansion in powers of this angle. He does this for the case of a wing lying 
entirely within its vertex Mach cone, and in a later paper (5) he succeeds 
in obtaining values for the force derivatives on a delta wing which at ,) 
Mach number of 1-5 has a 60° sweepback. These values, which are fo; 
dihedral angles less than 30°, are in good agreement with the results give 
here. 


2. Notation and basic equations 

The angles of dihedral and sweepback are denoted by « and A respectively 
Cylindrical polar coordinates (r,@,z) are defined with the origin at th; 
vertex of the wing leading edges such that the z-axis is in the direction of 
undisturbed flow and the wing leading edges are given by z = rtan) 
6 =a, 7—a. The wing is assumed to be symmetrical about the plane 
6 


d 





We will consider small deviations of the wing from the above neutral 
position and investigate the forces acting upon it. The only deviation that 
‘an produce such forces is a combination of the following displacements 
Firstly, a displacement of the line of symmetry of the wing through a small 
angle in the plane @ }7 (an angle of incidence). This has been considered! 
in Part I. Secondly, a displacement of the line of symmetry through 3 
small angle in the plane @ = 0 (an angle of yaw). This will be considered | 
in the present paper. It has been sufficient to investigate each of thes 
displacements independently, for on a linear theory they may be super- 
imposed. Another consequence of this is that the case of yaw need only 
be examined for a wing which consists of two plane surfaces intersecting 
on the line of symmetry. If the angle of yaw is denoted by e, this line of| 
symmetry is given by r = ez, 0 = 0. | 

For the remainder of this section we refer briefly to the theory of Part! 
since the initial steps of the two boundary value problems are essentially 
the same. 


Let the velocity, pressure, density, and Mach number of the incident 


potential 4 is introduced such that the perturbation velocity has com- 


ponents of 
V; | ,,e¢ \" or (] 


‘ > ~ 
cr r 00 Oz 


y ped 


r 


On a linear theory, the potential function ¢ satisfies the equation 


D rp re 14. ma "hoo : Bd... 


this equation. Furthermore, as the problem is of the cone-field flow type 


flow be U, po, pp, and M respectively and write 6? = M?—1. A velocity | 
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e, the boundary conditions involve no fundamental length) it is possible 
to make the following transformation into the polar coordinates p, 6, 

2p 

1+ p? 

The Mach cone with vertex at the origin (henceforth called simply the 
Mach cone) becomes the unit circle in the (p, #)-plane, and inside this circle 
I f 
the wave equation transforms into Laplace’s equation with p, @ as polar 


ordinates 


3, The leading edges outside the Mach cone 
The perturbation of the flow by the wing will be confined to the region 
enclosed by the envelope of the Mach cones which emanate from each 
point of the wing leading edges. Upstream of this region the wing will have 
no influence upon the flow and the pressure will be the free stream pressure 
When the wing lies totally inside the Mach cone from the vertex, the 
region of perturbation will consist solely of the interior of this cone. When 
the wing leading edges lie outside the Mach cone the region will consist 
of the interior of this cone together with two wedge-shaped regions, each 
formed by two plane surfaces which are tangential to the cone and intersect 
long a leading edge. It was shown in Part I that due to the hyperbolic 
nature of the flow outside the Mach cone, the pressure is discontinuous 
cross such plane surfaces. Across the Mach cone the pressure is assumed 
to be continuous 
By either the first order Prandtl-Meyer expansion theory or the weak 


] 


shock theory, the pressure inside each of the wedge-shaped regions is found 


i », Ue ; 
to differ from the free stream pressure p, by +? > Secyssin a, where ¢ 
2 
denned 
sin ob tana. (3) 
p 
U2 
rm ) € ° 
Lhus 1f —Py Do Po — sec & sin «x, (4) 
p 
the pressure in those parts of the wedges above the wings 0 x, 9 = 7—a 
Is p, and 2) p, respectively. As in Part I, the wing divides the region of 
disturbance into two independent regions and it is only necessary to con- 


sider the upper region, the flow field below the wing being determined from 


the antisymmetry of the problem. The boundary condition on the surface 
of the wing is obtained from the condition of zero normal velocity. This 
implies : — . 

VV Vesin « (5) 


on all the wing surfaces. It is sufficient to apply this boundary condition 
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on the planes 6 = a, @ = 7—«a, although in yaw the wings are slightly dis. 
turbed from these positions. The error in so doing will be of the secon 
order in e. 

If the field is mapped onto the polar coordinate plane (R,#) by the 
transformation (2), the wing leading edges and the z-axis are represented 
by the three points £, F, and O (see Fig. 1) whose coordinates are respec. 
tively (cosec i, x), (cosecys,7—x«), and (0,0). The perturbed region then 


P 



























Fic. 1. The configuration in the (R, @)-plane. 


becomes the area enclosed by the wings OF and OF, the are BC of the| 
unit circle, and the tangents EB and FC to this circle. The pressure out- 
side the whole area is py, whilst the pressure inside those parts of the ares 
outside the unit circle, viz. BEA and CFD in Fig. 1, is p, and 2p,—p, 
respectively. If the pressure inside the unit circle is denoted by pg, then 
Po = 2py—p, on CD, py = py on BC, and py = p, on AB. If « is sufi 
ciently large, the Mach planes EB and FC may intersect and be reflected 
at the wing surfaces. This was considered in more detail in Part I, wher 
it was shown that the boundary conditions on the unit circle can b 
generalized to hold for all values of «. In terms of the non-dimensional 
pressure 

Ds P2—Po (6), 

Pi—Po 


the antisymmetry of the boundary conditions gives 


Oo+71 g Oo+Ye : 
" Cp. sint ; Cp. sint is 
lim P3 do — lim 3 de ; ( 
so J ob sin 7 so | 86 sin t 
—-O+Y1 —8+ v2 ) 
on the unit circle, where y, and y, are such that 
apy > a / _ » 
y 174 nod(z—2 
V1 5 r—wy 1 ( x) ) (8 
Vo Ia—aty mod(z7—2x) J 
, , a [ r— 2s 
and 7 is defined by T 5 : (9 
o|7 om 


The boundary conditions on the lines 6 = a, 6 = 7—« are obtained by 
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tly dis replacing in equation (12) of Part I, the value of Y% on the surface of the 
© second 


wings, given by (5). Thus 


Cp 


bv ¢] 0 ong X, 7—X. (10) 
* ty 
resente 
res} The region OABCD of Fig. 1 is mapped on to the upper half of the 
spec. | 
on then} {Plane by the transformations 
ie ; : 
Ze p—ian|(7 2a)( pet9)= 720) (11) 
, l i 
{= t+in= - s(a2+ ). (12) 
\s was shown in Part I, this transformation maps the points OA BCDO 
of the p,@ plane onto the points aD, i. cos7T!, +!cos7!, +1, and o 
n the real axis of the ¢-plane. 
We now define a complex function W(¢) by 
Wit) CP3 | i CPs (13) 
of the| on c€ 


ire out: and rewrite the boundary conditions in terms of W. Then on the real 





he are axis W is real for | 1 and wholly imaginary for || > 1. Equations 
Py] 7) suggest that W possesses simple poles at the points € + |\cosT|, 7 0, 
1, then whose residues are given by these equations. Moreover pg is bounded 
Is sufi- everywhere so that VJ O(¢) as f + 00 where t - 1. A function which 
‘flected’ satisfies all these conditions in the plane is 
wher . 
oe ' sin T l l 
ean W (C) _ ~ - +S — . (14) 
or 7 C“)?1¢-+cost 6 ¢—COST. 
S10! . ‘ , ° . eo, 
where the branch of (1—£?)! is chosen to be that which is real and positive 
a for |C l on the real axis. On the surface of the wing @ = a, 
Jeo 
‘P3 _ 0 
c 1) 
and equation (14) can be integrated to obtain the pressure p, at the surface. 
> The value of p. on the line of symmetry is zero because of the antisymmetry 
of the bound ir'y conditio1 Ss Hence 
2 sint 7 
p aresin! —— J (15) 
j 7 \(€2— cos?7)?} 
8 where the range of inverse sine is taken from —4z to +42. Equation (15) is 
valid f ind us, including those values at which shock reflection occurs. 


} 


logether with equations (4) and (6) it gives the pressure distribution on the 
suriace of the wing. 
When the angle of dihedral is small, the surface pressure distribution 


can be expanded as a power series in a. Taking the first term of this 
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expansion, we have the result previously obtained by Robinson and 
Hunter-Tod (2) 
Po— Do “Po _~ aresin | ali | 

Ba cos ys (z2?—r? tan?A)!} 

Another special case occurs when the leading edges lie in the surface of 
the Mach cone. The pressure distribution is found from equations (4), (6), 
and (15) by taking the limit as 7 > 0 
P2—Po = a» ll a : (17) 
B(ar—2a) (1— p*7k7—*0)) 
4. The leading edges inside the Mach cone 

The situation in the physical plane has already been outlined in the 
previous section. All disturbances are confined to the interior of the Mach 
cone upon whose surface the internal pressure p, takes the value py. Unlike 
the case considered in § 3, the whole region inside the Mach cone is inter- 
connected and the upper and lower surfaces of the wing must be considered 
together. The configuration of the wings is entirely symmetric about the 
plane 6 = $7, and from equation (5) the boundary conditions upon them 
are antisymmetrical. It follows that the dimensionless pressure 


Da P2 Po (18) 

Py U? 
will also be antisymmetrical about the plane 6 = $7. Thus p, = 0 on the 
lines 6 = 4z, 0 }7 and it is only necessary to consider one-half of 
the Mach cone, viz. that half given by —}a7 < @ < 4x and containing the 


wing 0 x. 

In the (p,@) plane given by equation (2) the configuration is identical 
to that obtained in § 6 of Part I. The z-axis becomes the origin (0, 0) and is 
designated E or B according to whether we consider the upper or lower 
surface of the wing at that point. The leading edge becomes the point 
A having coordinates (/,«) where 

9] 
= Beota (19 


and the half Mach cone (—3z7 < @ < 47) is transformed into the unit 


semicircle CBEDC (see Fig. 2). Considering the circuit A BCDEA, we 
denote the derivative along the circuit by ¢@/¢s and the derivative along 


the inward normal to the circuit by ¢/en, so that the boundary conditions 
for Pp, are 


Ps _( on BC, CD, DE 
cs (20) 


Cp, 


on 


) on EA, AB 


(16) | 





EFFE 


Consi 


It was 


ABCD 





A, B, 
) tively 
to be 
and é 


In 


Aew 
Wié 
Mor 


on 1 


- 


COTT 
sing 
as | 
bot 
sur 














EFFECT OF DIHEDRAL ON THE AERODYNAMIC DERIVATIVES 431 
(Consider the conformal transformation 


Aly? | (2 1) i 1) 1 -+-2a)/27r {(3° 1) (C2 1)! ya $1)" 2a)/27 


"| tty | oo ie | 
(21) 


t was shown in Part I that this maps the interior of the configuration 


{BCDEA onto the upper half of the € = €+in plane, with the points 





ir om, An 
\ . 
Pk. \ 
_— 
/ 
J 
ol a, ey 5 : 
a aii ic T T ——»- 
« j 4 e /) Dp E Q = 
S—plane 
ia. 2 
{1 B.C. D. BE, and A corresponding to 0, —y,; 1, 1, 5, and & respec- 
tively (Fig. 2 In this transformation the branch of ({?—1)! is chosen 
to be that which is real and positive for ¢ real and greater than 1, and y 
ind § are positive constants defined by 
) | 
i (22) 
Za = | 
] : | , 77 +2a)/277} | (§2 1)! § |(7—20)/277 (23) 
In the ¢-plane we define a complex function W(¢) such that 
Cp -Cp 
1] f Ps ? t - (24) 
€ 1) c& 
Rewriting conditions (20) in terms of W, it follows that on the real axis 
W(£) is real for ¢ y and ¢ 6, and totally imaginary for —y < & < 6. 
Moreover, as in Part I, it can be shown that W(C) possesses no singularities 
on the re ixis at EH and B. At the wing leading edge, however, which 


corresponds in the ¢-plane to the point at infinity, a comparison with 
singulariti incurred in analogous problems suggests that Wm~ A C4 B 
as C00. This assumption is reinforced by a more detailed study of the 
boundary conditions at the wing leading edge. On both upper and lower 


suriaces of the wing, equation (5) shows that Vj Uesina, so that we 
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may use the theory of Part I (Appendix) to assert 
k 
tt cae . 
re) lim | Wd¢ 0, 


ko i: } a 


where the integral is taken around a large semicircle in the ¢-plane. 4 
function which satisfies all these conditions and is integrable at all finit; 
points of the upper half-plane is 


Wi) = 


—, (2 


um 
SY 

— 

—| 


’ 


u~ 


where A is a real constant and the branch of (f? L-yl—Sy)* is choser 
to be that which is real and positive for € > 5, n = 0. 

As @p,/€7 is zero on the wing surface, equation (26) may be integrate 
to give the pressure on the upper wing surface, 





QS \1 


6)'(€+y)!. (27 


uty 


DP, = 2A( 


There is no constant of integration since the pressure p, is zero at ¢ i 
The constant A is evaluated from the boundary conditions on the line o 
symmetry. An examination of the components of the perturbation velocity 
near this line reveals that 1) > —Ue as r > 0 in the plane @ Ler. Since 


¥ is zero on the Mach cone, equation (12) of Part I may be integrated t 
give , , 
_R ~1l—p? ep 


4 6 
dp. 8 
2p” 26 p 


« 


0 
Z : ' - , , sn OD, 
Upon transforming to the ¢-plane and substituting for on the rea 
: p cl 
axis we have 


») 


~ 
> 
~ 
no 
bo 
ir, 
LY 
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A Peete OE (x x @), (29 
I A p \to—s)y+s)*| - 


Lee) 


where the asterisk denotes the finite part of the integral as defined bi 
Hadamard (6). This is made necessary by the divergence of the integra 
when « <0. A detailed justification of equation (29) for this case i? 
given in Appendix II. 

When the angle of dihedral is small, the pressure distribution on the 
wing can be obtained explicitly in terms of the flow parameters and thi 
physical coordinates. To the first order equations (22) and (23) giv 
y= 5 (1+-/?)/2/ so that the transformation (21) can be simplified on thi 
wing surface to 


1 (1+/*)(1—p?) 
2 (l?— p?)4(1 p*l2)*’ 
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and equation (27) becomes 


= — (31) 


L (2 cot®A—r?)" 


To evaluate the constant A from equation (29) we need only consider the 


behaviour of the integral in the neighbourhood of the endpoint € = 8. At 


this point the integrand behaves as K(8—€)-“*-™/} for some constant K, so 


that to the first order deal? 


A oo 
Bar(1 —I*) 


. (32) 
Substituting for A in equation (31), the actual pressure on the wing surface 
is found to be , i 
2eap, U? r = 
os (33) 


Pe = Po tan 


:result which was previously found for this case by Robinson and Hunter- 
Tod (2). 

When the leading edge approaches the Mach cone we may write 1 = 1—t 
where ¢ is small. The parameters y and 6 are then given approximately by 


l at \* ‘ a art 2 : 
Y me =} ? 0 1+ 5(- 53) ° (34) 


| —f 27] 327 /(77 2a) 12 —_ 
é 1 gala on 6 = in, (35) 
2 \7— 2a 1 -+- per —<0) . 
| at 2 1+ 9277 7 —2x)]2 — 
é 14 | abl on 6 = a. (36) 
2 \7— 2a | — pom —20) 
To the first order in ¢t, equation (29) can be integrated by ordinary methods 
to give the constant A (c sin a)/7St. Substituting this in the expression 


27) for p, and taking the limit as ¢ > 0 the pressure distribution when the 
leading edge lies on the surface of the Mach cone is given by 
fep, U?sina = p7™7—2a) 


de Do : 
Be : 6(a7— 2a) 


] 27/(7—2a)* (37) 
This agrees with the result (17), obtained from the limit as the leading 
edge approaches the Mach cone from the outside. Thus all the aerodynamic 
derivatives in sideslip are continuous as the leading edges pass through the 
Mac 

Mach cone 


5. The aerodynamic derivatives in sideslip 

The derivatives are derived for a wing of delta planform. Thus if the 
maximum chord is ¢, the trailing edge is given by z = c, 0 = a, 7—«, and 
the semi-span and area of the w ing surface are ctanA and c* tanA respec- 
tively, 


092 F f 
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Let (x,y,z) be a left-handed system of cartesian coordinates such that} yhere 
the z-axis coincides with the z-axis of the polar coordinate system already 
defined and the x-axis is given by 6 = z = 0. We resolve the force and 
couple acting on the wings into force components X, Y, and Z along the 
axes x, y, zrespectively and moments L, M, and N about these axes. From 
a consideration of the antisymmetry condition it follows that 


Y=Z L 0, 








and 
and we define the non-dimensional derivatives m,, n,., x, as follows: Whe 
y , . , ass 2M tan?A to obt 
Non-dimensional yawing derivative m, poem 
, Po l *ec® 
y 9 i, 
: : — 2N tan?A . . 
Non-dimensional rolling derivative n, ——_____ (38 
) | 26¢3 
Po 
; : F : : , . 2X tana 
Non-dimensional sideslip derivative 2, ——> 
Po l “eC” 


(The yawing derivative is referred to the origin, i.e. the vertex of the wing | where 


leading edges.) Be 

When the leading edges lie outside the Mach cone, the pressure on the | | 
upper surface of the wing 6 = « may be considered in two parts: firstly, : 
the pressure p, on the whole wing, secondly, the pressure p,—p, on that Beot 


part of the wing inside the Mach cone. (This does not hold at values of » 
for which reflection of the Mach wedges occurs, i.e. «a > }a7+4u: there} 
will then be another term from outside the Mach cone.) The contribution | 

to X from the upper surface is thus 








¢ § B c zeota 

2 sin | | | (P.—p,) drdz + | | P; ards}. (39) } and b 
0 0 0 0 
In 
One integration may be completed by a substitution into the cone-field | j9° 4, 
coordinate R = fr/z. Furthermore, p, and p, are given by equations (4) and} gai, 
(6) so that with similar results for WV and N the derivatives are given bj , 
° : 6. M 
mM, 28x, sin ps | | M: 
4 sin2us sin2« . . wee ef 
a= a ” — | 2(cosee ys —l)- | P3(x, us, R) dR + | p3(—«a, us, R) aa | Mack 
0 0 | 

(40 its ti] 
the ¢ 
4 sin? sin « . ; | cale 
= — | cot’y-+ | ps(a.b, R)RAR + | ps(—a,, R)RAR| p 
3 cos us ° * aaa } | whic! 
(41)) Fort 






























EFFECT OF DIHEDRAL ON THE AERODYNAMIC DERIVATIVES 435 








th that | where 1 1 
ait i fsint é 
uroad p,dk | | = A : dp, (42) 
ce and | 7 — 20 (1+ p*)(€°—cos*r) 
a 0 0 
ne the 
From § tsinz [ Ep 
| psRdR : -_ — dp, (43) 
: m—2a J (1+ p*)*(€°—cos*r) 
0 
nd é 4(p” 7 — 20) _| pa" 2a 2. (44) 
When the leading edges lie inside the Mach cone, the results are simpler 
obtain, and integrating over a similar wing surface we have 
, cota Scota 
: fsinatan-A . : . ’ 
jv, tana = | | pyla,l, R)dR | py(—a,l, R) dR], 
g pia L 9 0 
(45) 
3 cot A 
t{tan-A ; - . : 
n, = | p,y(a,l,R)RdR | py(—a,l,R)RdR|, (46) 
p> | J 
¢ wing where 
ta 1 
P ; 8(a7+-2a)bA(1—b?c?) (1 »*) da a 
on teh =| dR —ar | le en 
— : 7(1—b | —c*) J 2(1-+ p*)?(1—a?b?)?(a*—c?)! 
" :' 
os at ) F &( 77 Za 1(] hc") r *(] »*) da 
Sotap | py RdR : f C———__, (48) 
4] 7(1—b2)(1—e2)! | (1+ p2)8(1—a2b2)i(a2@—c2) 
‘ | l|—al “71 a—c |(%—-2aen 
| | (49) 
L1+-al la+c] 
= Be. ite 
} nd b and ¢ are constants such that Y 59 O : 
b? | —c? 
In Figs. 3 and 4, values of Bn, and Bx, are given for dihedral angles of 
lO’ and 30°, at varying values of (tanA)/8. The yawing derivative can be 
} 


isily obtained from the values of Bx. 


6. Modified derivatives when the leading edges are outside the 


Mach cone 
7 Ifthe wing has a finite span, the pressure on the wing surface outside the 
Mach cone is no longer constant but depends on the shape of the wing at 
{))| Itstip. However, apart from a factor sin «, this pressure is independent of 
the dihedral angle. Thus the force derivatives on such a wing can be 
dR calculated from the values at zero dihedral by superimposing those forces 


which are due to the modifying effect of dihedral inside the Mach cone. 
For this purpose modified derivatives are defined which are based only on 
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Fic. 3. The variation of the rolling derivative with sweep-back for dihedrals 
of 10° and 30°. 
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Fic. 4. The variation of the sideslip derivative with sweep-back for dihedrals 
of 10° and 30°. 


the forces produced by the modified pressure (e.g. on the upper surface o 
the wing @ = a, the modified pressure inside the Mach cone is p.—?;): 
If X’, Y’, Z’and L’, M’, N’ are the force and moment components generated 
by the modified pressures, the modified derivatives of sideslip force and 
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EFFECT OF 


yvawing and rolling moments are defined as 



































; 2X’B ; 2M'B? , 2N’p? - 
v, 3) Mm, é. =» é i ll (50) 
Py l “eCc* Po l 2e¢c3 Po l 2ec3 
. 
a 30 40 50° 
Fic. 5. The variation of the modified rolling derivative with dihedral. 


respectively (from the antisymmetry condition Y’ Z' L' = 0). Hence 


by an integration over the wing surface inside the Mach cone, 


; 1 
BBx, ; sec ys sin® | 2 ( p3(a,%, R)dR — [ pz(—a, b, R) ar (51) 
0 0 

2 
an, fsecu sina] | ( P3(a, Jb, R)RAR | Ps x,%, R)R ar|, (52) 
0 0 
where the integrals are given by equations (42) and (43). The variations 
f Bx, mi, and Bn’. with dihedral are given in Figs. 5 and 6 for 4 = 0°, 20°, 


40°. 60°. SO 


7. Conclusion 
When the wing leading edges lie inside the Mach cone and 
ScotaA 0-645497, 


solutions for the sideslip, rolling, and yawing derivatives have been given 
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by Nocilla (5). Expanding the derivatives as power series in «x, he obtains 
values from the first two terms that agree very well for « < 30° with the 


results given here. 
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Fic. 6. The variation of the modified yawing and sideslip derivatives 
with dihedral. 


At dihedral angles of less than 50°, the results suggest that for a given 


Mach number and sweepback, the rolling derivative is proportional t 


sina and the sideslip and yawing derivatives are proportional to sin*s.! 


Clearly, such a formula would enable the derivatives at any dihedral t 


be easily obtained as the results for small dihedral are already know 


explicitly (e.g. see 2). The greatest deviation from this rule occurs whe 
the leading edges lie inside and close to the surface of the Mach cone 
At this point the true derivatives for a dihedral of 50° are approximately 
10 per cent. less than the values given by the above rule. For dihedra 


angles less than 30° the error is less than 3 per cent. 


In conclusion I would like to thank Dr. W. Chester for suggesting th 


subject of this paper and for his help in the preparation of it. 
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EFFECT 


APPENDIX I 


The pitching moment derivative 


In Part I and the above paper all the derivatives except that due to the pitching 
ment |} been determined for a delta wing. This last derivative is considered 
Jecau f the antisymmetry conditions, no perturbations of the kind considered 
regoing pages produce any contribution to the pitching moment. The only 
» the moment mes from a displacement of the line of symmetry of 
wing through an angle e’ in the plane 6 har; such a displacement was investi- 
lin Part I. Denoting by Z the pitching moment for a delta wing of maximum 
rdc,t verodynamic pitching derivative 1, (with respect to the vertex) is defined 
2Ltan?A 
i. ———= (53) 
Po l "ec 
An integration over the surfa f the wing shows that for a delta wing 
) 
=—~C,tand, (54) 
se 
re ( the lift coefficient defined in Part I. Values of C; are adequately given 


APPENDIX II 


Consider the circuit DPQ in the (p, @) plane, where QP is the are p = h,a < @ < 4a, 
DP he line @ wT, | p h. Ifh is small the are QP is represented in the 
eircuit I; g, 2) where 
4 6+ ke O(k?) (55) 
1 fi the transformation (21), 6, and k are such that 
7 — 2a = 
A x — 6;, (56) 
i les (57) 
j - - voi 
| | yo+ | 2(07— 1) 
In t] plane the line DP represented by that part of the real axis for which 
. This lin lenoted by T, so that I, +T, refers to the whole circuit 
DPO 
Now ler the follow ntegral, whose path of integration in the ¢-plane is 
* 1—p? . = 
| W d@. (58) 
Pp , 
I 
S quation (24) and rewriting the integral in the (p,@)-plane, 
a | { 1 OMe « Cp Cp cp " 
W di | Et dp — ga + 6 i +i Ys d0|. (59) 
{ p cp cp ob 
On t C.@ col t and p, Is zero, so that by substituting from equation 
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where [ ]r, here indicates the change in value upon completing the circuit T,. Ina 
like manner equation (59) can be simplified for the circuit I, by a substitution for 
Op,/ep by equation (11) of Part I: 


1 l p® ’ » 2 , ’ i ' 2 , 9 > 
| p ili a | tea jiPalt. + BG Volrs + O(h*). (61) 
I; PQ 
But from the boundary conditions on the wing surface V, ~ 


Uecos@ as p > 080 
that equations (60) and (61) together give 


. 


l—p* ... .. Qesina 2e 7% a —- 
| ; - W d¢ 3 —~ Bt n Pay or, t BG ely, 41, +O). (62) 
Tn+Ts 


At the Mach cone boundary the perturbation velocity Vj is zero and on the wing 
the condition of zero normal velocity gives yj = —Uesina. py, is zero on the Mach 
cone and its value on the wing surface is given by equation (27). 


W is given by 
equation (26), so that 





- l—p? i(2¢+y—8) ,, 2e  2Aik*(y+8)! 
A ea EE Hein SA Pen 1 O(h). (63) 
F p (C—0)2(C-+ y)? B h 
r+ 
Taking the limit as h > 0, it follows at once that 
5 

2e *f ] 2 fF —§ +4 

. -A f 2a Ye dé (x #0). (64) 
p , p (o—€)*(y+€)? 
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OF SOURCE-FLOW FROM 
AEROFOILS 


By L. C. WOODS (University of Sydney, Australia)t+ 
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SUMMARY 
An exact theoretical treatment of the influence of sources of fluid on and within 
two-dimensional aerofoils is developed. In the case of a source placed within an 
foil the fluid is ejected through a duct, the shape of which is taken into account 
he theory. The theory is given for incompressible flow but it could be extended 
leal wit l arized subsonic compressible flow. Formulae for the thrust, lift, 
1 moment are derived fr slasius’s theorem. The case of the ‘directed’ source 
s attention, and tl ience of the angle of ejection is determined. 


List of symbols 


], 0 
1 
yy 
} 
() 
‘ 
( 
: p 
} b, 
ae 
( 7 | 
1» 
A 
[Quart 


Journ 





velocity vector in polar coordinates. 
undisturbed velocity magnitude (at infinity). 
lc | l l q) 


f YT) 
QO ly 


26. 

the complex stream function. 

a subscript to denote values relating to a particular fictitious 
flow (not satisfying Joukowski’s condition). 

analytic functions of w,, defined by equations (1), (4), and (6). 

the real and imaginary paris of €. 

aerofoil chord length. 

incidence measured from the no-lift attitude. 

2a < do < 2a, py = 0. 


the real physical plane. 


the aerofoil lies in 


the source strength. 
(Uc, the mass coefficient of the source. 

convenient functions of the source strength, defined by 

(3), (5), and (12). 


defined by equation (18). 


equations (2), 
coetticients 
density. 
drag, lift, and moment forces. 
drag, lift, 


distance mea 


and moment coefficients. 

sured around the aerofoil surface. 

distance of source from trailing edge. 

as subscripts to denote quantities relating to the real non- 


circulating flow and the real circulating flow respectively. 


the trailing edge is at y a—B. 
values of y defining the positions of the stagnation points. 
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1. Introduction 

Source and sinks on aerofoils have long been considered as possible means 
of controlling boundary layer separation and the circulation. The mathe. 
matical theory for the influence of sinks on the lift of aerofoils has been 
given in various forms by a number of authors (see, for example, Ehlers, 1), 


but to the author’s knowledge existing accounts ignore the effect of duet 





: 











The _z-plane 


Fic. 1. 


shape, and usually consider only the simple case of a sink right on the 
aerofoil surface. 

In this paper we present an exact theory for sources within aerofoils 
the fluid from which is conveyed by a duct of general shape and emitted 
into the main stream. The special cases of (a) a source right on the aerofoi 
surface, and (b) a parallel walled duct ejecting fluid at a given angle into th 
main stream, receive detailed attention. 

The fluid is supposed to be inviscid and incompressible. Also the source 
fluid is assumed to be at the same total head as the main stream fluid, so 
that both the pressure and the velocity is continuous across the strean- 


lines separating the source fluid from the main stream. 


2. The basic conformal transformations 

Consider the flow shown in Fig. 1. There are three stagnation points on 
the surface, namely those at b, D, and D’, the positions of which are fixed 
as follows. First the Joukowski condition at the trailing edge F is satisfied 
by movement of D into coincidence with F, and then the positions of D 
and B are uniquely determined by the strengths of the source and thi 
circulation. 

The first transformation we introduce can be regarded as a fictitious non- 
circulating flow about the aerofoil such that the points D and D’ are at the 
same velocity potential. This flow plane is termed the w,-plane, wher 
Wy = bo +t, and is shown in Fig. 2. The aerofoil is assumed to lie in the 
interval —2a < ¢, < 2a inthis plane. The actual relation between 4a and 
the aerofoil chord length in the z-plane ¢ will be determined later. The 
source output is Q units per second, so we can assign ys, — 0 on D’G", and 
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Q on DG,,. The leading edge streamline is a line of discontinuity of 


but this causes no confusion in the following 
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The basic conformal transformations 
Fic. 2. 
\lso shown in Fig. 2 are three further planes derived by conformal 
transformations from the w,-plane. These mappings are readily calculated 
by the Schwarz—Christoffel mapping theorem. First the wp -plane is 


mapped into the upper half of the ¢-plane by 
ss () 
1a P(1—t?) +- — Int—2a, (1) 


Q e 
#2 In(#2)’ 


Sat | t 
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_ 1—(Q/47a)In ty 





where r _— ; 3 
1—# (3) 
and the points D’, D map into t = t), t = —t , respectively. 
The ¢-plane is mapped into the semi-infinite strip —o <7 <0, 
0<y < 27 in the ¢ (= 7++7y) plane by 
t —cosh 3€. (4) 


The points D’, D map into y = 7+6, and y = z—6, on 7 = 0, so that 
from (4 as 
(4) t, = sin $6. (5) 
Finally the o-plane shown in Fig. 2 is derived from the ¢-plane by 
G 2(1— 2#?). (6) 
Most of the detailed calculations will be made in the ¢-plane, but the 
o-plane has been introduced for the purpose of applying Blasius’s theorem 
to the problem. The symmetry of the mappings means that in this plane 
the characteristics of the real flow are continuous across A,, B, and hence 
closed contours enclosing the aerofoil surface will be closed contours 
enclosing the strip —2 < o < 2 in the o-plane. 
The following equations, which are readily derived from (1) to (6), will 
be required later: 
| a ° O r 


——- 2% = 1—-cos6, (7) 


or if Cg is the mass coefficient of the source defined by 


Co = Q/cU, (8) 

where c is the chord length, then 
Co 2n(-)sin* 38. (9) 

cl 7, 
Also es — Af | | + Oo-4)\, (10) 
a o | 
and My _ gp} 4.404 740.4 O(9-3)| (11) 
da | o o | 
Q i ie 

where 5 4sin*}6. (12) 


3. The general solution in the ¢-plane 
Consider any circulating flow past the aerofoil. Let (q,@) be the velocity 
vector in polar coordinates, and take lim(q, #) (U,«), then if w is the 


complex stream function, 


f in(U =) in )+i0, 
‘ dw q 
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where 


Y OF 


In the ¢-plane, although y 
leading edge streamline, the flow conditions at a given value of 7 on these 


ines are identical. Now / is an analytic function of z, and since z : 


0 and y 
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lim i (14) 


ia. 
27 no longer correspond to the 


z(f), 


fis also an analytic function of ¢. The function satisfying these conditions 


is (see Woods. 2) 


where 6() 


15) yields 


in( | 


. 


. 
0 


) is the value of @ on 7 


A(y*)cot (y* ! if) dy*, (15) 


0 at y = y*. On the aerofoil surface 


| O(y*)cot 4(y* —y) dy*, 


« 


on 





or In| In|sin }(y*—y)| d6(y*), (16) 
q T 
0 
on integrating by parts. 
At large negative values of 7 equation (15) yields the expansion 
f = by+b, e& +b, 5+ O(e3s), (17) 
where b HAy* je" dy*, n 1,2 
oe (18) 
by : O(4 ) dy 
2n 
Now © implies o ©, and from (10), 7 -~ —oo; thus we deduce from 
14), (17), and (18) that 
O(y*) dy* X. (19) 
9 
Equations (10), (13), (17) to (19) now give the expansion 
dw L 9 l « } 
Ue-ia) 471 — (b,—4b}3) + Ole 3)| (20) 
dz | o a | 


As a particular value of the theory consider the flow characterized by W». 


We denote the velocity vector in this case by (qo, @,) and set 


» Mg). 


lim (qo, 9) = (U 
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Then corresponding to the above equations we have 


In 


U lf ; , 
In| = — | In|sin $(y*—y)| d0(y*), (21) 
Yo 7 
0 
lw ed l 
GW _ Ueivl 1 4 %r0 - (op — 40%) +, + O(o-%)|, (22) 
dz | oO = i o } 
bro res Oo(y*)e"™ dy*, (23) 
0 
and XM => | Oy(y*) dy*. (24) 
0 


4. The closure conditions 
Let C be any contour enclosing the aerofoil, then the aerofoil is closed if 
| dz 0. In the o-plane this condition can be written 


* ({ dz\ (dw, 
| (7 )( 7 heals 


Hence from (11) and (22) 


aP ¢ b || a 
ido 10 L Olo-2 eb Je 
— | l (o-*) 1 


ox oO 


L O(o-2)\ do 0, 
C | 
which requires that a) = b,) from the residue theorem. 
It now follows from (12) and (23) that 


or LZ wae , . ” 

A 4sin*$6 | 6y(y*)sin y* dy*, (25) 
0 

l c 7 ke a i9e 

and 0) = | Ay (y™ Jee s y* dy*, (26) 
0 


which are the required conditions. 


5. The application of Blasius’s theorem 
In a general circulating flow the drag D, lift L, and nose-uwp moment 
(about the origin in the z-plane) .V, are given by 


e-ta1(D—7iL) lip | (ae) a (27) 
old ,. ft. dw\? - 
and M—iN bo (7) dz, 23) 
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fow is at an angle «a to the x-axis. 
2] To find VW we need the expansion for z in terms of co. Now 
dz be \(S 
da dwy doa 
dz aP ° l 
erefore ( 1 + (b,,>—4a?2 + 2a,) — + O(o-*) (29) 
| do | - - Oo 
a n using (11), (22) and b,, = a). Hence 
aP } 9 ‘ 9 
¢ y) | (Ds las + 2d) - Olo 3), (50) 
{ | 2 7 ge” | 
24 
| where the constant of integration has been made zero. This is equivalent 
| to fixing the origin of the z-plane at a certain point, termed the ‘centre’ of 
the profil 
S¢ Lil cae : , 
lo find this centre consider the surface integral 
| [ = | zdy 
0 
a 
We can write / i ~ : do, 
J do 
hen from (10). (30) and the residue theorem we find that J 0. Hence the 
| centre of the profile is calculated from 
ly 0 | y dy 0. (31) 
0 0 
| when x. y have been determined as functions of Y: 
25) | We can now evaluate (27) and (28) by the residue theorem. Equation 
27) ean be written 
: * (dw\2 {dz 
D—iL Lipe!® | } | do. 
6 * J \dz) do 
tf 
ind henee from (20), (29) and the residue theorem 
D L 2rzapPUe' . "b,. (32 
Similar from (20) and (238) to (30), 
ies M—iN 2rrip(a P)2e2(%-(b, — 3). (33) 
\ more explicit result for the drag can be obtained as follows. From 
7 ?0) ind () I 
- dau } aP \ 
( 1] 4 1 git + O(z-2)!, 
d | Uz | 
28 bh,aP 
so tha uw ES, ul 1__ g%]n(z) + O(z | 
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‘term being introduced to allow for the fact that the undisturbed 
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Recalling that the real part of the coefficient of In z in this expansion equal 
the strength of the source (or sink) at infinity divided by 27, we deduce 


Q = 2naP refeX~-%,}, (34 
or, from (12), a, = 4s8in?45 = refet%-%H,}. (35 
But, from (32), D = —2mnaPpU refeX%-%,}, 
and hence D = —pUQ. 


In terms of the drag and mass coefficients defined by Cp, = D/}pcl™, and 
equation (8), this reduces to the well-known result 
Y BY Al 20 
Cp a. Q: (36 
Equations (32) and (35) enable us to rewrite the expression for the 
moment in a more convenient form. From (33) 
M 2p(aP)? im{e?"%-%h,' — 2ap(aP)? imfeX%~-%H, }*; 
but, from (32) and (35), 
41 sin*46 


im{e%% ~b,}? - Zim{e%% Ib, }refel(%-h,} = ——_* , 
mpaP 
and hence 
ah a 4aP' acate - 
M 2rp(aP)? imf{e2(%-%H,} > 7) Losin’ ys. (37 
: 2 


6. Real non-circulating flow past the aerofoil 

We now consider the flow past the aerofoil which satisfies Joukowski’s 
condition, and is such that no lift force acts on the aerofoil. Let 6, denote 
the value of @ on the aerofoil surface in this flow, and let lim @ = «,, so that 


Z x 


a, is the no-lift angle. Referring to Fig. 1 we notice that in this flow the 
rear stagnation point must move to coincidence with F, while the other 
stagnation points must move to new positions to maintain zero circulation 
and unchanging source strength. The difference (@,—64,) is thus due only 
to the shift in the positions of these stagnation points. 

Let the trailing edge be at y = 7—f, and let the stagnation points 
originally at D’ (y = 7+6) and B’ (y = 27) move to y = 7+6,, and 
y = 27—A,, respectively, then 

- B a A a es 5) 


(T-1-0, << 7] 5) 38) 


O, - Go 7 
lice A, <y < 2n) 


From (19) and (24), Xy— Xp | (0,—0,) dy 


hence 
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from (38). In this equation a, and £8 can be regarded as known, a, following 
from (24), and £, which defines the position of the trailing edge, being 
letermined as follows. 


Let s be distance measured from B along the aerofoil surface, then 


, [ de 
. Yo 
0 
* U [dd,\ [c 
Therefore o] l | ldy dv. 
J Yo dt dy ; 


0 
On the aerofoil surface (1) and (4) yield 


dd ) 
z Ge PS 1. =, t cos dy. 
dt at = 


These relations and equation (12) enable us to deduce that 


| faP P U . 
s 5 | 7 (cos y+-cos 0)tan ly dy, (40) 
#\< / Yo 


the value of Ug, in this being determined from (21). Suppose now that the 
urce is placed at a distance / from the trailing edge, then it follows from 


l | ja P g U . 
: af 7 “heed +-cos 6)tan dy dy. (41) 


It remains to calculate (4aP/cU’). Let the total length of the aerofoil 

perimeter be expressed as 2kc, then k can be calculated from the known 
erofoil shape; it will be approximately unity for thin aerofoils. Then 
irom (40 on 

| of J (cos y+-cos 6)tan dy dy. (42) 
ta P tk |} do 


lhe parameter 8 can now be calculated from equations (41) and (42). 


In addition to (39) there are two further relations between «,, A,, and 4,, 


vhich determine them uniquely. First we have equation (35), which 


expresses the invariance of the source strength. From (18), (25), (26), and 
5) we find that for the flow considered in this section 
2r(] os 0)(1 COS| Xp 
COS(a%)—a,) | (8,—%)sin y dy—sin(ap—a,) | (6,—4)cos y dy, 
0 0 


Gg 
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which reduces to by (9 
2(1—cos5) = cos(ap—a,)+ C08(a%y— a, +A,)—C08(8, — avy + a4) — repre 
silat depe: 
cos(B+ a )— 1), P 
4 7 ; - ; =he ow Fr 
by (38). Equation (39) enables us to rewrite this in the concise form 
sin?45 = sin $f sin $6, cos $A). (8)) oy 
Secondly the zero lift condition is found from (18), (25), (26), and (32) t 
require that ae 
27 sin(ay)—,)(1—cos8)-+sin(vy—a,) | (0,—4)sin y dy 
0 
+- COS(a)— a) (0,—@,))cos y dy = 0, 
0 
which by (38) and (39) reduces to 
sin'}(A, +8,)cos 48 = sin 48 cos 4(3,—A,). (44){ The 
. . ‘ . . .* can | 
Equations (39), (43), and (44) determine a,, A,, and 6,, so the flow is 
uniquely defined 0 
u 
7. The lifting aerofoil 
Next we consider the case of a real circulating flow. The rear stagnatior 
point remains fixed at the trailing edge, but the other stagnation points| 
move with increasing incidence. Let 6, be the value of @ on the aerofoi 
surface in this circulating flow, and let lim @ = a, so that « is the angle of| 
incidence. Suppose the stagnation points originally at y = 7+-8, y = 7-3 
and y = 27 are now at y = 7+56,, y = 7—B (Joukowski’s' condition), and\ ‘I 
; ~ wher 
y = 27—Ag, respectively, then ‘ 
‘ ? tne © 
a ee ae 
be p i 0) possi 
0,—9, 7{ (7-+05 < y < 7+0) (45 
ce Ag < y < 22) 18. S 
0 (elsewhere). Th 
The calculation now proceeds just as in section 6, Corresponding to (39 and ¢ 
43), and (44) we find : shap 
(45) (44) 2(a— ay) B L-A.—dz, (46 : 
eee Pua Te _| to th 
sin*46 = sin $8sin $6, cos $Ag, (4i)7 
2 2/ 2°2 as influc 
and L = 4rpaP U{sin $(A,+6,)cos }8—sin $8 cos $(5, —Ag)}. conv 
By equations (46) and (47) we can eliminate A, and 5, from the equatiol) flow: 
for L. Then replacing L by the lift coefficient C,;, = L/3pcU*, we find that | give 
:' 4aP\.. Pana As 
C7 277| —— |{sin(a—a,—fB)-+ 2 cot 48 sin?48!, (45 
, cl I 2] an ere 
, on th 
: 4aP\ . ; ; 
or Cr 2a sin(a—a)—B)+2Cy cot $8, (49) the s 
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by (9). The second term on the right-hand side of (49) does not necessarily 
represent all the effect of the source on the lift as both a) and B may also 
lepend on Co. 


1 
1) . ee 1 : : : 
From (18) and (37) it follows that the moment is given by 








orm 
4 » ) | . ») { h ») 
M 2p(al “sin 2 (ap } (6,—6,)sin 2y dy 
| (32)4 | 0 
COS 2(a%)— a) #,—6,)cos 2y dy 
0 
27 27 
. sin 2 x) | 6, sin 2y dy COS 2(a)—«) | 4, cos 2y dy\- 
0 0 
14] The integrals involving (#@,—4,) are evaluated from (45), then A,, 5,, and L 
g in be eliminated from the result by (46), (47), and (48). It is found that 
row ° 
7 (4a.P\*| ) ) Baa 5 1A\cin215 
| ; cos 20 sin 2(a Xo) 16sin 3P COS(A— Xp -18)sin®36 
t\ ¢ | i 7 % 
ve nacion |} 6 cot $8 sin*}6 sin 2(a)—«) | 6, sin 2y dy 
n py 7 . 
0 
ngle off ] : P ‘ , 
COS 2(a9 x) | 6, cos zy dy » (00) 
- p 
‘ vhere C\, is the moment coefficient M/3pc?U?. This equation applies to 
} the moment taken about the origin of the z-plane. Further reduction is 
possible only when the shape of the aerofoil and duct is known. 
45) 
8. Sources on thin aerofoils 
| The theory so far developed is exact, regardless of the source strength 
- to (39),| and of the aerofoil thickness. To make further progress we need to know the 
, ' 
pe of the aerofoil and duct. For simplicity we will confine our attention 
‘ the flat plate aerofoil, although the results obtained relating to the 
luence of the source on the aerofoil characteristics will be applicable to 
nventional aerofoil shapes with but little error. The effects of the source 
equatiol) flow found below could be linearly superimposed on the thickness effects to 
id that | give a reasonably close result for conventional aerofoil shapes. 
\s far as the duct is concerned only two simple cases will be considered 
+* rr Z . . . 
ere. These are (a) no duct at all, that is the case of a source placed directly 
n the aerofoil surface, and (b) a parallel walled duct of such a width that 
} 


he stagnation points occul the corners of the duct exit (see Fig. 3). If 
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the duct width at exit is too small, or too large, the flow will be required to 
turn the sharp exit corners, a physically unrealistic case that we wish to 
avoid. 





Physically_realistic source flow 





Fic. 3. 


8.1. Source on the aerofoil surface 














Fig. 4 shows the fictitious flow introduced in section 2. Assuming that 
the leading edge C is at y = 27—A we have the following table of values 
Point | B F D EB Db’ S B’ 
‘: | 0 a—B a—d 7 a+é8 2a—dA 2r 


Hence from Fig. 4 it is apparent that 
(‘ B<y< 7-8) 
7 (7 - 


y : 
| (on A<y < 27) 
0 (elsewhere). 
oo 





Vp 
sate mF 


D\r ey 
WH 


A em a 


Source on a flat plate 





Fig. 4. 
With these values of 6, equations (24), (25), and (26) yield 2a)+A+/ 
cos 8 = cosA, and sin8 = sind, respectively. Thus 


A B, (5. 
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“ = —B. (53) 


The next step is to find the relation between the perimeter distance s and 


y from equations (16), (40), and (42). With the step function defined in 


51) the Stieltjes integral in (16) degenerates to give 


and therefore 


, jsin 5 | 


a 
In( ae 


Yo = 
ii 


1 (2 —S—y)sin 3(7+6—y)sin 4(27—y)} 


Cc 


7 —B—+)sin }(7—y)sin }(27—A—y)| 





er 
os (0+ y)cos 4(0—y) ~ 

as : “tan Fy, (54) 
sin(b--y) 


n making use of (52). With this value of q,/U equation (40) gives 


s being measured from 


that is, from 


If Sis the distance me 


56) give 


where « 


| 


») 


Bb. 


52) and (55), 


asured from the leading edge (y 


ta P " 
= }{eos 8 cos(y+-)}, (55) 
The chord length is ¢ = s(7—B)—s(27—A), 
ta P (56) 
Cc ° v0 

rp 


—f) then (55) and 


S le(1—cose), (57) 


3 is the angular displacement measured from the leading 


> 


edge. Finally, the source is at a distance s(7—f8)—<s(z) from the trailing 


edge, that is. from (55) 


| 


and (56), 


l }c¢(1—cos 8). (58) 


The lift coefticient now follows from (49), (53), and (56); it is 


( j 


ae on 0 " 
27 sin a+ 2Ug cot LB, (59) 


which is identical with that for a sink on the upper surface ate = B(y = 9). 
The moment coefficient is easily deduced from (9), (50), (51), (53), and 


56). The result is 


where from (31) and (55 


hord point 


The centre 


of 


pressure 


. Y : Q . { 
2C, sin $8 cos(a+ 4; 


© 


(2 
3) ——“ cot $8, (60) 


this coefficient is for the moment about the mid- 


is at a distance & = cC,,/(C;,cosa+Cpsin a) 


ipstream of the mid-chord point. From (36), (59), and (60) this distance 


Is found to be 


2C'y sin 38 cos(a+438)—(CZ/)cot $8) (61) 


1 24+ 2C, cosec $8 cos(a+ $8) | 
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For Cg = 0 this gives the classical result ¢ = jc, while for a = 0 we have 
- * 9 C 
= c{siney—521. (62 


We note here that if a source is placed on the aerofoil surface in such a 
way that the centre of pressure is at the mid-chord point, then the lift due 
to the source alone has the coefficient 
(63 
which follows from (9), (59), and (62). 


8.2. The ‘directed’ source at the trailing edge 
Fig. 5a shows the real non-circulating flow about an aerofoil with a 


straight walled duct discharging near the trailing edge; Fig. 56 shows the 


y=0 \ 


/ 7 ~ , 
28\5 ~ A 
C/ 











The directed source. 





Fic. 5b. 


details of the flow from the duct. This flow is also the ‘fictitious’ flow intro- 
duced in section 2—-we have simply put 8 = 8, but bearing the closure 
conditions in mind we have compensated for the loss of this parameter by 
introducing the parameter 45, defined in Fig. 5b. Later in the calculation 
we allow 6, to tend to 6 making J and F coincide, and then we have the case 
of a directed source discharging right at the trailing edge. 

The fact that the z-plane in Fig. 5a is doubly-covered is only aconsequence 
of our using a thin aerofoil instead of one sufficiently thick to accommodate 
the duct. In practice the duct would be curved to lie within the aerofoil, 
but the principal characteristic of the duct as far as the external flow i: 
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meerned will be the direction of flow at its exit, and provided the final 
gction of the duct is straight and several times longer than it is wide it 
3 unlikely that our model will introduce significant errors. 


From Fig. 5 we find that 


7 dar A . y 2a 
os ‘4 Y ) (64) 
\—r (7—0, y < 1+0) 
so it follows from (24), (25), and (26) that 
San a ~~ th, 4-8)—A, (65) 
1-+-cosaA * cos 8)+ (2 — “Joos 5, (66) 
sind = ~(sin8-+sin 8,). (67) 


Equations (66) and (67) show that 6, cannot in general be chosen arbitrarily. 

However, if we restrict our attention to a theory correct only to first order 
in ,/(Cg) it follows from (9) that 6 and hence 65 (89 < 8) are smail first order 
quantities, and equation (66) can be ignored. In this case (65) and (67) give 


@=—-- (3 T So). (68) 


Now (4aP/cU) will differ from unity only by a term of order 56, so to first 
rder (48) and (68) give 
C; 27a-+27(5+5,), (69) 


on putting B 6. In the limit 6, > 6, (9) and (69) then give 


ony 
/ 209 


C,, = 2ma+47 i ) (70) 


N\ 7 
Similarly equation 50) reduces to 
Cy = dm (71) 


whatever value is given to 8,. For « 0 it follows from (70) and (71) that 


he centre of pressure is at the mid-chord point. 


It is interesting to note that if a source is placed on the aerofoil surface 


so that the centre of pressure is also at the mid-chord point, the lift agrees 


to first order with that given in (70), provided we give 7 its ‘average’ value 
of 7. This follows by a comparison of equations (63) and (70). This result 


is what we would expect on physical reasoning, and gives some confirmation 


of the approximations used in the derivation of (70) and (71). 
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9. Concluding remarks 

The main results of some practical interest given in this paper are equa- 
tions (70) and (71). It may seem that a much shorter method of deriving 
these results could be given by simply specifying the flow direction across 
the duct exit D’/ (see Fig. 5b). It could, for example, be taken to be constant 
and equal to 7 over this exit. This procedure would obviate much of the 
algebra given in this paper, but unfortunately it is not a valid approximation 
for source flows. First it is apparent from Fig. 56 that the average value 
of 6 on D’J, say @,, would be somewhat less than 7. The approximate 
theory (see below) actually yields 

Cy, = 40,4 (2CQ), (72) 

and comparing this with (70) we see that 6, = 7/7. Secondly, it is difficult 
to give physical meaning to this artificial boundary condition. It could 
perhaps be regarded as being produced by a uniform downward motion 
of a solid surface placed across the duct exit. 

Equation (72) can be derived as follows. From (69) we see that if 6 = 6, 
over an interval Lis 
44,5, but here 6 will not be related to Cy by (9) as in this theory the source 


m—d8 < y < 7+6, then the corresponding value of C 


has no effect on the mapping of the z-plane on the ¢-plane. If we put 5 =0 











in (40) we find that for a thin aerofoil s 3c(1—cos y), and hence if h is the | 


width of the duct exit it is related to 6 by h/c = 1—cosé = 482, to first 
order. Assuming the velocity in the duct exit'to be U, we have 
Co = AU el hic 


and hence C,; = 40,5 = 40,,/(2CQ). 
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SUMMARY 


This paper presents a new and convenient method of calculating the rate of heat 


transfer from a plane wall into a laminar incompressible boundary layer, for an 
arbitrary distribution of surface temperature and a main stream velocity distribution 
of the form [ cx", where x is the downstream distance from the leading edge of 
the momentum boundary layer. It is assumed that the temperature difference 
between the surface and the main stream is small enough to ignore any buoyancy 
effects. A power law representation is first constructed of the calculated exact 
velocity profile for various values of m, with deviations from the exact profile, up 
to U/U, 0-95, of less than 4 per cent. for 0 m 1 and 7 per cent. for m 4, 
The ensuing temperature equation is transformed into a von Mises form and solved 
by the method of sources. When the surface temperature is independent of x and 
0<n 1, the numerical estimates of heat transfer rates are found to be within 


3 per cent. of previously calculated values, based on the exact velocity profile, and 
the errors are considerably less than those involved in Lighthill’s method, which 
is based on a linear approximation to the exact velocity profile. The method of 


sources leads to a particularly convenient form of solution when the surface tempera- 
ture is given by a series in x, and the numerical results are in agreement with those 
of Chapman and Rubesin who used the exact velocity profile. 


This source method of solving the von Mises form of the temperature equation is 


extended in a following pa 


er to treat the problem of calculating heat and mass 


transfer in turbulent boundary layers. 


1. Introduction 

AN approximate solution of the problem of forced convection in a developing 
laminar boundary layer, for a free stream velocity given by U, = c,2” and 
for an arbitrary temperature distribution on the surface of a semi-infinite 
flat plate, has been obtained recently by Lighthill (1). The Heaviside 
operational method was employed, and the analysis based on a linear 
approximation to the exact velocity profile in the region of flow very near 
the surface of the plate. In the particular case when the surface temperature 
of the plate is kept constant, a formal solution, based on the exact 
velocity profile, is easily obtained and is described by Lighthill (1); the 
associated numerical values of rate of heat transfer are tabulated (1, p. 367) 


[Quart. Journ. Mech. and Applied Math., Vol. IX, Pt. 4 (1956)] 
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for various values of m, together with the numerical values calculated on the 
basis of the linear approximation to the velocity profile. For the Blasius 
case (m = 0) the numerical value given by this approximate solution is 
within 3 per cent. of the exact value, but the error increases as m increases, 
and at m = 4 the error due to the linear approximation is about 22 per cent. 

The analysis given in this paper provides an alternative and convenient 
method of approach. It is based on power iaw approximations to the exact 
velocity profiles for various values of m: these may be constructed so that 
they deviate by not more than 4 per cent. for 0 < m < 1, and 7 per cent, 
for m = 4, from the exact profiles (computed from the tables given by 
Hartree (2)) up to U/U, = 0-95. The temperature equation, based on 
this type of approximation to the velocity profile for a given value of m, 
is first transformed into a von Mises form. The appropriate boundary con- 
ditions are specified by regarding the heated plate, assumed to be infinite in 
the lateral direction, as an assembly of continuous infinite line sources, each 
source being normal to the main stream. A solution of the temperature 
equation is then obtained for a single source, and the distribution of source 
strength (or local rate of heat transfer) is then determined, by solving an 
appropriate integral equation, so that the prescribed surface temperature is 
achieved. In the constant temperature case and when the values of m lie 
between 0 and 4, the numerical values of rate of heat transfer computed 
by this method are within 3 per cent. of the values computed on the basis of 
the exact velocity profile; the errors involved in the application of Lighthill’s 
linear approximation are therefore considerably reduced. 

We find that the source method is particularly useful in dealing with the 
general case when the temperature of the plate is an arbitrary function of 
the distance x from the leading edge of the plate. If, for example, the 
surface temperature of the plate is represented by a series in x, then 
the coefficients of the series describing the temperature distribution in the 
flow over the plate are obtained as simple quadratures, and are in exact 
agreement with the numerical results obtained by Chapman and Rubesin 
(3), when m = 0, by numerical integration of an ordinary differential 
equation for each term in the series. 


2. Solution of the temperature equation in a laminar incom- 
pressible boundary layer over a fiat plate for an arbitrary 
surface temperature distribution and Ll’, = c, 2” 

Let (U’, V) denote the components of velocity, in a laminar incompressible 
boundary layer, relative to coordinate axes Ox, in the plane of the plate 
and parallel to the direction of the main stream velocity, and Oy normal 
to the plate surface. Denoting the temperature by 7’ and the thermal 
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diffusivity by «x, the equation describing the distribution of temperature 
(cf. Howarth (4, p. 760)) is 
oF . we oT 
7 : = . (x ). (1) 


Cx cy cy\ cy) 
We consider the case in which the main stream velocity is given by 


U, = ¢,2", (2) 


z being measured from the leading edge of the momentum boundary layer. 
The distribution of U’ and V needed to integrate equation (1) may be 
calculated from the well-known solution (cf. Goldstein (5, p. 140)) of the 
momentum boundary layer equation. This solution is expressed in terms 


of the similarity variable €, given by 


E = 3[U)/(va)}¥y = 3(c,/v)tamDPy, (3) 
where v denotes the kinematic viscosity of air. The stream function y is 
given by bs (Uy x)! f (é) (ve, iam +DI2F (£), (4) 
and consequently, using the relation U = é/éy, the x component of 
velocity is given by 1) LU, f’ Le, af’; (5) 


the derived function f’(é) has been tabulated for various values of m by 
Hartree (2). We now suppose that numerical values of the coefficient 5, 
and the index 8 can be determined to give a good power law representation 
of the exact velocity profile (5), in the form 


I y C zd, EB (6) 


the accuracy of the approximation being discussed in § 3. The form of the 
function f (and hence of ys) which is associated with the approximation (6) is 
then given by (4) and (6), namely 


> 


f = 2b,(1-+4B)1eH8, (7) 
since sy = f = Oat & = 0; the velocity U at a given 2 value is then expressed 
in terms of the stream function & by the relation 

17 b,¢ : (] B)b{2b,(ve,)*a m+1 2} lich p) (8) 


Using the von Mises method (5, p. 126) we now simplify the left-hand 
side of equation (1), replacing x and y as independent variables by x and wb. 


The basic equation (1) is transformed into 


¢ ‘i Cc «it é fy 
c l : 9 
Cx c I c us | ( 
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when « is indepetu.. » of temperature. Substituting (8) into (9), the 
temperature equation becomes 
oT a ae oT 
i Cy vil? +p)m pj(2 =p) of - 4 3 (10) 
Cx a ou\" ous 
where 


Cy b,c, Ki (c, v)'(1+)(25,¢,) 1}8/4+B) and f B (] 1B), 


Writing 2° = X, where s = (1-++-m)(2+8)(2+28)-!, equation (10) simplifies 
to ny om 
Cc 1 Cc 1¢ 1 

; 1] 

dale al =a) Mm 


We now suppose that a continuous infinite line source, normal to the 


in which c = ¢,/s. 


“-axis, is stationed on the surface = 0 at x = ay and emits Q units of 
heat per unit time, per unit length of source. We write 

n w—x X—Xo. (12) 
With the source stationed at 7 = 0 the basic temperature equation (11) be- 


comes arn OST 
a. (ys! 7} (13) 


on Ob\" Ob 
We seek solutions of (13) in the form 
T —T, = An’g(u), (14) 
where u =n and 7, is a constant, which satisfy the conditions 
associated with a line source, namely 


T +> 2, 7 > O, ub 0. (15) 
T + T,, n > 0, ys > 00, (16) 
and aT => ll n> VU, b> 0. (17) 


The condition that the downstream flux of heat from the source, per unit 
width of any section normal to the plate, is independent of 7 is expressed 
in the form ib 


pe, | (T—T) db = Q, (18) 
where p and c,, denote appropriate values of the density and specific heat 
of air. Substituting (14) into (18) shows that 


A- je 0, (19) 


x 


and A = (Q/pc,) | g(u) du. (20) 


0 
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Substituting (14) and (19) into (13) then leads to the relations 


p 1 /(2—t) = (1+ )/(2+-8), (21) 
,d (dg dg os 
- BoA's) ae? = (22) 
where 
B = (2—t)e = 2b, ¢, «(1+m)-[(e,/v)4(1+B)(2b, ¢,)-2]P!4+8), = (23) 


A suitable exponential form of solution of (22) exists; we find 
g(u) = exp{—u?®" B-1(2—#)-“}, (24) 
leading to the expression for temperature 
T'—T, Ay-Veexp{—u? B-1(2—#)-"}, (25) 
which satisfies conditions (15), (16), and (17). Substituting (24) into (20) 
shows that the quantities A and Q are related by the expression 
A Qf pc,, BUe-(2—t)(-HePf(2—#t)-4J-1. (26) 
Consider now an assembly of continuous sources (of infinite extent in a 
direction perpendicular to the a-axis), the strength Q@ of each source being 
a function of 2, spread along the surface 0, 2 > a,, say, and suppose 
this arrangement leads to a temperature 7\(x) on %¢ = 0, x >a,. The 
temperature excess 7'—7T, at a point P at a downstream distance x > x, 
on & = 0, due to a single source at x = 2, is, by (25), equal to 
A (a8 — ap) “Ne, (27) 
and consequently the temperature excess at P due to the whole assembly 
of sources is given by 
z 


U0 


( (a®—a8)-U2YA (29) dxy = T;(x)—T. (28) 


When 7\(x) is prescribed, this is an integral equation for A(x) and hence, 
by (26), for Q(x), the local rate of heat transfer. Writing 


a Xx, ae d, 1/(2—#) HL, (29) 

equation (28) is transformed into the Abel form (6, p. 229) 
| (X—d)-Hu(d) dé = T,(X)—T, (30) 
where u(d) A(d)h%—-8s-1, (31) 


The solution of (30) is known (6) to be 





T,(X)—T,] dX, (32) 


u(d) (d—X )#-! 
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provided » < 1; the form of equation (21) shows that this condition jg 
satisfied. In § 3 we consider the form of solution in the special case given 
by a, = 0 and 7;,(X) = constant, and in § 4 the case 

y % 1 3 


T,—T, = ¥ a,x". 


3. Forced convection from a semi-infinite flat plate through a 
laminar incompressible boundary layer when the surface is kept 
at a uniform temperature and U, = c, x” 

In this problem a, is zero, J, is the constant free stream temperature, 
and 7; is the constant temperature of the plate. The solution, expressed 
in terms of the u(¢) function, is given by equation (32) to be 


u(d) = (T,—T))7 sin po d4-1, (33) 


and hence, using (31), 


A (7, —7))sa-! sin po xr— V2, (34) 


0 
Equations (26) and (34) then give the local rate of heat transfer as 


Q (2p) [U% (vq) HF, (35) 
where the coefficient 


BP : (7;- Ty) kot? B),, 1p} (24 B)/ ee m)} (2+B)/ l +f) -1/(24 3) , 
< (1+ B/2)2 PEP sin[(1+-B)(2+-B)- aU [(1+B)(2+8) I; 
k is the thermal conductivity and o the Prandtl number. Using equations 
(25) and (34) the associated temperature distribution at a point (a, ~) is 
expressed in the form 


T—T, (7; T),)sz 1 sin pr | (x8 —a§)-e oa 1/2 y, 


0 


x exp{ — Pf? B-1(2—t)-1(a8—a§)-1} day. (36) 


As in § 2 we write x = ¢, a = X, 1/(2—t) = p, and in addition we write 
Pp?) B-1(2—t)-! = y; equation (36) then reduces to 
x 

6 = (T—T%)/(1,—T) = 7 sinpr | (X—¢d)-#d#-1 exp[—y(X —¢)-!] dd. 

0 al 
(04) 
It was noted by Calder (7, p. 170) in another context, that this integral is 
a degenerate case of the confluent hypergeometric function and it may be 
expressed, in terms of Pearson J functions (8), in the form 


0 ] I(y A, ft). (38) 
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Expressed in terms of the original quantities this may be written in the 


form 


6 1 — 1{2v(1+m)b,(1+8)-1(2+B)—'1é?2+P, —(1+8)(2+8)-}. (39) 


We see that Q(x, 
Howarth (4, p. 794). The total rate of heat transfer from a plate of length 


) is proportional to (Up/z)!, i.e. to 2-2, as given by 
is then given by 


| U(x dX, 2(1--m) (ce, y)sal+ml2 Fr, (40) 
0 


Expressed in terms of the Nusselt number Nw | Q(ay) dxy/k(T,—T) 
0 


and Reynolds number R = a, U;(#,)/v = ¢,2}7*"™/v, this finally becomes 


Nuk 1 2U(2+P)(] + B)-N2+B( 2+ B)1+P2+P) sin[(1+B)(2+ 8)-7] > 

P[(1-+B)(2+f)-2]b2@+P(1 +-m)-1+Bx2+B)'GN2+P), (41) 

This expression has been calculated numerically for o = 0-72 and m = 0, 
|, 4 using values of b, and f which were obtained by fitting a power 

law to the appropriate exact velocity profile for each value of m. These 
profiles were computed from the tables given by Hartree (2). It was found 
possible to choose the values of b, and f (given in Table 1) so that the power 
laws deviated by only small amounts (about 4 per cent. at most for m = 0, 
,4, land 7 per cent. for m = 4) from the exact velocity profiles for values 
of € from zero up to the values corresponding to U’/U, = 0-95, i.e. over the 
significant part of the boundary layer. The results based on the power law 
ipproximations and Lighthill’s (1) linear approximation are shown in 
Table 1 together with the values based on the exact velocity profiles. The 
error in the calculated values of heat transfer rates based on the linear 
approximation increases with the value of m, because the deviation of the 
actual profile from the tangent at the origin increases with the value of m. 
Using the power law approximation we see that the error is reduced very 
considerably. The choice of a good power law fit is, of course, to some 


extent arbitrary, but the definite result which Table 1 illustrates is that 


for 0 m t power law representations can be constructed with a 
maximum deviation from the exact profile (up to U/l, = 0-95) of 4 per cent. 
for 0 n 1 and 7 per cent. for m 4, and calculations based on these 


are within 3 per cent. of the heat transfer rates calculated on the basis of 


the exact velocity profiles. 


In the case m 0, the temperature distribution function (39) was 
calculated, for o = 0-6, 1-1, 15-0, using (i) the Lighthill linear approxima- 
tion, and (ii) the power law approximation (b, 0-61, 8 = 0-88). These 


are displayed in Table 2, together with the values calculated by numerical 
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TABLE 1 


Calculated rate of heat transfer from a flat plate at constant temperature 








m , . ; 0 3 i 1 4 
b, ; , : , : , 7; . | 0-61 0-81 0:99 1-35 1-78 
Bp 0-88 0-86 0-80 0-76 0-66 


Per cent. maximum deviation of power law pro- 
file from exact profile (up to U/U, 0-95) . 3 3 4 3 7 
NuR- based on exact profile P m = 0-585 0-596 0-576 0-495 0-325 
NuR-* based on linear approximation to profile 0-601 0-648 0-653 0-587 0-398 

NuR-} based on power law approximation to 


profile ‘ ‘ ; , : ‘ 0-569 | 0-610 | 0-587 | 0-512 | 0-331 
Per cent. error due to linear approximation 4 +3 +9 +13 +183 | 99} 
Per cent. error due to power law approximation —3 +2 +2 +3 2 





integration of the formula, given by Howarth (4, p. 781), which is based 
on the exact Blasius velocity profile. The differences in the temperature 
distributions are seen to be very small. 


TABLE 2 
Calculated values of (T —T))/(T,—T)) = 0, the temperature distribution func- 
tion, when m = 0, for o = 0-6, 1-1, and 15-0: 0,, 43, Ac refer to calculated 
values based on the exact velocity profile, the power law approximation 
(b, 0-61, B 0-88), and the linear approximation (b, 0-664, B 1-0) 


respectively 














G 6 Il 15 

é 04 6; 6, 04 Op A 64 6; 6, 
0°25 0O°139 | O'I4C O°143 | O'I7I 0°173 | O'175 | O°406 | 0°423 | O42 
o"50 0°275 | 0°279 | 0°25 0°340 | 0°343 | 0°345 | O°752 | 0°706 | O°755 
o°75 o°'410 O'414 O°422 0502 505 O°512 0°946 0948 0°945 
1°00 o'54I 0°545 554 | o650 | 0650 | 0659 | O'995 | O°995 | 0°995 
1°25 0654 0055 o°O7I 0769 o'77!1 o°750 

1°50 | O'751 | O'757 | 0-773 | 0862 | 0-864 | 0°873 

1°75 0°535 545 0°560 0°905 0°960 | O°972 

2°0 o'dsgI 5 O°gI2 





4. Forced convection from a semi-infinite flat plate through a 
laminar incompressible boundary layer when Ll’, is independent 


of x and 7,—T, is given by 5 a, x” 
0 


The method of analysis given in this paper is found to be particularly 
convenient when the surface of the flat plate is kept at a temperature 
which is a function of distance x from the leading edge. The only available 
calculations based on the exact velocity profile are those by Chapman and 
<ubesin (3) for the case when 7; — 7, is given as a series ¥ a, x” and m = 0. 
n=0 
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It is interesting to compare our results with those based on the exact 
velocity profile in this case. 
Writing os 
om n 9 
T; : (i > a,2 (42) 


in equation (32) we obtain 


na 1 g ”(n/s) 
u(d) N\ @,, sin par d AX — 


Lund T dd J (d— X)! “4 

n=0 ; 

7 sin ur ) a,(n/s + pw)p ste-DB (n,n s+1), (43) 
n=0 


where B,, is a beta function (6, p. 253) and, following the analysis given 
in preceding sections, the local rate of heat transfer at a point 2 is given by 
ox 


V(x) so (var) |? > a, x2 1+BX2+8) ‘w(1 | B)(1+Bx2+B) '(2+-B)-ue+B) 4 


n=0 
sin{a(1+)(2+8)—P{(1+B)(2+f)-Hoye+P) oil2+B( 1+ 2n) x 
B,{(1+B)(2+B)-, 14+2n(1+8)(2+ 8) 3], (44) 


where U), is the uniform free stream velocity, and the other quantities have 
been previously defined. Using the linear approximation to the velocity 
profile (8 = 1, 6, = 0-664) this reduces, as expected, to the formula given 
by Lighthill (1, p. 369) using his operational method; for o = 0-72 this 
has been shown by Lighthill to yield results which are within 3 per cent. 
of those obtained by Chapman and Rubesin. A calculation of (44) based 
on the power law approximation (8 = 0-88, b, = 0-61) reduces this error 
to less than 1 per cent. 

The temperature distribution in the flow over the surface, following the 
analysis described in § 3, is then given by 


'—T (277)! sin[(1+-8)(2+8)-17] > a,(1+2n) > 


n 


B,{(1+f)(2+8)-, 1+ 2n(1+f)(2+8)-} | (a’—a§)-Veay 


exp{—?-(2—#t)-1B-1(as—a8)-l\ day. (45) 


Writing 2x dh, x X, and 4/X u, this reduces to the form, discussed 


by ( hapm in and Rubesin, 


(46) 
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where 
4 l . - “7 2 i) i) g 
X,,(€) = 771(1+8)(2+ 8) sin[(1+)(2+8)—7](1+ 2n) > 
¢ Q\/< : 1 2n) 
3, {(1+B)(2+B)-4, 14+-2n(1+B)(2+B)-} | (1—w)-Humt+2m— 
X< eXp} 26, vU 5 ly B) 1(2+-B) Ip~-1(] w) 1g2 P} du. (47) 
The first term (n = 0) in the series reduces, as in § 3, to a Pearson J fune- 
tion, and we obtain 
r ¢ Q é 28 ¢ 
X 9(€) 1— 1{2b, ve-1(1+-B)-1(2+-B)-1€2+P, —(14+-8)(2+8)-!}. (48) 
The X,,(€) functions were obtained by Chapman and Rubesin by numerical 
integration of an ordinary differential equation in the cases n = 0, 1, 2, 3, 
4, 5, 10. A considerable simplification of the calculation is thus obtained 
by using a power law approximation to the exact velocity profile. For 
n 0 we find from Pearson’s J Tables (8) that the expression (48) gives 
numerical results in exact agreement with those given by Chapman and 
Rubesin (3, Fig. 5). It has not been found possible to reduce the form of 
the integral in equation (47) for n > 0, but it has been evaluated numeri- 
cally for a number of values of € for n = 1, 2, and 10, and the ensuing 
numerical values of X,,(€) are found to be in agreement with those given by 
Chapman and Rubesin. 


5. Conclusions 

The method used here of calculating the rate of heat transfer from a flat 
plate into a laminar incompressible boundary layer, for an arbitrary distri- 
bution of surface temperature and a main stream velocity distribution of 
the form U, = ¢,x™, is based on a power law approximation to the exact 
velocity wcilte for various values of m and leads to convenient formulae 
describing the rate of heat transfer and the temperature distribution. 
These power laws may be constructed so that they deviate by not more 
than 4 per cent. for 0 < m < 1, and 7 per cent. for m = 4, from the 
exact velocity profiles up to U/U, = 0-95. We find that in the case of uni- 
form surface temperature, the errors involved in Lighthill’s linear approxi- 
mation, giving rate of heat transfer for m = }, 4, 1, and 4, are considerably 
reduced by using the appropriate power law. 


«o 
The case in which m = 0 and 7,—T, = a,x” has been discussed by 
0 
Chapman and Rubesin (3), and the tedious method adopted by these 
authors of computing the X,,(é) functions in the solution 


T- T) ya ay Ly X,,(€) 


n=0 
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; considerably simplified by our method of approach. In particular the 
{,(¢) function is expressed in terms of a tabulated Pearson J function (8) 
nd the X,(€) functions, for x 0, are expressed as simple quadratures. 
Their numerical evaluation replaces the Chapman and Rubesin step-by- 
step numerical integration of an ordinary differential equation for each 
Jue of n. The numerical results are found to be in agreement with those 


f Chapman and Rubesin 
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ON THE CALCULATION OF HEAT AND MASS 
TRANSFER IN) LAMINAR AND TURBULENT 
BOUNDARY LAYERS 


Il. THE TURBULENT CASE 
By D. R. DAVIES and D. E. BOURNE (University of Sheffield) 
[Received 1 December 1955; revise received 17 May 1956] 


SUMMARY 


The method of analysis described in Part I (7) is applied in this paper to obtain 
solution of the problem of heat transfer from a flat plate through a fully turbulent 
boundary layer, when the surface of the plate is kept at a steady constant temper 
ture ; the temperature difference between the surface and the main stream is assumed 
to be small enough to ignore any buoyancy effects, and the main stream velociti 
is independent of downstream position. The equality of heat and momentum edd 
diffusivities is assumed, and the analysis is based on power law representations ¢! 
(i) the experimentally determined velocity profile, and (ii) the distribution of eddy 
viscosity, deduced from Townsend’s experimental work concerning the structure o/ 
the turbulent boundary layer on a plane wall. The ensuing theoretical temperatur 
distribution over the plate and the distribution of local heat flux along the surfa 





of the plate are found to be in close agreement with the results of experiments | 
Elias. 

On the basis of equality of momentum and vapour «ddy diffusivities the meth } 
is extended to treat the problem of calculating the rate of evaporation from 
saturated liquid surface into the developing turbulent boundary layer adjacent 
a smooth wind-tunnel wall. The analysis takes into account the effects of tl 
transition and laminar sub-layers over the evaporating surface, and the theoreti 
values of rate of evaporation are found to be in reasonable agreement with the result 
of evaporation experiments in a wind-tunnel. 


1. Introduction 

MEASUREMENTS have been made, e.g. by Pasquill (1) and Davies at 
Walters (2), of the rate of evaporation from saturated liquid surfaces int 
developing fully turbulent boundary layers, and theories of turbule! 
diffusion in boundary layers, see, for example, Sutton (3), have been disp 
cussed in relation to experiments of this type. However, these theor 
are strictly applicable when (i) the thickness of the momentum boundat 
layer is independent of downstream distance so that, if the downstreat| 





pressure gradient may be ignored, the Reynolds shearing stress is indepe} 
dent of position in the boundary layer, and (ii) the contribution to t 


(Quart. Journ. Mech. and Applied Math., Vol. IX, Pt. 4 (1956)] 
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total downstream flux of vapour within the laminar sub-layer and transi- 
tion layer near the surface may be ignored. In large-scale evaporation 
experiments in the open atmosphere the downstream flux in these sub- 
yyers may be neglected, and it is also possible that condition (i) is a good 
pproximation in the atmospheric case. However, in the developing 
turbulent boundary layer on a plane wind-tunnel wall this condition is 
not valid, as Townsend’s recent observations (4) show, and we shall 
lemonstrate in this paper that the downstream flux in the laminar sub- 
yer and transition layer can be quite considerable. A theoretical treat- 
ment is thus required which does not depend on the validity of conditions 

and (ii) 

The method of obtaining a solution of the diffusion equation in terms of 
simple similarity variable, applied recently by Davies (5) to the problem 
fheat transfer from a flat plate at constant temperature, is not suitable 
when the leading edge of the vapour boundary layer is at some distance 
lownstream from the origin of the momentum boundary layer—as it is 
in general in wind-tunnel evaporation experiments. We introduce an 
ternative method in which we regard the evaporating area, assumed to 
be infinite in the lateral direction, as an assembly of continuous infinite 
ine sources, each being normal to the main stream velocity. The distri- 
bution of source strength (or local rate of evaporation) is determined so 
that the surface vapour concentration is zero upstream of the leading edge 
nd a constant (the saturated value) downstream of this edge. The analysis 
is based on a power law representation of the eddy viscosity distribution, 
given by Townsend (4), and is applied firstly to the problem of turbulent 
flow over a flat plate kept at a constant temperature. It is found to lead to 
theoretical results for rate of heat transfer which are within about 3 per cent. 
of the results of experiments by Elias (6) and of the theoretical results 
obtained by Davies (5), using the similarity approach. The method of 
sources, described in Part I (7) and incorporating in addition an approximate 
treatment of the flux in the transition and laminar sub-layer, is then applied 
to the evaporation problem and we find that the theoretical results for rate 
of evaporation are in reasonable agreement with the results of experiments 
by Davies and Walters (2). 


2. The temperature equation in a fully turbulent boundary layer 
over a Semi-infinite flat plate 
Let (U,V) denote components of mean velocity, in a fully turbulent 
boundary layer, relative to coordinate axes Ox, in the plane of the plate 


ind parallel to the direction of the free stream, and Oy, normal to the plate 
suriace. Denoting the mean temperature by 7’, the thermal diffusivity by 
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bution of temperature, given by Howarth (8, p. 821), is 


oe eee e[ eT oT" 

l ; + ; ; K - én - ° (] 
Cx cy cy\ cy cy 

constructed by generalizing Townsend's results (4) and writing the momen- 

tum eddy viscosity in the form 


e = KU, xh(é), (2 
where x is measured from the origin of the momentum boundary layer and 
the similarity variable € is given by 

ra y (kat): (3 
the parameters k, A, g, and p are derived from the observed velocity profiles 
in the particular conditions to be considered, and h(é) is obtained from 
Townsend’s results as in previous work (5): Uy is the main stream velocity, 
We now suppose that Townsend's values of h(€) may be adequately repre- 
sented by a power law form 

h(é) ag, (4 
a detailed discussion of the problem of evaluating a and «and the dependence 
of the theoretical transfer rates on the numerical values of these parameters 
being given in sections 3 and 4 of this paper: in this connexion we note 
also that it has not been found possible to carry out the analysis on the basis 
of a polynomial representation of the function h(é). A mathematical 
treatment of the complex differential equation obtained by substituting 
equations (2) and (4) into (1) has not yet been achieved but, apart from the 
extremely thin region in the neighbourhood of the laminar sub-layer and 
transition layer, €,; is numerically considerably larger than « which we 
therefore ignore in the fully turbulent region of flow. An approximate 
correction, of the type used previously by Davies (5), to allow for the effect 
of the sub-layer is later incorporated in the analysis described in sections 
3 and 4. In this paper we write our basic diffusion equation in the form 

yer per é («7 (5) 
Cx cy cy\ oy 
using the result obtained by Davies (5) that «,, = ¢ in forced convection. 
Following the Reynolds turbulence technique of taking mean values over 
a sufficiently long time-interval at a point, the continuity condition for an 
incompressible fluid applies to the mean flow, i.e. 
oU. av 


Cx cy 


= 0 


> 


x, and the eddy heat diffusivity by ¢,,;, the equation describing the distri. 


It has been shown by Davies (5) that a consistent diffusion model may be 
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and hence we may write 
ow . ow : 
| and U: ; (6) 
Cx cy 
where ys is the stream function associated with the mean flow. Using the 
yon Mises method (9, p. 126), as in Part I, we now simplify the left-hand 


side of equation (5), replacing x and y as independent variables by x and w. 


The basic equation (5) is transformed into the equation 
oT — cof ,or 
K Uae él ). (7) 
Ca Cw Cds 


We find that Townsend’s mean velocity profile may be represented by a 


power law, U/U, |-72€°, quite accurately (to about 1 per cent. of the 
ictual velocity profile) over almost the whole thickness of the boundary 
layer, and very accurately in the crucial inner zone € < 0-01. We write 
| cre a | | , a ~ 2 

n genera U/U, bEP, (38) 


Following Davies (5) we then write 


wb = xf (€), (9) 
ind, using equations (6) and (8), we find that 
kU, b(1 + B)2€1+8, (10) 


Substituting expressions for € and U in terms of ys, obtained from (8), (9), 


ind (10), into equation (7), and introducing a new variable defined by 
a > (11) 
where s »—q(a+B 8) eads to the equation 
l+p—q |+-8)-}, leads to tl juat 
( {ig C oT 
. c (u' : (12) 
cx Cus ous 
where f (a 3)(1 3)-1 and 


c KabU2(1 + B)%+P0+8)"*(kbU,) +P 3) *g-1, 

The temperature equation (12) is identical in form with the basic equation 
of Part I (7, equation (11)), and the analysis of the source method described 
in Part I (from equations (11) to (32) of that paper) is therefore applicable 
in the physical context of this section. For a prescribed surface tempera- 
ture, the distribution of source strength (or local rate of heat transfer) can 
then be computed from equation (32) of Part I. The particular case of 
constant surface temperature on a flat plate is considered in section 4, and 
in section 5 the analysis is extended to treat the problem of evaporation, 
when the leading edge of the evaporating area is positioned downstream 
from the leading edge of the momentum boundary layer. A correction is 


made in both cases for the presence of laminar and transition sub-layers. 














472 D. R. DAVIES AND D. E. BOURNE 


3. Forced convection from a semi-infinite flat plate through a 
turbulent boundary layer when the plate is kept at a uniform 
temperature 
In order to test the validity of our analysis based on a power law repre- 

sentation of the eddy viscosity distribution, we first apply the method of 
sources to the problem of heat transfer from a semi-infinite flat plate 
through a turbulent boundary layer, when the plate is kept at a constant 
temperature and there is no appreciable completely laminar layer near the 
leading edge. This case has been previously analysed by Davies (5) using 
the actual values of viscosity given by Townsend (4) over the whole thick- 
ness of the boundary layer. Comparison of the results obtained by the two 
methods serves as a check on the method of analysis described in this paper 
and provides support for its application to other cases. 

We first solve the problem with the sub-layers neglected. The tempera- 
ture of the plate is taken to be independent of x and is denoted by 7° in 
this case; the main-stream temperature is denoted by 7). With a, = 0 and 
using the appropriate values for s, t, and c, evaluated in the preceding 
section, the solution of the integral equation (32) of Part I (7) may be 
written in the form 

A (7; —T),)s2-(sin pr)xy tt 'O + p)(1+B)—a(a +B)K2+B—a)~* (13) 


where pe 1/(2—t) = (14+-B)/(2+B—a). (14) 

The variable A is related to Q by equation (26) of Part I to give the local 

rate of heat transfer at a downstream distance x,, and the parameter B 
involved in this equation takes the new form 
B = KabU}(2+-B—«a)(1+f)@+Pa+B)* x 

x (kbU)) “+ PAP (1 + p)(1+-B)—g(a+B)]*. (15) 


Following the analysis discussed by Davies (5, p. 330) it may be shown 


that the representation of h(é) by aé* leaves the result 1+-p = 2q, obtained 
previously (5), unchanged in form and consequently (13) simplifies to 
A = (T—1%)sa- x9-! sin pz. (16) 


The derivation of the temperature distribution function then follows as in 
section 3 of Part I, and this distribution may be written in the form 
(T—T)/(7,—T) = 1-1 (y/X, —p), (17) 
where J is the Pearson J function (10) and y = y2~“B-1(2—1#)-!: using 
equations (9), (10), (11), (14), and (15) the ratio y/X in the turbulent case 
may be written in the form 
y/X = k*gbé2+8-% Ka(1+p)(2+B—a)] = u, say. 


We now introduce an approximate treatment to allow for the presence 
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of a laminar sub-layer and transition layer following the method, intro- 
duced by Davies (5), in which use was made of empirical mean velocity 
forms. These empirical formulae (cf. Howarth (8, p. 824), Goldstein (9, 
p, 334), and Preston (11)) are based on measurements over a wide range 
f Reynolds numbers and are expressed in the forms 

U/U, U.yjv=y, (y, <5) (18) 


in the laminar sub-layer, 


U/U. 3°05+5log,y, (5 < y, < 36) (19) 
in the transition layer, and 
T/T 5-5+2-5log.y, (y, > 30) (20) 


in the fully turbulent layer. lis the ‘friction velocity’, which we determine 
isa function of x by fitting the empirical form (20) to the measured velocity 
profile (’(y) at various distances from the leading edge of the plate, and v 
is the kinematic viscosity of air. Equations (18), (19), and (20) are known 
to give a good description of the mean velocity profile very near a wall. 
As indicated by Preston (11), there is some variation in the values of the 
coefficients given by various experimenters, but we find that this scatter 
ultimately leads to no significant variation in the predicted rates of transfer 


or temperature distributions. 


TABLE | 
Calculated heights of laminar sub-layer, transition layer, and fully turbulent 








laye r over a flat plate when U, 3,500 cm./sec. 
x (cm.): 10 20 30 40 50 

Laminar sub-layer height (cm.) calculated from 

equation (18 : : 0-004 | 0-005 | 0-005 0-005 0-005 
Transition layer height (cm aleulated from 

equation (20 7 . 0-02 0-03 | 0-038 0-03 0-03 
\pproximate transition layer height (em.) cal- | 

culated from equation 2] > 2 0-01 0-02 0-03 } 0-03 0-03 
Approximate height (cm.) of fully turbulent 

layer taken from Elia experiments . ‘ 0-3 0-5 0-7 0:8 1-0 





To calculate the distribution of temperature, and hence of downstream 
heat flux, at a given downstream distance 2,, say, we first calculate, from 
the known values of U(x), the height y, = 30v/U, of the transition zone 
at this particular x value. The corresponding value of €, say €, = y,/(k24), 
is then caleulated and, following the method used by Davies (5), we take 


the values of Yr (x) satisfying 


Pe Yr (kar) (21) 


+ 


to be our approximation to the top of the transition zone for x < 2,. As 
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seen later in this section, this procedure yields a value of 7’; which is inde- 
pendent of a—this is necessary to ensure the validity of the analysis leading 
to equations (16) and (17). Taking 2, 50 em. and Ul, = 3,500 em./sec., 
Table 1 shows that these values of y, approximate very closely to the 
values of y,(x) given by equation (20), for values of x between about 20 em, 
and 50 em.: for values of x less than about 20 cm. there is a small error. 
but the yp values are certainly well above the top of the laminar sub-layer, 
The error involved in taking y,(x) to be the height of the transition layer 
is thus localized to values of x less than about 20 cm., and represents only 
a small percentage of the total thickness of the turbulent boundary layer; 
this is not likely to influence the theoretical values of temperature and flux 
at downstream distances of the order of 2 = 50 cm. In fact it has been 
shown by Davies (5) that the theoretical values of temperature based on 
this method, when applied to a plate of length x, = 50 cm., for example, 
are in very close agreement with the experimental results of Elias at 
x 30, 40, 50 em. (5, Fig. 2). It was also shown (5) that the calculated 
values of heat flux from a section of the plate are in agreement with the 
experimental values given by Elias. 

The flow is now taken to be fully turbulent in the region € > &, and, 
if 7* denotes the temperature at the assumed lower limit of the fully 
turbulent layer, then 


(T$—T)|(T,—Ty) = 1-H (ty, — 1), (22) 


2) 
where u, = k®qbé2+?-~[Ka(1+f)(2+8—a)]-l. Using (17), we obtain 
finally ‘yy /J* 4 bl yi* 

. 6, = (T—T)/(%—T7) 
eal , Tk mM 92 
1—(T—T,)/(T7—N) = (1 (u)—L(u,,)/—1(u,)], (23) 
which is the result obtained by Davies (5, p. 330) by a different method. 

We now assume that the mean shear stress and vertical heat flux in the 
laminar sub-layer and transition zone are independent of y. Denoting the 
mean shearing stress by 7, and the vertical heat flux by q) we have, following 
Howarth (8, p. 823), rom 

(J ) Yo ply l T le (24) 
where 7" is the ‘friction temperature’. At x = 2,, on the surface € = &,, 
an expression for gj may also be derived from equations (17) and (22). 
We obtain 


ee wien 
Vo —pC, KU, ax} al: *), 4 
pe, KU, k-2?-“(T%—T7)N, (25) 


where 
a|k?qba 1k -1(] +B) ol x2+B—a) (2+ 8 y)(1+812+B—-a)-? 


N —= ers — a 
x)“ }[1—L(w,) Je" 


; (26) 





P{(1—a)(2+B. 
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An expression is given by Howarth (8, p. 831) for the temperature in the 
transition zone. This is applied here to evaluate 7 at the top of the 
transition zone, and is written in the form 


(T7,—T7f)/T. = 5log,(50+-1)+-5e, (27) 


where 7; is the actual (constant) temperature of the plate and co is the 


l 


Prandtl number v/«. By eliminating g, and 7. from equations (24), (25), and 
" 


{= é 


we then obtain the result 

T,—TY)/(%—T7) 5U, KNkOAU Sa? Ulog,(50+1)+0}, (28) 
giving 7% 

We now find by numerical evaluation of the friction velocity that, in 


the particular conditions of Elias’s experiments (Davies (5) gives p = 0-78 


ind q 0-89), U_ is proportional to 2°", i.e. to x?-4, and 7, (or pU?) is 
proportional to x~°*? which is in approximate agreement with the result 


given by Goldstein (9, p. 362). We note that the factor .V, given by (26), is 
independent of x, if we take the values of y, given by (21) as our approxi- 
mation to the height of the transition layer, and consequently 

ry ry riy ry 

(7, —T)/(1j— 77) 


is independent of x. The expression 


M—T  (%—T)[,_ (h—-T]7 
1—T, (TIT (M-TH. 
1—J(u,)|1+5U, KNk1U 71 x? -Ulog,(50+1)+o}]-1, (29) 


for 7’,, obtained from (22) and (28), is also independent of x: as mentioned 
earlier in this section this is a necessary condition for the validity of the 
analysis. 
The temperature distribution function is now given by 
(T—TT) (T,—T%) 
T—T, (%—TT) (%—T7) 
1, —T, (7,—T77T) 
(1)—T7) 
I(u)—I(u,) , 5KU,2 
Las * kU_ag-? 
5AU, Nilog,(50+1)+c} 


TT) ~A-—p 
kU, x4 





Vilog.(50+1)+-c! 


: (30) 





] 


{n expression for the total flux of heat, Q,(2,), across the top of the 
transition zone into the fully turbulent zone, obtained by integrating 


25) with respect to x. 1S 


Qala) = Kpe, Ul TI—Ty)N (hq) 2h (31) 


1 
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we note that expressions (30) and (31) are identical with the corresponding 
expressions obtained by Davies (5) if N-! is replaced by the integral 


x 


F dé calculated by Davies. 





é 
An expression for the downstream flux of heat, at x x,, across the 
transition zone, per unit width of a plane section normal to the flat plate, 
can be obtained by using the empirical velocity and temperature forms 
(8, p. 831). This contribution to the total flux is given by the integral 
300/U; 


pe, | U(T—T) dy 


30r/L 
pe,,(T; T)U, | 3-05-4 5log,(U, yv 1)] ’ 
5v/U;r 
a log. [1—o+oU,y (5»)] bo dy. (32) 


| log (50-4 1)+o+kU,24-! (ONKU, yy 
Similarly, it can be shown that the downstream flux hes the laminar 
sub-layer is given by the integral 
5v/U; 
pC» | U(T—T,) dy 
‘ aril l25ckl / ; ; 
Ae tn is — 3kU, af PNIALI5KI “lop, (50+1)+o}]" " 
The sum of (31), (32), and (33) then gives the total downstream flux of 
heat across x = 2}. 

An alternative and very convenient method of calculating the flux may 
be obtained by first integrating equation (16) with respect to x), noting 
that (16) neglects the presence of the sub-layers. We obtain, using also 
(29), the form 
Q'(x,) 

pc, (T,—T)[1—L(u,) a 1+5NKU, x?-%(kU,)—-{log,(50-+-1)+o}]-?> 


1 U7 ] p(l—a\1-4 | Kl ane : 
a (= ) a -B 3 xa \ =|" | es Bx1+8) | (x)x4 sin(7p), 
7 T ( 

P d (34) 
where pu (1+-8)(2+8—.«)-!. On the same basis the flux across the region 
from y, = 0 to y, = 30 is given by 

30v/U; 
PCy U(T - T)) dy 
0 si P 4s 30v/ U7 
(7; Ty)PCp bl n( (ka, )-P i 


‘ — = > dy. 
I(u,) [1 5N KU, 2?-“(kU,)“log,(50-+1)+0)] Ay 


(1 
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The result obtained by subtracting (35) from (34) then gives the flux into 
the fully turbulent zone, and the total downstream flux across x = 2, is 
obtained by adding to this the flux in the sub-layer given by the sum of 


(32) and (33). 


Numerical evaluation of these expressions for the flux and 6 distributions 


depends on a determination of appropriate values of the parameters a and a. 


Fig. 1 illustrates that a power law may be fitted fairly well to the « distri- 
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Fic. 1. Distribution of eddy viscosity given by 
Townsend. h = e/(#U,), where e = eddy viscosity, 
’ distance from virtual origin of momentum 
boundary layer, U, = free stream mean velocity; 
é y/x, where y height above the surface; 
power law h 0-00065€9,  ------ power law 


h 0-00088E°48, 


bution given by Townsend for € < 0-01, but for € > 0-01 there is clearly 
a buffer region between the intense turbulence in the region immediately 
above the transition zone and the smaller degree of turbulence in the main 
stream, and it is not possible to represent by a power law the actual e 
values over the whole thickness of the boundary layer. However, the 
arguments given by Fage and Falkner (12) and Lighthill (13) in the laminar 
case suggests that the crucial part of the boundary layer in the heat transfer 
problem is the inner zone, € < 0-01. We find that the calculated distribu- 
tion of 6, based on a power law representation of ¢ in this inner part of the 
boundary layer, is in agreement (apart from a small deviation of, at most, 
2 per cent. for values of € greater than 0-025) with the values of @ calculated 
by Davies (5, Fig. 2), on the basis of the Townsend e distribution taken 
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over the entire thickness of the boundary layer. This result shows that the 
outer layer, € > 0-01, is only important in controlling the rate at which 6 
approaches its asymptotic value 1-0 in the free stream, and the error 
incurred in using a representation of the inner boundary layer over the whole 
thickness of the layer is small. 

In order to obtain a reasonably good power law representation for 
€ < 0-01 we find that the numerical value of « must be chosen to lie between 
0-40 and 0-46 and the associated value of a must lie between 0-00065 and 
0-00088 respectively. However, we find that the differences in calculated 
values of rates of heat transfer discussed in this section, and of evaporation 
rates discussed in section 4, are not significant when based on 


(i) h = 0-00065€% and (ii) h = 0-00088E%°, 


These are the extreme power laws (shown in Fig. 1) which lead to a reason- 
able fit for the inner boundary layer € < 0-01. Using both power laws (i) 
and (ii) we now calculate the factor V, having determined appropriate 
values of g, b, and £8 from Elias’s measurements (6), following the method 
discussed by Davies (5). The numerical value of the parameter k is chosen 
so that Elias’s velocity profile fits the Townsend profile U/U, = 1-72€°, 
where in Elias’s flat plate experiments € = y/(kx®). Since the velocity 
profile forms then agree and the flows are over a plane smooth surface with 
no buoyancy effects, we then assume, as in (5), that the values of h(€) given 
by Townsend apply in the conditions of Elias’s experiments. We find that 
the numerical values of N-!, based on power laws (i) and (ii) are within 
: 


4 and 3 per cent. respectively of the equivalent factor 
p 


F dé which was 





evaluated by Davies (5) on the basis of the actual €« values over the whole 
thickness of the boundary layer. Consequently, taking o = 0-72 for air, 
the numerical values of the flux expression (31) for Q,(2,), using (i) and (ii) 
respectively, are within 4 and 3 per cent. of the value calculated previously 
by Davies (5). The numerical values of Q(x,)/(7,—7)), where Q is the total 
flux, calculated from ¢ uations (32) to (35), using power law forms (i) and 
(ii), are 0-58 watts per U. and 0-575 watts per °C. respectively for a length 
of 50 cm. of plate and a free stream velocity of 35 m./sec. These results 
are both within about 3} per cent. of the experimental value of 0-56 watts 
per °C. given by Elias (6). 

This calculation can also be used to illustrate the importance of the 
presence of the laminar sub-layer and transition zones in calculating the 6 dis- 
tribution, although it has been calculated by Davies (5) that, at the higher 
speeds involved in Elias’s experiments, only 5 per cent. of the total down- 
stream heat flux takes place in the sub-layers. We find, for example, that 
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the value of @ given by equation (30) at the top of the transition layer is as 
high as 0-56; we note that this value is also indicated approximately by 
Elias, although there is considerable scatter of his experimental points in 
this region. Moreover, we find that, if a solution is computed which leaves 
out the laminar sub-layer and transition zone in equation (30), the caleu- 
lated 6 values are about 25, 7, and 2 per cent. less than the experimental 
) values at é 0-01, 0-02, and 0-03 respectively, thus emphasizing the 
considerable error involved in this case if we leave out the sub-layers. 
The analysis described in this section therefore indicates that the method 
of approach based on a power law representation of the eddy viscosity in 
the inner part (é 0-01) of the fully turbulent boundary layer, together 
with the approximate treatment given of the laminar sub-layer and transi- 
tion zone, gives numerical results for 6 and Q in close agreement with 
those derived by a previous method, which was based on the actual values 
of e over the whole boundary layer thickness. These results are also in 
close agreement with the experimental values given by Elias. This suggests 
that we may extend the treatment to the problem of evaporation into a 
developing boundary layer from an area source whose leading edge is at a 
distance downstream from the leading edge of the momentum boundary 
layer; a solution in terms of the similarity variable € is not attainable in 


this case. 


4. The problem of evaporation into a developing fully turbulent 

boundary layer 

The rates of evaporation calculated in this section are based on the 
conditions recorded in the experiments of Davies and Walters (2) at 
Sheffield. In these experiments the rates of evaporation from various 
rectangular saturated surfaces (of aniline) were measured in order to deter- 
mine the dependence of evaporation on the dimensions of the area at a 
given wind speed. In conformity with previous notation we replace pc,, T 
by x, the mean concentration of vapour, in the analysis of section 2. If we 
measure x from the origin of the momentum boundary layer and suppose 


the area source is positioned with its leading edge along x = a,, the con- 


ditions of the problem can be expressed in the form 


" Xo (the up-stream value), z< mh, y > 9, (36) 
y, (the saturation value), z>4,, y = @, (37) 
ind xX > Xo 2 ay; and Yy —> 0. (38) 


We now assume that vapour and momentum diffusivities are equal; the 
ippropriate solution of the integrai equation (32) of Part I (7) then 
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becomes 





u(h) = (xi—xXo)77! (sin pa 


F ¢ 
6 y , 
_— | GX ax 
dd . 
aj 
’ = . s a ‘ 
= (xi—Xo0)7 1 sin ur(d—aj)e-}, (39) 
where x; is taken to be the constant value of y on the area when the sub- 
layers are at first neglected in the analysis. The expression for A becomes 
" rae — 1 ae Pd 8 8\(u—1)>(s—1 
A = sz—sin ua(x1— Xo)(73§—af)4-Pa—Y, (40) 
and consequently the equation for @, given by (37) in Part I, is replaced by 
1 . | 2 : I 


(x—x0)/(x1—Xo) 
(X— A) 
=atsinuz | (X—A,—¢)-¥d4—-1exp{—y(X—A,—¢)—} dd, (41) 
0 
where A, = aj. As in Part I the integral can be expressed in the form of 
a Pearson J function and we obtain 
(x—Xo0)/(xi— Xo) = 1—HyX “01 —(q/2")7, — 4}, (42) 

where p = (1+ 8)(2+8—«)-}. 

Now in the wind-tunnel evaporation work of Davies and Walters (2) the 
numerical values of » and q were not calculated; it was only possible to 
measure the mean velocity profile in the neighbourhood of the evaporation 
plate and no measurements could be made of the mean shearing stress 
distribution. However, the experimental investigations carried out by 
Townsend and by Davies and Walters were both over a smooth wind-tunnel 
wall, and we assume that the result p = q = 1, shown to be suitable in 
Townsend’s work, also applies to that of Davies and Walters. The variable 
yX~—1{1—(a,/x)*]-* in expression (42) may then be written in the form 

k*b[y(ka) 1}(24 B—a) ae 
Ka(1+-B)(2+f—a)[1—(a,/z)*9]} 7° 

The approximate treatment of the transition and laminar sub-layers is 
now introduced. If we denote by x, the mean concentration at the top of 
the transition layer, calculated from (20), at 2 = x,, the downstream dis- 
tance at which we wish to evaluate the y distribution and flux, then 


(X%—Xo)/(x1—Xo) 1—I(v,, —p), (43) 
where v,, corresponds to the top of the transition zone, at x = 2,. The 
function corresponding to 6, may then be obtained following the method 
of section 3. 

We assume that heat flux and mean shearing stress in the sub-layers are 
independent of y. Following the method described in section 3, we equate 
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the expressions for vapour flux at the upper and lower limits of the sub- 


layers at x = 2,, and obtain 


(Xa XI) AXo Mas 
SKU (kU,)—[log.(50,-+ 1)+-0, {xj (aj —ajq) 43" aX2+B—a)"* MP, (44) 


where 
al k2ba-1.K -1(1 + B)-1]4-o2+B-2)1(9 4 B— (1+ BX2+B-a)~* m 
M | ‘i Mn so. tesco 
M{(1—a)(2+B—a)—}[1—I(v,)]e” 
and o, = v/D, D being the molecular diffusivity into air of the vapour 


concerned, replacing v/« the Prandtl number. The x distribution may 
then be obtained as in section 3. We note at this stage that the value of 
y, and v, are first obtained from equation (20) and the measured velocity 
profile, at a point over the downstream edge of the evaporation plate, and 
for other values of x, v,(x) must be chosen numerically so that the expression 


(x1—Xo)/(x1— Xo) [1—1(v,)} 01 — (x1 — x0) (xo— xX? (46) 
is independent of x, this being a necessary condition for the validity of the 
analysis. Takinga, = 20cm.,2, = 50em.,and U, = 676cm./sec. (therelevant 
values in the evaporation tests), Table 2 shows that the ensuing approxi- 
mation to the top of the transition zone is in close agreement (apart from 
a small error near the leading edge of the evaporation area) with the 
position of the top of the transition layer, calculated from equation (20). 


TABLE 2 
Calculated height of the laminar sub-layer, the transition layer, and the fully 
turbulent layer over an evaporation plate, 30 cm. in length, positioned with its 
leading edge 20 cm. downstream from the origin of the momentum boundary 


layer: U, 676 cm./sec. 
Downstream distance from leading edge of plate (cm.): 0 10 20 30 
Laminar sub-layer height (cm.), calculated from (18) . 0-02 0-02 0-02 0-02 
Transition layer height (cm.) ca ulated from (20) ; 0-12 0-13 0-13 0-14 
\pproximate transition layer height calculated from (46) 0-10 0-12 0-13 0-14 
Approximate height (cm.) of the fully turbulent zone 
educed from the experimental results of Davies and 
2 2-0 


Walters ; . | 08 1-2 1-6 


A convenient expression for the rate of evaporation from an area of 
length / may now be obtained in two ways. We first choose a value of k 
so that the measured mean velocity profile is in agreement with Townsend’s 
velocity profile which may be represented by the power law 

U/U, = 1-72E08, 
The expression for the downstream flux of vapour in the fully turbulent 
li 
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zone at a given station x = 2, = a,+/ may then be written in the form 


x x 


My , TOT ly . one \—0°15 i __ A\9,0°15 
Y (21) (x—Xo) dy l 72U(x, Xo)(A2,) ‘ (1—@)y*® dy, 
300/U; 300/U; 

(47) 
where 6 is the vapour concentration distribution corresponding to equation 
(30). The flux in the transition and laminar sub-layers is given by the sum 
of expressions (32) and (33) with pc, 7' replaced by x, NV by M, and with 
p=g¢@ 2 

An alternative expression may be obtained, as in section 3, by inte- 
grating expression (40), assuming the sub-layers are absent, and using 
equation (26) of Part I. The total flux from an area of length /, when its 
leading edge is at x = a,, is then given by 


Q’()) = d{(a,+1)’$—aj]}"*, (48) 


where 
d (xi —Xo)7 [1 —L(v,)]- > 


x [1 +5MKU,(kU,)-Ylog,(50, + 1)+0,}(1 —ajz/ayz) shite zs S. 


T7.b\Ha x1+8)~! se Soa : 
x | 0 3} [ak U,(2 } B- x) Hk pia+BX1+B)~*] (.)SIn(jz77). 


1+f 
On the same basis an expression for the flux in the region from y, = 0 to 
y, = 30 is given by (35) with xy and M replacing pc, T and V,andp = q = 1. 


Following the method of section 3, corrections for the presence of the sub- 
layers are then applied and the total flux across 2 = x, may be calculated. 

In order to evaluate these flux-expressions numerically appropriate 
values of a,, k, K, ete., are first deduced from the results of experimental 
work by Davies and Walters (2). Since the velocity profile was measured 
at one value of x only, the exact position of the virtual origin of the 
momentum boundary layer is not known, and consequently an exact deter- 
mination of a, is not attainable. However, an approximate value may be 
deduced, and we find that the calculated values of rate of evaporation EF are 
not very sensitive to our numerical choice of a,—at least in the conditions 
investigated. The experimental results obtained by Davies and Walters (2) 
show that the dependence of F on length / may be represented as a power 
law proportional to 1°76 for plates of lengths between 8 and 30 em., and 
we now choose the value of the parameter a, so that a power law 18 may be 
fitted to the values of the expression Q’ = d[{(a,+1)*—aj]"*, given by (48), 
for values of / between 8 and 30 cm. We find that for a, < 10 cm., 


8 > 0-80 and for a, 45 cm., 8 < 0-70 and for 17 em. < a, < 23 ecm., 
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3 is approximately 0-76. We first take a, = 20 em. and later check the 
sensitivity of the calculated evaporation values on this choice, by taking 
Ly 15 em. and a, {5 cm. This application of the Q’ expression is 
justified by the result that its computed value is found to be within 2 per 
cent. of the value calculated later in this section for the total rate of 


evaporation. 























] 
= 
Hf 
co Ps — 
91 r 02 -03 
Fic. 2. Mean flow distribution U/Uj, measured by 
Davies and Walters at U, 676cm./sec. € y/(1-32), 
wher 85 em. and is measured from the estimated 
position of the leading edge of the momentum boundary 
lavel U/U, L-72e°%. 


The parameter k is now determined by considering the wind velocity 
profile measured at a point over the downstream edge of the flat surface 
of the faired holder used by Davies and Walters (2); these velocity measure- 
ments are shown in Fig. 2. The downstream distance of this velocity profile 
from the leading edge of the evaporation area was 65 cm. so that its distance 
from the virtual origin of the momentum boundary layer is found (by 
idding the a, value) to be about 85 em. Values of the similarity variable € 
it this downstream distance are then given by y/(85k), and by choosing 

|-3 we find that, except for the values very near the top of the turbu- 
ent boundary layer, the measured mean velocities are reasonably close to 
Townsend’s profile: the departures near the top of the boundary layer 
ie probably due to a different degree of turbulence in the free streams 
fthe two experiments. It is suggested that the value of & is greater than 
|-0 (the value found by Townsend) due to some thickening of the turbulent 


ayer induced by the slight curvature near the leading edge of the faired 
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holder used in the experiments. Since the measured thickness of the 
turbulent boundary layer at this downstream position was about 3 em., 
the approximate value of € corresponding to the top of the boundary layer 
is now seen to be given by 3/(1-3 x 85) = 0-028 which is roughly the value 
of € measured by Townsend (4) in this region, and, since the conditions of 
flow were similar, this approximate agreement provides further support 
for the numerical value of a, assumed in this calculation. We also note 
that, since the measured distance L from the tunnel entrance to the position 
of the leading edge of the evaporating plate was 190 cm. and a, is taken 
to be 20 cm., the boundary layer on the tunnel wall becomes fully turbulent 
at a Reynolds number given by (L—a,)U,/v = 170 x 676/0-15 = 7-6 x 105 
which is within the range of possible Reynolds numbers given by Goldstein 
(9, p. 326). The velocity profile now coincides with the profile given by 
Townsend and we assume, following the arguments given by Davies (5), 
that K is 1-3 and the values of h(€) given by Townsend are applicable. 
We work, as in section 3, with two power law representations between 
which lie these values of h(€); we take 


(i) h = 0-00065€°# and (ii) A = 0-00088E%%*, 


Appropriate numerical values of several other parameters required in the 
calculation of rate of evaporation are now easily obtained. In the conditions 
of the evaporation experiments by Davies and Walters (2), x is zero and 
x, is calculated from a knowledge of the molecular weight of the liquid 
concerned, the temperature of the air stream, and a correction applied to 
obtain the surface temperature of the plate (this differs slightly from the 
air stream temperature due to the latent heat effect of the evaporating 
liquid). The value of D, the vapour diffusivity of aniline at the temperature 
(17° C.) of the plate, is also required and is given by Mack (14) to be 0-069, 
and v is 0-15 em. sec. units for air; this leads to a value for o, (= v/D) of 217. 

The expression (47), giving the flux of vapour and the associated sub-layer 
correction factors were first computed by numerical integration for 
a, = 20cm.,/ = 30cm. (the maximum length tested), U, = 676 cm./sec. (the 
measured free stream velocity), using firstly the power law h = 0-00065£°. 
The total calculated flux by this method is found to be 0-0470 g./min. fora 
width of 10-1 cm. (or 4 in.), the contributions from the fully turbulent zone, 
the transition zone, and the laminar sub-layer being 0-0326 g./min., 0-0134 
g./min., and 0-0010 g./min., respectively. This is in reasonable agreement 
with the experimental value for the total flux, given by Davies and Walters 
(2), of 0-0420 g./min. Calculations of evaporation from areas of width 4 in. 
are considered, as the dependence of evaporation on length of area was 
evaluated experimentally in detail by Davies and Walters (2) for this 
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particular width only. We assume in this paper that 4 in. is sufficiently 
large in the conditions of the experiment, and over the length range con- 
sidered, to permit us to neglect any effect which the enhanced rate of 
evaporation in the neighbourhood of the longitudinal edges may have on 
the index 8. In support of this assumption we note that measurements 
were also made by Davies and Walters to determine the dependence of EZ 
on width w of area, for lengths of area / = 5 in., 13 in., and 15 in., and 
from lines of regression obtained in each case for EF against w, we find 
that for w = 1-5 in., 2 in., 3 in., and 4 in. the associated 8 values are 0-71, 
0-745, 0-755, and 0°76, which suggests that 8 = 0-76 is close to the value 
of the index appropriate to evaporation from an area of very large width. 

The second method of calculation using expression (48), with the appro- 
priate correction factors, and the same power law h = 0-00065£°*° gives 
a calculated value of 0-0464 g./min., which is again in reasonable agreement 
with the experimental value and also in close agreement with the calculated 
results of the first method. This second method of calculation has the 
advantage of not requiring a numerical integration, and, due to our choice 
of a,, the dependence of theoretical values of E on length of area is in 
agreement with the results of experiment. 

Using the power law representation h = 0-00088£°4® and expression 
(48), for convenience, with the associated corrections, we now find that 
the total rate of evaporation in the same conditions, i.e. a, = 20 cm., 
1 = 30 em., and L, = 676 em./sec., is 0-0444 g./min., compared with the 
result 0-0464 g./min. obtained on the basis of the power law representation 
h = 0-00065€°4°. We note that these are theoretical values associated with 
extreme power laws enclosing the actual ¢ distribution given by Townsend; 
both are within reasonable agreement of the experimental value of 0-0420 
g.min., and similar results are obtained for other values of /. Finally, in 
order to test the sensitivity of the dependence of the theoretical results on 


the value of a, we find by the second method of calculation (i.e. using 


1 
expression (48)) at 1] = 30 em. and U, = 676 em./sec. that EH = 0-0460 
g./min. and 0-0491 g./min. for a, 15 em. and 45 cm. respectively, the 
former being within | per cent. and the latter within 6 per cent. of the 
result computed for a, 20 em. 

The results obtained in this section therefore suggest that, although an 
exact calculation was not attainable due to the absence of an experimentally 
determined a, value, no large error has been introduced into the calculation 
by (a) the approximate estimation of the parameter a,, and (b) the existence 
of a small range of possible a and « values, giving power law representations 
of the eddy viscosity distribution. The results show that the two-dimen- 
sional diffusion model described in this section gives results in reasonable 
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agreement with experiment in the conditions discussed. However, the 
problem of accounting for lateral diffusion and finite width of area in these 
conditions of a developing momentum boundary layer remains to be solved. 


5. Forced convection into a fully turbulent boundary layer which 

is preceded by an appreciable laminar boundary layer 

Since the method of sources is found to lead to a relatively simple 
approximate method of approach to the heat transfer problem in both 
laminar and fully turbulent boundary layers, it is of interest to consider 
its possible application to the important problem of flow over a flat surface 
for Reynolds numbers at which there is a significant laminar zone of flow 
near the leading edge followed by a transition zone with a final fully turbu- 
lent region of flow. We have not been able to treat equation (1) and conse- 
quently we have not obtained a solution for the transition zone. However, 
it is possible that an approximate treatment could be constructed by ignoring 
the transition zone and assuming the existence of a laminar zone up to the 
virtual origin of the fully turbulent layer at x = 2,, say, calculated from 
the observed profiles at various points in a fully turbulent zone. If 7, is 
the constant temperature of the whole plate and (as in section 3) 7", the 
surface temperature of the turbulent layer when the sub-layers are assumed 
absent, then the temperature excess at a point P on y = 0 for 2 > a, due 


to the known distribution of sources in the laminar zone up to 2, is 


X, 
(7, —T))7-1(sin p, 7) | (X—0)-#: 61-1 dd = f(X), (49) 
0 
where 
X , X, a, Sy (2+-8)(2+26)-1, fy (1+-8)(2+)-, 


and f(X) is obtained by a numerical integration of (49). The temperature 
excess at P due to an unknown distribution of sources along x > 2, is given 
(see Part I, section 2), by 

sz} A(d)d4 $2)/82( YY’ d) He dd. (50) 
where 


X’ = (x—2,)", po = (1+8)(2+B8—a)-, and s, = q(2- 


QO 
” 

P< 
LO 


and the integral equation to give the unknown function 
A(¢) 8 PH -Vis24(h) 


X’ 


is f(X') + | (X’'—¢)-Hu(d) dd = T,—7). (51) 


0 
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Its solution is 


, d 
uld) 7 SIN [Ly 77 
bi dd 


. 
0 


(6—X’e-7"—T.—f (X’)] dX’. (52) 


As in section 3 two further equations could be obtained by using an approxi- 
mate treatment of the sub-layers, such that 7; is a constant and expressed 
in terms of 7}. 

The numerical work involved, however, would be extremely heavy and, 
as the scatter in the experimental measurements by Elias does not allow 
a reasonably accurate determination of 2,, this has not been attempted in 


this pape l 


6. Conclusions 

A solution of the problem of forced convection from a flat plate through 
a fully turbulent boundary layer has been obtained by transforming the 
temperature equation into a von Mises form and using the method of sources 
to determine the distribution of rate of heat transfer on the plate and the 
associated temperature distribution over the plate. The treatment is based 
on (i) an application to flow over a flat plate of a power law representation 
of the eddy viscosity distribution, which was obtained by Townsend from 
his measurements of mean Reynolds shearing stress and mean velocity 
profiles over the plane wall of a wind-tunnel, (ii) the equivalence of momen- 
tum and heat eddy diffusivities, and (iii) an approximate treatment of the 
flow in the transition and laminar sub-layers. The agreement of the 
calculated results with those obtained previously by another method and by 
experiment then supports the application of the source method to discuss 
the problem of diffusion of vapour from a saturated surface into a developing 
boundary layer. The discussion of this problem includes the case when 
the leading edge of the vapour boundary layer is downstream of the leading 
edge of the momentum boundary layer, and incorporates an approximate 
treatment of the downstream flux of vapour in the sub-layers. The calcu- 
lated values of evaporation rates are found to be in reasonable agreement 


with measured values in a wind-tunnel. 
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SUMMARY 

This paper is concerned with problems of thin elastic plates with simply-supported 
boundaries. Solutions are obtained for plates in the form of half-planes and strips, 
including cases of isolated interior loading. A uniqueness theorem is established to 
ver both cases, 
1. Fundamental equations 
THROUGHOUT this paper the mid-planes of the plates are chosen as the 
plane Z 0 of a rectangular cartesian frame O(x, y, Z), and the notation 
used is that of A. C. Stevenson (1). The displacement w of the mid-plane 
in the direction k = OZ satisfies the equation 

Viw F(x, y), (1.1) 

where F(x, y) is zero if the faces of the plates are stress-free. Corresponding 


toa right-handed orthogonal triad fi, §, k the mean stresses nZ, sZ and the 


mean stress couples ns, sn, 88 ( nn) are given by the equations 
nZ-+-isZ 2Pe- _ V2 x (1.2) 
Cz 

ns sn =_— P( 1+ n)Vi Ww, (1.3) 

Cw 
ns—- SN 27188 —4P(1—n)e 2i8 =e (1.4) 

cZ* 
where z = a+iy, fi = (cos@,sin 8,0), cos@ = i. ft, 7 is Poisson’s ratio, and 


P is a constant depending on the thickness 7' of the plate. The general 
solution to (1.1) may be written, in the case F(z, y) 0, as 


w = Q)(z)+-202(Z)+-w(z)+- (2), (1.5) 


where (2(z), w(z) are functions of z which are regular in the region R occupied 
by the plate. From (1.3)—(1.5) 


ns/P 2(1+ y){Q'(z) + 2'(Z)} —(1— yn) [eP{2" (Z) +- @"(Z)} + 


j 


”" 


+ e?!420"(z)+-w"(z)}]. (1.6) 


lhis paper is concerned with plates in the form of half-planes and infinite 


(Quart. Journ. Mech. and Applied Math., Vol. IX, Pt. 4 (1956)] 
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strips, so that over the boundaries, chosen as y = 0, y =a the angk 
6 = +47 and vs —#x. From (1.6), 


ge/P = 2(1-+ y){Q’(z) + Q'(2)} —(1— n)f2zQ" (2) +20" (z)+-w"(z) +6" (2)}. 


¢ (ln) ew ge, (I—n) (2, ay 1a) Bw 
ange Sp: : Pe aes BS 2 alee ) '(z)+- "(2 es 
"4 Uo oP a CY ? edie, 


2. Boundary conditions. Uniqueness of solutions 

Problems of clamped elastic plates have been considered by many writers 
In particular, Muskhelishvili (2) has given an extensive bibliography 
Fewer attempts have been made to satisfy the conditions at free or simply- 
supported boundaries, e.g. Timoshenko (3). The problems considered her 
involve specified values of w and 7s along the boundaries. For half-plane 
and strips this requires the specification 


w = f(z), yx = g(x) (2.1 


along y = 0, with corresponding equations along y = a in the case of the 
strips. If these boundaries are simply supported or hinged, then 


f(x) = g(x) = 0, ete. 


From (1.8) it is clear that the boundary conditions (2.1) are equivalent to 
specification of ois 
Oo-w 


w= f(x), = = h(x). (2.2 


0202 
To establish uniqueness of solutions let w,, w,, with corresponding pais 
of potentials Q,, w,; Q,, w, be distinct solutions to (1.1), having the same 
singularities in the closed region R defined by the plate and satisfying the 
same boundary conditions (2.2), ete. Also, let 
&2, Q,—Q, and ay W1— Ws. 
Then x = w,—vw, satisfies V} y = 0 in R and on the boundaries 
a= Sig =O, 

Since # extends to infinity it is essential to postulate orders of magnitude 
at infinity. It will be assumed that the demands made at infinity do not 
involve terms of the orders 

Q = o(1), w = O(1), Q’,w’ = o(2-), Q” = o(2z-"), (23 
where accents denote derivatives, but that all larger terms are specified. 
Comparison with the corresponding theorems for clamped plates (4) indi- 
sates that more information at infinity is required here in order to specify 
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,unique problem. From (2.3) it follows that, at infinity, 
Q jue | h 


52 23 
av Czy C ‘ 
Xx o(z). X oil). x 0(z a), x o (2-?). (2.4) 
. Oz C202 0202° 
Let C denote the real axis from x -Ltox = Land asemicircle of radius 
L above the real axis or the segments from x L to x L of the lines 


;= Oand y = a, together with straight joins of corresponding extremities. 
In both cases the integral J defined by 


T= | |%—Xaz+y CX di (2.5) 


| C2 CZ02 Ox 


tends to zero as L tends to infinity, since y, V? x are zero on the straight 
sections, the integrands are 0 (z~!) at infinity, and the lengths of C are O(z) 
ut most. Assuming continuity of Q,, Q5, Q5, we, wo in R, (2.5) may be 


converted to the surface integral 
° ; Fa a 2 y > 
I 21 x | ds. (2.6) 
p Cc 
FS 
Thus V?y vanishes throughout R and, remembering the conditions at 
infinity, it is easily seen that Q, = 0. The problem has now been reduced 


to one in which y = w,(z)+,(2), i.e. V2 = 0 throughout R, x = O(1) 
and @y/éz = o(z~*) at infinity, and y vanishes on the straight boundaries. 


ie Atal ; ; i 
The integra 7 { [v9(z) + a9(Z)]o4(z) dz (2.7) 


vanishes when L tends to infinity, and is equal to the surface integral 


. g , =? - y . ? 
J 20 | [ wo (z)@_(Z) | ds. (2.8) 
R 
Thus w)(z) = 0 throughout RF and from the conditions at infinity and along 
the real axis w, = ik, where k is a real constant, and hence w, = wW,. 


3. The simply-supported half-plane 
From (2.2), (1.5) it is evident that if w and yx are specified along the real 
axis, then the complex potentials must satisfy the conditions 


2 refxQ(x)+ w(x)! Fie), 2 re Q’ (x) h(x). (3.1) 
These equations are satisfied by 
| 2... i fide 4, 
42 (2) on ; ; 20)(z) T wz) on ee 9 (3.2) 


which are analytic in y > 0 provided that f(t), h(t) are sufficiently small at 


infinity (5) 
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4. Isolated loads in a simply-supported half-plane 
The potentials 


Q, = A(z—2p)iog{(z—z)/(z—Z)}, Wo(2) ~Z9 Q(z), 
A = —W/16nTP, (4.1) 


correspond to an isolated load W at z = zy = a+iB (8 > 0) but are other- 
wise free from singularities in the upper half-plane. From (1.5), (1.8) the 
corresponding displacements are 


Wo : A(z—2)(7—2Z,)log{[(z Z%9)(2— 29) ]/[(z—2)(Z—Z»)]}, (4.2) 





and 
G(x,y) = yey i. (1 7) ~ — Hog @ F020), % 70 0" “Zo. 
4} 4 dx" \ ~~ (z—Z)(2—z)  z—Z% Z—z,| 
(4.3) 
When y = 0, 
(We)y-0 =9, Ga, 0) = 4AB2/{(e—z9)(e—%)}. (4.4) 
Thus the boundary y = 0 will be simply supported if potentials Q,, «,, 
obtained from (3.2) with f(x) = 0 and h(x) = —G(a,0) are added. From 
(3.1), (3.2) . 
, } q 2/ 2 t =_ DQ. 
ep Ec te ee (4.5) 
a J (t—2)(t—Z)(t—z) z— 
and hence 
Q(z) — 2iAB log(z—Z,), w,(z) —2zQ,(z). (4.6) 


Combining (4.1) and (4.6), potentials giving the isolated load and simple 
support along the real axis are 
Q(z) = A(z - ,)log{(z—Zz)/(z—Z,)}— A(z)—Z,)log(z—z,), (4.7) 
w(z) = —AZ,(z- z,)log{(z—z9)/(z—Z,)}+ Az(z,—Z,)log(z—zZ,). (4.8) 
At infinity the above solution is of the form 
Q(z) —A(z)—2Z,)log z+ A (Z)—2 )+0(1) | (4.9) 
w(z) = A(z,—2,)z log z+ O(1) J 
For the uniqueness theorem to apply in this case it must be supposed that 
the terms mentioned explicitly above at infinity have been specified along 
with the load W at z) and the boundary conditions along the real axis. 
The reason for this is that equilibrium is impossible without stresses or 
couples at infinity to balance the moment of W at z, about the real axis. 
The above terms are supplying the equilibrating forces. To this extent the 
above solution is non-physical. A physically satisfactory solution may be 
obtained by combining loads W; at z; (i = 1,..., n) such that 
> W(z;—zZ,) = 0. (4.10) 
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The solution in this case is 


Q(z) = > A,(z—z,)log{(z—z,)/(z—Z,)}— > A ,(z;—2,)log(z—zZ,), (4.11) 
w(2) > A: 2;(z—2,)log{(z—z,)/(2—%,)}+-2 > A,(z,—2,)log(z—z,). 

(4.12) 

At infinity Q(z o(1), w(z) O(1) (4.13) 


and the solution is unique 


5, Alternative treatment of the half-plane 

The above simple method of solving half-plane problems is not the most 
convenient for extension to consideration of infinite strips. An alternative 
method is therefore given. From (1.5), (1.7) it follows that 


w 20)(z) + 202(Z)+- w(z)+ @(2Z), (5.1) 


2 ox 2r 
Problems involving specified values of w and yx along lines parallel to the 
real axis may be solved by obtaining the corresponding values of v and w. 
Consider potentials Q,, w, related by the equation 
w,(z) 202, (2). (5.3) 
fv,, w, are the corresponding functions given by (5.1), (5.2), then 
(?,) a0 4re Q)(2x), (w,),a9 = 9. (5.4) 
Similarly, if potentials Q,, w, are related by the equation 
w,(Z) zQ,(z)+2 | Q,(z) dz, (5.5) 
then the corresponding functions v,, w, satisfy the equations 
(v,),,.9 = 9, (w),-9 = 4re | Q,(x) dx. (5.6) 
As in the previous method the problem has been reduced to the determina- 
tion of functions which are analytic in the upper half-plane and have 
specified real parts along the real axis. The boundary conditions 





= p(x), yg = g(a) (5.7) 
may be obtained by potentials satisfying (5.3), (5.5) provided that 
é l ; iul y F p(t) dt xo 
$02, (Z) : eu du pity dt “ | 7 : (5.8) 
4 {2,(z) dz | e' du | q(t)e iut dt a | We, (5.9) 


may be expressed as Fourier integrals. 
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6. The infinite strip 


Extension of the above theory to include a second boundary requires 
potentials which do not disturb the conditions obtained along the real axis, 
From (5.3), (5.4) it can be seen that the potentials Q,, ws, with corresponding 
functions v3, wz, where 


w,(z) = —2zQ,(z), 02,(z) - | {(yy+e 


~~ 


/, jet=u4 (y, iyo)e isu} du 
(6.1) 


give (¥3),-9 = (W3),-9 = © for all real functions y,, y, of the real variable u 


(6). A similar statement holds for the functions Q,, w,, v4, wy, where 


w,(z) = —20,(z)+2 [ Q2,(z) dz 
, , oe ee (6.2) 
Q,(z) dz = - | (5, +18,)e"— (8, —18,)e-#) du 
, | 


Thus, after obtaining specified conditions along the real axis by either of 
the methods given in §§ 3 and 5, the above pairs of potentials may be 
added and the functions y,, y2, 5,, 5, determined to give specified con- 
ditions along a second boundary. This section is concerned with specified 
values of v, w along the line y = a. From (6.1), (6.2) 


s ‘ : 
(Ug +4)y=a {| —acu*(y,+-ud,)—y,us]cos ru+ 
+[acu?(y.+ud,)+-y, us|sin xu} du (6.3) 


(Ws + W4)y, a 


76 
ag 


{[ac(y,+ud,)—8, s]cos xu 





+ [—ae(y.+ud,)+6, s]sin xu} du 


where c coshau and s sinhauw. It will be assumed that the conditions 
along y = a can be expressed in the form 


L 
. 


5 . 
Vy=a | (7, COS xu-+-7,sin xu) du 


0 : 
° (6.4) 
x 
ee ; 
a (a, cos xu+-oysin xu) du 





. 


The above functions are chosen to eliminate the unwanted terms obtained 


along y = a when getting the required conditions along the real axis as 
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vell as establishing specified conditions along the second boundary. Com- 


parison of (6.3) and (6.4) gives the equations 


acu*(y,+ud,)—y, us 7, | 
ACU"(Yo+Uds)+yot $ T2 . 
: z (6.5) 
ac(y;- uo,) 0,s O71 
aAC(Ys US») ; 52S = Oy 
Hence 
y, us T,(acu—s)+acuo, 
0; USs* acT,—(acu-+s)uo, a 
(6.6) 
Vo US" T,(acu— s)—acu*o, 
0. US* ACT,+-(ACU-+- S)Ucy 


From (6.6) it is evident that the differentiations employed in deriving 
6.3) require restrictions on 7,, T., 0), 2 at infinity. It will be sufficient if, 
at infinity, 7,, 7. o(u-3) and o,, o, = o(u->). Further, if 71, 72, 01, 02 are 
bounded at u 0, then y,, y. are also bounded, but 8,, 5, are O(w-*) and 
modifications are necessary before the solution is valid. However, the 
formal pair of potentials 


LF oe 67 
Q.(z) dz (8 10 l+-azu—z*u*/2) (6.7) 
(5, —18,)(l—izu—z?u?/2)} du 


give v. = w 0 everywhere and may be subtracted from (6.2) to give 
satisfactory integrals without changing the boundary conditions or making 
changes in the conditions at infinity. Indicating the modified potentials 


by a star. (6.2) becomes 


Ww; ) OF ) Z QF (z) d 
; lf ' 6.8 
() / s(d, r ] zu+2°u-/2) (6.3) 
2 
0 
«.* 
(6, —18,)(e-&"— 1+ izu+27u?/2)} du 


The solution to the strip with specified w and jx along y = 0, y = a is 


herefore given by the potentials (6.1) and (6.8), T1, Ta, 84, Fo being deter- 
mined from the boundary conditions (6.4) and "i> Ye 


5,, 5, obtained from 


6.6), 
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7. Extension to cases in which boundary functions are not 
expressible as Fourier integrals 
In many cases of interest the boundary functions are not expressible as 
Fourier integrals. In this section a method is given for dealing with an 
infinite strip under the same conditions as above except that 
x 
r(x) dx exists (7.1) 


—@ 


Wy=a r(x), where 





but 7(2) is not expressible as a Fourier integral. This case arises in the next 
section when dealing with an isolated load in the interior of an infinite strip 
with simply-supported boundaries. It is well known (7) that r(z) may be 
expressed in the form 


| - Lfdu f 0: 
r(~)—7(0) | r’(t) dé = —- | i | r’(t){sin u(a—t)+-sin ut} dt. (7.2) 
7 u 
0 0 ae 
sf , 
Thus r(x) r(O) + | {p,(cos ru—1)-+ pysin xu} au (7.3) 
7 u 


0 
where p,, pz are effectively the sine and cosine transforms of r’(a) respec- 
tively and are assumed to be finite but non-zero at u = 0. Writing p, = 0,4 
and p, = o,u, equation (7.3) becomes 


= x 
{> r(()) — 5 ] | 5 sf IOS 1797 LL 73 2 id " 
r(x) r(O) — o,du+ 10, COS TU-+-o, SIN ru} du, (7.4) 


0 0 
where a, o, are O(u-!) at wu = 0. It is clear that, owing to the form of 
(7.4), direct comparison between boundary conditions as in (6.3), (6.4) is 
not possible. Consider the formal pair of potentials Q,, w., where 


w¢(Z) 2 Q,(z) 
Q4(z) = ir(0)/(4a)—(2i/z) | {ey,-+ (aeu—s)8,/a} de (7.5) 

Clearly, 0 
~ ty\r(0) a I hak | [cy,+(acu—s)d,/a] dul. (7.6) 


0 


Thus the boundary conditions along y = 0 are unchanged by these poten- 
tials but 


x 


(We)y-a = 1(0)—(8/z) | {aey,+(aeu—s)d,' du (7.7) 
0 


r(0)—(8/7) | o, du (7.8) 


0 
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from (6.5). Thus these formal potentials give the first two terms in the 
expression for w,_, in (7.4) and the deduction of equations (6.5) is again 
possible. From (6.6) it is evident that 5,, 5, are O(u-*) at wu = 0, as before, 
but y,, ye are now O(u-") at u ) and further modification to ensure 
convergence at this point seems at a sight necessary. However, addition 
of Q5, ws in (6.1) and Qg,, we in (7.5) yields 





( 


we Ws We yd Qe 
O 01. 
a ——/- 
: 7.9 
ur(O)/(4a) (2/7) | typ tty )(et c)+ » ( ) 
0 
+- (vy —ty2)(e-*"— ¢)]—2(acu—s)5,/a} du } 


and it can be seen that no further modification is necessary, since the only 
divergent term in the formal potentials Q,, w, cancels the divergence at 
wu = 0 in the potentials Q,, w,. The complete solution in this case is given 
(5.8), (5.9) which yield the required con- 
OF, wt from (6.8), (7.9) which give specified 


by Q,; WW; hee. 5 


Wo in (5.; ».5), 
ditions along y = 0, and oF ot { 


conditions along y = a. 


8. Isolated load in the interior of a simply-supported infinite strip 
The potentials given in (4.7), (4.8) a to an isolated load at 

Z= 2 r.+iB (8 > 0) and give w = vs = 0, and therefore v = 0, along 

the real axis. Along y = a they give 

4A (a+)? 4A(a—B)* 


7 aan 8.1 
ee (a—a)?+(a+B one 





2 (w—a)?+ (a—B)’ 
x)?+ (a—B)* 





—4A Ba log|(x—«)?+(a+f)?}. (8.2) 


Thus, bearing in mind that the above boundary conditions have to be 
eliminated, the functions 7,, 7, defined in (6.4) are given by 





7, e% (7A COS at) T, €%/ (7A sin au) 8 cosh Bu+asinh Bu. (8.3) 
The function w,_, is not expressible as a Fourier integral but 
, | x 2 (a+ B)? 

yr’ (x) ; 1w 2A (x—a«)log \ 3 ( By (8.4) 
dx (x—a)*+(a—B)* 


can be so expressed. Using the method given in the preceding section, the 

functions o,, o, defined in (7.4) are, in the present case, given by 

o, we (7A cosau) = o, u%e™/(7A sinau) = Bucosh Bu—(1+au)sinh Bu. 
(8.5) 

Kk 
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Also 
9 9 x? (a T ‘ 9 9 > 
r(0) = Af{a?+ (a—B)?}log{ — By’ + 4AaB log{a?+ (a+ )?}. (8.6) 
x?-+- (a —B)? 
Since o,, a, are O(u-!) and 7,, 7, are O(1) at w = 0, and oj, oo, 71, 7. are 
O(e“8-®) at wu = o, the solution given by OF, wk, OF, wf is valid. 
Conclusion 


Comparison of the above with corresponding investigations on clamped 
elastic plates (4) shows the influence of boundary conditions on problems 
involving the biharmonic equation. For example, the conditions required 
at infinity in the case of the half-plane differ in the two cases, although 
these are sometimes dismissed by merely demanding ‘infinitesimal stresses 
at infinity’. The same approach has been employed here as in the paper 
quoted on clamped elastic plates, but it does not appear possible to apply 
directly the methods employed by Muskhelishvili (2) and Green (8) to con- 
sideration of simply-supported boundaries. 
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SUMMARY 


Timoshenko’s theory of transverse vibrations of beams predicts the existence of 


two possible types of wave motion of given frequency. The first type is progressive 
for all frequencies, whilst for the second type there is a critical frequency above which 
the waves are progressive and below which they are stationary. The reflections of 
an incident wave at a beam end are determined, and it is shown that at a simply 


supported end there is a single reflected wave, whilst at a free or clamped end 
reflected waves of both types are produced. On the basis of this theory a brief 


discussion is given of the propagation and reflection of a flexural pulse. 
1. Introduction 

BEFORE discussing Timoshenko’s theory it is useful to consider the reflection 
of waves according to the elementary equation for flexural vibrations of 


beams of density p, and mass m per unit length, viz. 


{ a2 
EI ‘ Yim J 0. (1) 
Cx? ot 


In this equation put 


where / is the radius of gyration of the beam section about the neutral 


axis, and ¢, /(E/p). Equation (1) now becomes 
} 2 
onu Cru 
0 a 
ve : (3) 
8s C7 


which has solutions in the form 


u exp| i(wr ps)|, (4) 
where w is assumed real and positive, and p may have any of the values 
Vw \w. In these solutions, which represent transverse waves of 


circular frequency Cy, w/k, the real values of p give progressive waves having 
phase-velocities « Cy \w, Whilst the imaginary values give stationary 
waves. Because of their exponential shape the latter waves cannot exist 
as free waves in an infinite beam. 


[Quart. Journ. Mech. and Applied Math., Vol. IX, Pt. 4 (1956)] 
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Consider a semi-infinite beam extending from s = —oo to s = 0, and let 
u = exp[i(wr—svw) | (5) 
be an incident wave at the end s = 0. Three types of end condition may 
be considered, namely, 


i Cu 

(a) Simply supported: u = 0, — = 0. 
os*~ 

29 23 

. e*u Ou 

(b) Free: = = 0, —— <= 0. 
ds? ds* 
‘ ou 

(c) Clamped: “u= 0, — = Q, 
os 


The end conditions (a) are satisfied by taking a reflected wave in the form 

wu’ = —expli(wr+svw)]. (6) 

In this case the incident wave is reflected with equal amplitude but in 

opposite phase, and the amplitude of the shear force at the end is twice 
that of the incident wave. 

With end conditions (5) and (c) it is necessary to introduce both a reflected 

progressive wave and a stationary wave, given respectively by the equations 


wu’ = —texp[i(wrtsvw)], (7) 


” 1 


u" = +(1—i)expliwr+svw], (8) 
in which the positive and negative signs refer to conditions (b) and (c) 
respectively. For these end conditions the incident wave is again reflected 
with equal amplitude, but with a phase change of 37. In the case of the 
clamped end the total amplitudes of both bending moment and shear force 
at the end are 2v2 times those of the incident wave. Again, because of their 
exponential form, the stationary waves will be comparatively insignificant 
except near the end of the beam. 
Green (1) has given a method of solution of flexural vibration problems 
in terms of the progressive and stationary waves considered above. 


2. Timoshenko’s equations 

The elementary bending theory gives an approximation which is satis- 
factory only for low-frequency waves in which w <1. A more accurate 
solution may be obtained by including Timoshenko’s (2), (3) corrections 
for shear deflexion and rotary inertia, which lead to the following equations 


of motion 


> ae (Uy, T Us) 0, (9) 


(10) 
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Using the notation of a previous paper (4), ku, and ku, are the deflexion 
components due to bending and shear deformation respectively, and y is 
a constant given by EI 


i Ch’ (11) 


where C is the effective shear rigidity of the beam section. The calculations 


in this paper are based on a value y 4. 








aes 
——— 
a ’ 
= T T 
5 10 Sw 
Fy 1. Variation of uw, and py with w. 


The deflexion components are related to the bending moment M and 
shear force F by the equations 
El cu, 


M F 
hk és? 


(12) 


PF=c™%. (13) 


cs 


: i ; ee O7tt, a6 1 Ou, 

[t is interesting to note that if the ‘coupling’ terms —” in (9) and 
OT? y cs 
in (10) are omitted, the two equations would each have the form of a simple 
wave equation, and would yield respectively a wave of pure shear deforma- 
tion with a velocity of propagation c,y~!, and a wave of pure bending 
deformation with a velocity cy. 
With the coupling terms present, solutions of (9) and (10) may be ob- 


tained in the form . 
' Uy, exp| i(w7 ps)| | 


, (14) 

u pwexp|t(wr ps)|} 
whe re pA y(p* w*) (15) 
and p may have any of the values +p,, +p., where p, and p, (p, Pe) 


are the positive roots of the equation 


w* p*)(w? p) —qw* 0. (16) 
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There are thus two types of wave motion possible, each having a dispersive 
propagation and consisting of a mixture of bending and shear deformation. 
The first type, corresponding to p /p;, iS a progressive wave since p, 
is real for all values of w. The associated value of ., is plotted in Fig. 1 for 
the particular case when y = 4, and it is seen that the wave consists mainly 
of bending deformation at low frequencies, and shear deformation at high 
frequencies. 
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Fic. 2. Phase and group velocities. 


For the second type of wave the value of p, is real or imaginary according 
yo? = 1. (17) 
This type of wave is therefore either progressive or stationary according as 
the circular frequency is greater or less than a critical value c)/kvy. Again, 
the wave consists mainly of bending deformation at low frequencies, but 
at high frequencies . approaches a limiting value —y/(y—1), and both 
bending and shear deformation are present. At low frequencies the two 
wave types approximate to the progressive and stationary waves predicted 
by the elementary theory. 

When » is real, equations (14) represent a wave having a phase-velocity ¢ 


and group-velocity c, given by the equations 


go 


c Fe >. (18) 
Co p 
Co ” dw (19) 
Co «dp 


The phase and group velocities of the two types of wave are shown in 


Fig. 2. 
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3. Wave reflections 

In order to determine the reflections of waves of the form (14) at the end 
s = 0 of a semi-infinite beam, it is necessary to express the end conditions 
in terms of the components u, and u,. By application of equations (12) and 


(13) the conditions may be obtained in the following forms 


1: C7Uy 
(a) Simply supported: w,+4, 0, — 0 
, os“ 
; Omu ou 
(h) Free: = 0. '.. §}, (20) 
cs” Cos 
, ou 
(c) Clamped: U,+u, = 9, ’—0 
; Cs 





The last condition is obtained by noting from (12) that 6u,/és is continuous. 
However, the clamped condition does not impose restriction on @u,/és, since 
(13) permits discontinuity in this quantity. 

At a simply supported end it may be shown that an incident wave of 
either type is reflected with the same amplitude but with a phase change 
of 7, as in the elementary theory. At a free or clamped end, however, a 
single incident wave gives rise to two reflected waves, one of each type. By 
putting p = p, in (14) we obtain an incident wave of type 1, and if yw? > 1 
the reflected waves may be written in the form 

Uy, = 4, eXp{i(wr+p, s)|+ a9, exp[i(wr+p,s)] (21) 


U., = fly Oy, EXpli(wr+p, 8)]+po ag, exp[i(wr+ pp 8)] 
Similarly, for an incident wave of type 2 the reflected waves may be written 
U) x12 EXp| i(wr- P,8)| | Xoo EXP| i(wr } P28)| ) 
U., = py Xo EXPli(wr+t p, 8)|]+ py ag exp[i(w7+p,8)] J 
The reflection coefficients a,,, a1, 2, % may be determined from the 
appropriate end conditions. Thus for a clamped end with an incident wave 
of type 1, conditions (20) (c) lead to the following equations : 
b+ py ayy (1 +p,)+a9;(1+ p19) 0, (23) 
Pit %y1 Pit % P2 = Y, (24) 
giving 
Pi + pe) + Poll + py) 
Py(1+-p2)—po(1 +p) 


__ ite) (26) 
Po(1+py)—py(1+ pH) 
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Similarly, it is found that 
a Sa (27) 


2po(1-+-49) 
Py(1+ pe)—po(1+ py) 
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Fic. 3.. Reflection coefficients for a clamped end. 
Similar results may also be obtained for a free end in the following forms: 


_ M2 Pita Po 


sea ? (29) 
hi P2—F2P1 
» 
a eee «|: (30) 
P2 2 P1— Pi P2 
e 2heDo 
412 = eS art. (31) 


Pr Ha P2—P2 Pr 
The numerical values of these coefficients are shown in Figs. 3 and 4. 

When yw? < 1, the waves of type 2 degenerate into stationary waves, 
and only reflections due to an incident wave of type 1 need be considered. 
At a free or clamped end there is now a reflected progressive wave and a 
stationary wave, and their forms may be determined by substituting 
ip, = q (where q is real and positive) in (21) and in the appropriate equations 
for a,, and a,,. It is seen that these coefficients are now complex, and also 
that |a,,| = 1. The incident wave is therefore reflected with equal ampli- 
tude but with a change of phase. 

When an incident wave of either type is reflected at a simply supported 
end, it may be shown that the total amplitude of the shear force at the 
end is twice that of the incident wave. At a clamped end the ratios of the 
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total bending moment and shear force to those of the incident wave may 
be written ¢,,; and ¢,, respectively for an incident wave of type 1, and 
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Fic. 4. Reflection coefficients for a free end. 
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Fic. 5. Magnification of bending moment and shear 
force at a clamped end. 


bm, Pr. for an incident wave of type 2. These ratios are given by the 


following equations 


de\* ‘i 
Dm P72 : 1+ X11 T xai(7 ‘ > (32) 
P 
bn = bua = 1—ay—O%g, 242. (33) 
My Py 


The numerical values of the ratios are plotted in Fig. 5. 
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4. The propagation and reflection of a pulse 

The propagation of a flexural pulse according to Timoshenko’s theory 
has been discussed by Davies and others (4-7). 

If the beam is given a transient disturbance (such as a transverse impact) 
at some point, the subsequent distribution along the beam of wave com- 
ponents of various frequencies will be determined by their group velocities, 
provided a sufficient interval has elapsed for the dispersion to become well 
established. It may be shown (4) that three wave-fronts are propagated 
in each direction from the point of disturbance, and these wave-fronts are 
associated with stationary values of the group velocity, i.e. 

—# = 0, (34) 

CW 
The fastest wave-front (A) has a velocity cy and contains high-frequency 
components of type 2. The second wave-front (B) has a velocity somewhat 
greater than cy y~', and contains components of type 1 having a frequency 
given by w = 0-354 when y = 4. This wave-front differs from the other 
two in that it only exists in an asymptotic sense, and is not apparent except 
at large values of 7 after the disturbance. The third wave-front (C) has a 
velocity cyy~? and contains high-frequency components of type 1. This 
wave-front is inverted, so that the disturbance associated with it lies ahead 
of the ‘front’ and not behind it. The first and third wave-fronts are of 
particular interest in problems of transverse impact, since the maximum 
stresses are likely to occur at these points. 

The ultimate response at any point is composed of low-frequency waves 
having a small group velocity, and this response may be determined with 
sufficient accuracy by the elementary beam theory. 

By considering the reflections of the individual wave components, it is 
possible to determine the nature of the reflected pulse due to the arrival 
of a flexural pulse at a beam end. When the end is simply supported it is 
easily seen that the reflected pulse is of the same form as the incident pulse, 
but of opposite sign. At a free or clamped end, however, the end conditions 
introduce an additional coupling effect between the two types of wave 
motion, and the reflected pulse is more complicated than the incident pulse. 
In particular it may be shown that the reflected pulse contains five wave- 
fronts. 

The first wave-front (A) of the incident pulse will bring high-frequency 
components of type 2, and the reflections set up at a free or clamped end will 
contain high-frequency components of both types. The first part of the 
reflected pulse will therefore contain two wave-fronts, of types (A) and (C), 
travelling at velocities cy and cy y~! respectively. The first waves of type | 
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arriving at the end will be those of frequency given by w = 0-354 (for 

4), and the reflections of these waves will produce a corresponding 
wave-front, of type (B), in the reflected pulse. Since progressive waves of 
type 2 do not exist at this frequency, the reflected pulse can contain only 
one wave-front of type (8). Finally, the arrival of the third wave-front (C) 
will bring high-frequency components of type 1, and, as in the case of the 
first wave-front, will result in a further two wave-fronts, of types (A) and 


(C), in the reflected pulse. 


5. Conclusion 

It should be noted that Timoshenko’s theory, on which the analysis and 
discussion in this paper is based, is only an approximation, although more 
accurate and physically more satisfactory than the elementary beam theory. 
Fundamentally the approximation lies in the one-dimensional character of 
the equations of motion (9) and (10), which give two types of wave motion, 
corresponding to the two roots of the frequency equation (16). In an exact 
two-dimensional solution of the problem (3), using the equations of elasti- 
city, a transcendental frequency equation is obtained, giving an infinite 
number of possible types of wave motion. 

It has been shown (8), (9) that the lowest phase velocity according to 
the exact solution agrees well, at all frequencies, with the velocity of waves 
of type 1 in Timoshenko’s theory, but there is no close correspondence 
between the velocity of waves of type 2 and any of the higher velocities 
predicted by the exact solution. It seems likely that Timoshenko’s theory 
will be satisfactory in problems where only the lowest velocity wave motion 
is important, but in other cases the accuracy of the theory requires further 


investigation. 
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THE CLOSED-FORM SUMMATION OF SOME COMMON 
FOURIER SERIES+ 


By CHI-CHANG CHAO? (New York, N.Y.) 
[Received 1 December 1955; revise received 27 March 1956] 


SUMMARY 
A method is presented for the summation in closed form of Fourier series whose 
coefficients are ratios of polynomials of certain types frequently encountered in 
practice. The result is obtained as the solution of a linear differential equation with 
constant coefficients, which can be solved by elementary methods. 


1. Introduction 

THE solutions of many problems in applied mechanics are most readily 

obtained in the form of Fourier series, and it is often not recognized that 

some of these may be summed in closed form, thus simplifying the results. 
Consider, for example, the problem of the simply supported vibrating 

beam shown in Fig. 1 whose steady solution is (1): 


_ Ag LE ot S ( sin(tza/L) | 


= op (1) 
Ayn 7 | (¢47r4a?—w?L4)i| 


which is a Fourier series whose coefficients are ratios of certain polynomials. 
A series of such nature is easily summed in closed form by the method 
described below. In this particular case the result is (see example at the 
end of the note) 


qu w , w , w ‘ w 
me , —2+ese /—Lsin /—(L—«x)+esch /—Lsinh /—(L—2z)+ 
2Ayw*| a a a a 
/ 
Ww , oe. w ‘ Ww ‘ 
ese -Lsin /—x+esch /-—Lsinh x|sinwt, (la) 
a a a 3 


It may be of interest to note that the summation of Fourier series in 
closed form has attracted considerable attention in a number of recent 
papers (2,3). Some simple series of the type considered in the present note 
may be treated by the methods of the papers listed in (2); the present 
method is advantageous in these cases since it employs only elementary 

+ The results presented here were obtained in the course of research sponsored by the 
Office of Naval Research, Contract Nonr—266(20). The author also wishes to express his 
appreciation to Professor B. A. Boley for his invaluable advice. 

{ Research Assistant, Institute of Flight Structures, Columbia University, N.Y. 


[Quart. Journ. Mech. and Applied Math., Vol. IX, Pt. 4 (1956)] 
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techniques. The method of (3) is applicable to a wider class of series than 
the one described below, but is more cumbersome in the present case since 
it requires quite complicated integrations. 

| y 


We sinat 


VTTTTTTTITTTT tii d 


==> X 























Fic. 1. 


2. Summation of sine series 


If a Fourier sine series is of the formT 


fe) = > iPM) ™), 2 
J\&X) Zz, Q(i2) in| L (2) 


where P and Q are polynomials, its summation can be easily performed as 


follows. Rewrite f(x) as 





? l m—1 
2(m—k)-1 j 2(m—k)-1 
o 2. + al *-*+(—IF 2 oa 
f(x) = 7 YY #0 —— —2=2__ gin a, & (2a) 
“ {=3 > a, af"-© 
a, A 
k=0 
9 - 
> J, sin X; 2, 
L ee 
t 1 
where x, = tn] L (2b) 


u 


and the coefficients a,, b,, and c, are known in any particular case. 
Without loss of generality one may take 


Ay - (2c) 


Note that the denominator of (2) is a polynomial of degree 2m, while each 
summation of the numerator is a polynomial of degree (2m—1) at most; 


thus convergence is assured. 


+ The summation procedure described in this paper is valid for the sine series of form 
(2) and the cosine series of form (8) only. Thus, for example, those series derivable from 
by this method. 


log cot @ are not manageable 
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From (2), a0 “fla )sin a, x dx. (3) 


With the aid of two iain by parts, one obtains 


I, = -—(1 a) f (x)eos a; &]0+-(1/a; )*[ f (x)sin a; a 
L 
—(1/a;)? | f"(x)sina;x dx, (3a) 
0 
where primes indicate differentiation with respect to x. Introduce now a 
function f,(7) such that 
f"(x) = a,f(x)—fi (2). (4) 
Substitution into (2a) gives, after two additional integrations by parts and 
some rearranging of terms 


(ao?-+a,)1; = —a,[ f(x)cos a; x] x; )[ f(x) cos a; ly 
[3 ‘)sina;2 dx. (4a) 


Terms containing (sin«;x) to be evaluated at x = 0, L vanish and have 
therefore been omitted. Introduce now a function f,(2) such that 


fi(x) = az f(x)—f,(2). (5) 
Repetition of the above operations gives, in general, 
fj-1(%) = a;f(x)—f,(x) (J By Bpn m). (5a) 


The final result is 


L 
m é zy, _ a 
["S fi, (x)o2?™-®—1¢e0s 2) + | fi»(x)sin a; x da 
ce =0 5 


I, = — a ee eet (6) 


Pp: a, afm—") 
The denominator of this expression is identical with that of (2); the numera- 
tors will also be the same provided that 

fn(x) = 0 (6a) 
and f,(0) = 6,, f,.(L) = —e,. (6b) 
Equation (5a) is the typical equation of a set of m linear simultaneous 
ordinary differential equations of second order with the functions 

a a? 

as dependent variables. Elimination of all functions f; from (5a) yields the 
following single equation of order 2m for f(x), namely 


m 
> (—1)"-*a, D*"-f(xz) = 0, (7) 
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where the operator D = d/dx. This equation may be solved under the 
boundary condition imposed by the 2m equation (6b). 

The form of (7) is virtually identical with the denominator of the Fourier 
series coefficient ((2.a) or (6)), in fact (7) may be immediately written by 
replacing a?"-" by (—1)"-*D%"-f(z) in the denominator of these 
equations. 

Thus, for a given Fourier series of form (2), the closed form summation 
of the series is obtained by the solution of an ordinary differential equation 
whose form (equation (7)) is immediately determined by the denominator 
of the Fourier series coefficients and with boundary conditions (equations 


(5a) and (6 b)) determined from the numerator. 


3. Summation of cosine series 
A Fourier cosine series of the form 


x 


<— P(7?) ima , 
f(x) — COS -constant, (8) 
— ()(7?) L 
where P and Q are again polynomials, may be similarly summed. Rewrite 
1(x%) as: 
1 ; m—1 ” k-1) 
2 c—1 - 2(m—k— 
KK 2 2 2» of™ +(—1) > cya 
f(a) mal \Y K=0 , 2 x; x 
LO" Lhe me i 
; >a" 
k=0 
K 2 a 
= S K. cos x; X, (8a) 
b' bike 
i=] 


in a notation similar to that of (2). Proceeding in a manner similar to that 


employed for the sine series, the equation analogous to (6) becomes, in this 


case, L 
n 4 . ™ : L - 
> fi(x)o?"—*— Neos a; x] 1 fin(X)Cos a; x dx 
k=0 A 
K — Rascsinesiieenienens (9) 
m 
s ay, yz(m k) (2 (0) 
k=0 
where the f,. functions are still defined by (5a) and a, 1. The denominator 


in this expression is again identical with that of (8); the numerator will also 
be the same provided that 

f(x) C = constant (9a) 
and f’ (0) b,., fi{L) Ce. (9b) 


The differential equation for the function f, analogous to (7), is 


})” Ka, Drx»Am WF (a) c. (10) 
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The constants of integration appearing in the solution of (10) are found 
with the aid of conditions (9b). The constant C is evaluated from the 
condition L 
[ f(x) dz = K,. (11) 
0 
4. Example 
As an example, consider the series 





Sin a; 2, (12) 


: a,[8,—(—1)'0,] 
f(x) = Z x8 ae 


i=1 


which would be identical with equation (1) except a common multiplier if 
0, = 0, = 3 


Then, in the notation of (2a), 

m= 3, a, = l, ay 

= b,=—¢, = ¢, = 0, 
The differential equation (7) is 
foro py” = 0 
and the boundary conditions are 
f(0) = f(L) = f"(0) = f"(L) = 0, 
fir(L) = 0 — f'*(0) = 

The final results are then 


9 
ap — 2+-ese BL sin B(L—2x)+esch BL sinh aL—)| + 


(12b) 


oy 
1. a,| — = + ese BL sin Bx+esch BL sinh Be 


With this result (1a) is immediately obtained by setting 6, = 0, = }. 
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